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ABSTRACT 


A semi-closed Rankine cycle vropulsion plant for the 
propulsion of a small deep diving submersible vehicle is 
analyzed. The propulsion plant is fueled by cryogenic hydrogen 
and oxygen. 


Weights and volumes, efficiencies of various plant com- 
ponents are estimated, and the plant is compared with competi- 
tive propulsion plants for a 50 kw, 1000 kwh mission at depths 
of 8000 and 20,000 feet. Weight of the plant in air is about 
5500 pounds heavier than a comparable Brayton cycle propulsion 
plant at the 20,000 foot depth and about twice the weight of a 
hydrogen-oxygen fuel cell plant. It is about half the weight 
of a comparable silver-zinc battery system. At the 8000 foot 
depth it is about 1000 pounds heavier than the Brayton cycle 
plant. 


A porous plug combustion chamber employing premixed hydro- 
gen, oxygen and steam is evaluated, and a novel three-four disk 
axial impulse re-entry turbine is designed using turbine inter- 
stage leakage equations developed. A computer analysis of film 
condensation with superheated steam in the presence of three 
non-condensable gases is presented for a horizontal tube and 
emell condenser. 
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CHAPTER I 


Introduction 


The need for new propulsion plants to power deep diving 
submersibles becomes apparent when one considers the power and 
endurance limitations on present battery powered vehicles. The 
associated heavy weight of batteries, even if they are exter- 
nally mounted, makes mission life short and limits available 
power levels to low values. 

The small submersible requires tending by a larger service 
vessel which must transport the submersible to the operating 
area, place it in the water, and remove it at the end of the 
mission. This launching and recovery process is severly ham- 
pered in any but the best weather conditions. While it has been 
suggested that such a vehicle operate independently and thus 
eliminate the launch and recovery problem, insufficient 
Emeurance levels seem to preclude this for anything but nuclear 
powered deep diving submersibles. To build a vehicle capable 
of extreme depth (8000-20,000 ft) operation and also capable of 
extended endurance would require a larger than normal submer- 
Sible crew (two or three) as well as a significantly larger 
portion of the vehicle capable of withstanding extreme pressures. 
In essence, this would require larger living space which, in 
turn, would place even more difficult requirements on pressure 
vessel design. The weight of a spherical vessel, the most 
efficient shape for extreme depths, increases in an exponential 
fashion as vessel diameter is increased. Thus, an indevendently 


Operating submarine would probably require a series of 


Wy 


—— be 






relatively small interconnected spheres for minimum weight. 

At present, the best solution to operatiig at extreme 
depths for missions of 24-36 hours duration appears to lie 
with the tended submersible. This choice has led to a requirement 
for small, high energy density propulsion plants. While the 
volume of such a plant is of concern, its overall weight in air 
is of primary importance, since it is this weight that the 
service vessel must lift. This vehicle weight, of course, 
includes sufficient flotation material to make the vehicle 
neutrally buoyant in water. Thus, in the comparison of various 
propulsion plants for the same mission endurance and the same 
full power level, one should include sufficient flotation 
material to make the propulsion plant neutrally buoyant. 

In this thesis a semi-closed, hydrogen-oxygen fueled, 
Rankine cycle propulsion plant is investigated as a power 
source for submersibles operating at extreme depths in the 
range of 8000-20,000 feet. The Bat antsacs of such a plant are 
the high energy density of the reactants and the virtual 
absence of non-condensable combustion procuete: 

The proposed propulsion plant is effectively an open 
steam cycle with recirculation of a portion of the condensate 
as feed water to dilute the high temperature combustion gases. 
A schematic diagram of the basic cycle is shown in Figure I-l. 

In order to determine the optimum operating conditions 
for the plant, thermodynamic studies were conducted. The results 
of these studies required analysis of the effect of these ovtimum 
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Analysis of such a propulsion plant requires detailed 
investigation of the most critical components in order to 
determine reliability, safety, efficiency, size and weight. 

In the course of the thesis, it was determined that these 
critical components were the combustion chamber, the turbine 
and the condenser. The combustion chamber was important from 
the aspect of reliability as well as safety and high combustion 
efficiency. The turbine was found to be of critical importance 
with respect to the overall plant thermal efficiency. The 
importance of the condenser lies in its size in relation to 
other plant components, and in the safe removal of unburned 
hydrogen and oxygen from the plant. 

Using information from the detailed analysis of the 
various plant components and their efficiencies, the overall 
plant configuration is developed and overall weights and 
volumes of the power plant and associated containment vessels 
are presented. The overall weights are corrected for sufficient 
flotation material to achieve neutral buoyancy. The propulsion 
plant is then compared with other existing and proposed propul- 
sion plants for comparable missions. 

In contrast to the normal thesis format, lists of 
references and symbols are found at the end of each respective 
chapter or appendix. This was necessary because of the large 
number of engineering disciplines encompassed by the thesis 


and the nomenclature peculiar to each. 
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CHAPTER =: 


Previous Work 


Chemical-dynamic power systems for deep submergence 
applications have been studied by several investigators. 
Reference 1, a study conducted by the Committee on Undersea 
Warfare, National Research Council of the National Academy of 
Sciences, compared various candidate propulsion plants including 
the fuel cell. This report recommended parallel investigation 
of fuel cells and chemical dynamic systems for a period of 
several years, followed by development of the most promising 
candidate. Sternlicht and Bjerklie (2) assessed a basic 
hydrogen and oxygen fueled open expansion cycle (discharging 
to ambient sea pressure) as well as other dynamic systems and 
determined that in the medium power, medium endurance range 
(about 50 kw, 1000 kwh), the closed dynamic cycle and fuel cell 
were comparable in an analysis of reliability, weight, service- 
ability, development time and development cost. 

Morrison, McCartney and Blose (3) assessed the depth 
insensitive heat engine using hydrogen and oxygen as reactants 
as well as other competitive propulsion plants, but presented 
no details on the specific cycles considered. For a steam 
boiler and turbine plant below 200 HP, they estimated a basic 
cycle efficiency of 20%. They also stated that for power 
levels less than 100 HP, the reciprocating steam engine appeared 
to be more efficient than a steam turbine. 


Catterson and Swain (4) assessed the semi-closed power 






cycle for operation with various fuels, among them hydrogen 

and oxygen, and compared them with H5-0, fuel cells, batteries 
and radioisotope-dynamic systems for a 30 kw power level. For 

a design depth of 6000 feet, the chemical-dynamic system was 
considerably lighter, including flotation material, than other 
systems in the range of 2-30 hours mission duration. They 
predicted an overall thermal efficiency of 20-30% for a Rankine 
cycle using water as the working fluid. These figures represent 
the results of an investigation of the semi-closed cycle of this 
thesis and the normal closed steam cycle with a boiler. Catterson 
and Swain recognized the weight penalty of encapsulation spheres 
for cryogenic hydrogen. 

Walter (5) proposed the cycle of this thesis using cryo- 
genic oxygen with cryogenic hydrogen or methane as fuel. The 
proposed power level was 400 HP. For a 1000°C (1832°F) turbine 
inlet temperature and a 29 lbf/in* condenser pressure, Walter 
predicted a 34% overall thermal efficiency for a plant burning 
stoichiometric methane and oxygen at 72.5 lbf/in*. The proposed 
turbine was a 5 stage machine. As in the cycle of this thesis, 
diluent fluid would reduce the combustion gas temperature to a 
Hevel acceptable for the turbine. The combustion products, 
water and carbon dioxide, would be stored on board. 

Balukjian and Rackley (6) investigated the closed Brayton 
Cycle propulsion plant for deep submerth ores for design depths 
of 8,000 and 20,000 feet for a 50 kw power level and at mission 
endurance levels of 1000 and 2000 kwh. They estimated a power 
cycle efficiency of 44.7% for a propulsion plant fueled by 
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cryogenic hydrogen and oxygen. This represents an overall 
thermal efficiency of about 33% (7). 


A supercritical CO. cycle has also been proposed for deep 


2 
submergence application by Karig (8) using a modified Feher 

cycle (9). Karig predicts an overall thermal efficiency of 

over 40% with this cycle. This cycle, essentially a semi- 

closed Brayton cycle, would operate in the same fashion as 

that proposed in this thesis, recirculating a diluent to reduce 
combustion gas temperatures to acceptable turbine inlet 
temperatures. 

Several other chemical-dynamic propulsion plants have been 
investigated and some have been built. The Alton Cycle (10) 
is one such plant. It burned diesel oil with hydrogen peroxide 
as the oxidant. 

The hydrogen-oxygen encapsulated fuel cell appears to be 
the most promising fuel cell candidate at present. It has not 
yet been installed on a deep submersible vehicle. The fuel 
cell plant has been discussed in many papers (11), (12). The 
estimated conversion efficiency of the hydrogen-oxygen fuel 
cell is in oe of 50% and may be as high as 70%. 

Several other fuel cell candidates exist. Ghormley and 
Harrison (13) discussed a hydrazine-hydrogen peroxide system, which 
would not require encapsulation. This system is now under 
development at the Naval Ship Research and Development Labora- 
tory. Overall weight for this system is reportedly about the 
Same aS the encapsulated hydrogen-oxygen fuel cell system. 
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It is apparent that the hydrogen-oxygen fueled semi-closed 
Rankine cycle may be expected to exhibit an overall thermal 
efficiency of 30% or perhaps greater. The degree of development 
of the major equipment associated with such an application is 
high. Only the combustion chamber, as noted by Morrison, 
McCartney and Blose (3), requires extensive development. While 
the achievable thermal efficiency is apparently far below that 
of fuel cells, it appears that the propulsion plant could be 
readily developed to provide reliable power for submersibles 
Operating at extreme depths. The cost of such development 
could be considerably cheaper than for a comparable fuel cell 
plant. It is the intent of this thesis to analyze the applica- 
tion of a hydrogen-oxygen semi-closed Rankine cycle for the 


propulsion of a submersible operating at extreme depths. 
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CHAPTER SLi 


Method of Analysis 


The design of a propulsion plant such as that proposed in 
this thesis must proceed on two paths. The first path is that 
of producing a feasible design. Each piece of equipment must 
be analyzed to determine if, in fact, it can be built to meet 
the constraints placed upon it, and the overall plant must be 
examined in the same way. The second path is that of producing 
an optimum design within the constraints of the mission. In 
this thesis, both feasibility and optimization have been 
considered. 

Safety and reliability of the plant are primary considera- 
tions which must govern the design of individual components as 
well as the plant as a whole. While absolute safety would be 
desirable, the probability of no accidents with hydrogen and 
oxygen cannot be guaranteed. The example of the space program 
in the safe handling of these potentially lethal reactants 
shows that it can be done if adequate precautions are taken. 

Reliability must be built into each component as well as 
the entire plant. The design goal for the semi-closed Rankine 
cycle propulsion plant is set at 4000-5000 hours time between 
overhauls. The common arrangement of a deep submergence 
propulsion plant is to house the plant and its reactants and 
products within several containment spheres (encapsulation) 
located within the hull form of the submersible. To open one 
or more of UiSeaene vessels because of an equipment 
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malfunction would require considerable effort and would not be 
normally accomplished in the field. Hence an adequate time 
between failures should be built into each component to meet 
the overall design goal. This goal may not be attainable. 

The analysis of a propulsion system requires certain 
known factors. Some of these factors are obvious from the 
character of the system itself. Others must be imposed by 
certain assumptions. For instance, it is readily apparent 
that the semi-closed Rankine cycle would not be competitive at 
shallow depths, since at these depths a plant exhausting com- 
bustion products to the sea is more easily built and less 
expensive. The use of noeeen as a fuel, moreover, incurs 
a great penalty in encapsulation weights over other fuels 
because it cannot be stored outside a pressure hull and because 
of its large specific volume, even as a cryogenic liquid. It 
is only at depths where the exhausting of combustion products 
to the sea requires excessive expenditure of pumping power 
‘that the depth-insensitive, combustion-product-condensing 
plant becomes attractive. This depth is less than 8000 feet. 

A depth of 8000 feet would permit operation in the area 
of the continental shelves and a large percentage of the 
continental slopes. Furthermore, it has been shown that a 
20,000 foot depth capability would permit access to about 98% 
of the ocean bottom. Thus two depths are of interest in this 


thesis, 8000 and 20,000 feet. 





A second factor, one which must be imposed, is mission 
duration and power level. The Navy's Deep Ocean Technology 
Program has established a basic 1000 kwh, 50 kw mission (1), 
for a small submersible with a search mission. The power 
profile for this mission is shown in Figure III-l. This mission 
represents a possible application for the proposed propulsion 
plant and has been used for detailed analysis of the cycle. A 
propulsion power curve versus vehicle speed is shown in Figure 
III-2 for a typical deep submersible for which the 50 kw full 
power level is applicable (2). A speed of about 7.5 knots is 
about the maximum attainable for a 70% transmission efficiency 
and a cluttered hull for the 50 kw useful power output at 
rated conditions. The hull in this case would have projecting 
search lights, manipulator arms, etc. In the analysis, the 
1000 kwh mission and 50 kw power level have been assumed to be 
Over and above any parasitic power loads peculiar to the prop- 
ulsion plant. 

Whiie the propulsion plant discussed in this thesis may 
be useful at other levels of power and endurance, no definite 
mission eeu be determined. In the interest of determining 
the behavior of the cycle at higher power levels, a 500 kw 
plant was also studied. For this power level, the major 
components are sized and component efficiencies are evaluated. 
The weight of reactants for such a plant is, however, unspecified. 
Hence detailed and overall plant weights have not been 


developed for the 500 kw plant. 


24 






J lia0dd dae NOMS! isan ols 


| SJONSYSS53Y 


ADYSANS NOISSIN IWLOL % 
OOl O08 O9 OV Od 


JSINYOD 


AEE OSE 


Ov 


= 
> 
TE, 
m 
= 
< 
m 
po) 


O08 





SHU 2° ID SAY IV LOL 


SAY OO*? INSOSV 
S4YQO°S = YSANSNVW 
SAU OO'GS .  ASINYD 
sau 2°] LN30S30 
NOISSIW YMY OOO! YOS NOILIGNOS LV SWIL 


ol. 


HAMOd GsSlVY % 






Q33dS JTDIHSA SA FAUYND YAMOd NOIS INdOUd ¢-I AYNSIA 





SLONH Ni Gs3d$ 
Ok 6 8 os oS vp 6S 





Au a!) 


@ fa[0) | GtalBIE|E =| 


TINH G3YSLLATIONA 


TINH GS¥3LLAND 


Oz 
Ob 
“L 
09 p> 
4 
08 AS 
O 
A) 
ool & 
5 
oa) we 
7 
A) 





Lo 
CN 








With the mission power profile specified and the design 
depth selected, the analysis of the propulsion plant was con- 
ducted. The first step in the selection of an optimum plant 
was to evaluate the thermodynamic characteristics of the cycle. 
From the thermodynamic assessment, one could then determine 
the most critical components, i.e., those which would be 
governing with respect to plant volume and plant efficiency. 
These components were then analyzed in detail, and the effect 
of varying the thermodynamic parameters on overall plant weight 
could be seen. For instance, it was found that an increase of 
condenser volume of about 1 cubic foot would result from a 
reduction in condenser pressure from 2.8 to 1.0 Ibf£/in* “for 
the 50 kw plant. This correspondingly would save about 6 
cubic feet of cryogenic hydrogen in a 1000 kwh mission. Com- 
parisons such as these led to the conclusion that the lowest 
attainable condenser pressure would yield the optimum plant, 
recognizing that condenser pressures less than 1.0 ibf/in- 
-were probably not attainable. 

Because of ne Lecaee number of variables and the difficulty 
of specifying many of them analytically, no overall optimiza- 
tion scheme was devised. An example of such difficulty arises 
in turbine design where the number of possible stage pressure 
ratio and stage combinations becomes very large. If the tur- 
bine could be completely analyzed by computer optimization 
techniaques this would not be so difficult, but such a scheme 
is difficult to devise. The approach used in the turbine 


analysis was to determine optimum stage characteristics by 


==) 






computer optimization and then to determine stage pressure 
ratios and rotative speed by hand calculatiois. While this 
approach of sub-optimization of system components in order to 
achieve an optimimum overall system.is not rigorous, it 
probably results in a plant which is not far from the actual 
optimum. The selection of the optimum thermcdynamic conditions 
assists in this process. Assessment of optimum thermodynamic 


conditions is subject to computer analysis. 


oo 






REFERENCES 
1. Thermal Chemical Dynamic Power Supply fox Deep Submersrble 
Vehicles, NAVSHIPS 0907-000-6010, August 7, 1970. 
2. Deep Ocean Technology Testbed Submersible Design and 


Engineering Data, Westinghouse Ocean Research and 
Engineering Center, January 6, 1968. 


29 





CHABTER Ly 


Thermodynamic Assessment 


The first step in analyzing the semi-closed Rankine cycle 
propulsion plant was to investigate the effect of various 
plant parameters of overall thermal efficiency. A schematic 
diagram of the cycle with regeneration is shown in Figure IV-l. 

An enthalpy-entropy diagram of the semi-closed Rankine 
cycle is shown in Figure IV-2. The state points in this dia- 
gram correspond to the points in Figure IV-1l and are represen- 
tative of the recirculated fluid. 

From these diagrams one may write the first law of thermo- 
dynamics for bulk flow across the control volume of the 


combustion/mixing chamber: 


Oo co} 
GE roan a A e - @ - a ci) e e - e a 
7 0 = 1e.0 m + 13 06 Mm 149 (m tm) + ¥ We (IV-1) 


where 
E = internal energy within the control volume, Btu 
0 Stagnation enthalpy of the reactants, including 
the heat of combustion, Btu/lbm 
ho om stagnation enthalpy of the recirculated diluent 
water, Btu/lbm 
Lao stagnation enthalpy of the steam leaving the 


combustion chamber, Btu/lbm 
m_, m= mass flow rates of reactants (or combustion 


products) and recirculated water respectively, 


ibm/hr 
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Oo; Ne = rate of heat addition and shaft power out 
respectively, from the control volume, Btu/7he 
By rearrangement of equation IV-1l1, one may then determine 


the recirculation ratio, 


Tm iy 7 iy Oo 
° a - 1 
es 4,0 370 


For the case of no regeneration, one may Simply replace 
i, a by 1, ae the stagnation enthalpy of the feed pump dis- 
Ul Ud 
charge. 


The first law of thermodynamics for bulk flow may also 


beswraitten for a control volume surrounding the regenerator: 


0 0 
dE - cuits e e e — e inna e e S-5f 
fe o> ig gfmotm)tls , my +350 ™r eo Mee +f x 
(IV-3) 
where 


stagnation enthalpy at expander exit, Btu/lbm 


of 
5,0 


stagnation enthalpy at condenser inlet, Btu/lbm 


of 
6,0 


Combining equations IV-2 and IV-3, one may now write 


axe: Porta) : (iso 46,0) 
—~ Si > a ae ee as eee 
m 27O@ 52,0 57 O) 26,0 


As one might expect, and is now evident from equation IV-4, 
the recirculation ratio increases with regenerator effectiveness, 
Since it is regenerator effectiveness which determines the 


a3 





magnitude of the term (i, Cae aie This results in an increase 
f f 


in pump work as regenerator effectiveness increases, but this 
is at little expense in overall thermal efficiency. The pump 


work is usually negligible when compared with expander work 


for a Rankine cycle. 


Regenerator effectiveness, as defined by Kays and London 
(1), is as follows: 


(m_+m_) c 
onan « 


_ ,steam (To 26) 


NR a 5 (Iv-4) 
ee “po, steam Bis Say 


where the product, (m +m) Cy steam in the denominator 
represents the minimum product of mass flow rate and specific 
heat. The minimum product of Me,» for this application applies to 


the steam side for all cases of practical interest. 


Equation IV-4 reduces to the following simple relationship: 


mc 


i - == (iyv-5) 
R T,-T, 


Thus one may determine the regenerator steam outlet 
temperature, if one knows the feed pump discharge temperature 


and expander exit temperature. 


= = | a - 
ee te tiene) + WR Fo (Iv-6) 


The assumed value for the enthalpy of the combustion 


products, 1. of is devendent upon the condition of the 
tA 


reactants when they come in contact with the combustion 
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chamber. The following values were computed for As ae 
a 


14,0,75 = 6869.8 Btu/lbm - Reactants at 75°F before 
contact with the combustion chamber 
a oyeey = 6559.8 Btu/lbm - Reactants at saturated liquid 
conditions at a pressure of 1 atmosphere before 
contact with the combustion chamber 
These enthalpy values reflect the same base enthalpy as 
reference 2, zero enthalpy of saturated liquid water at 32°F. 
The value of ae Ts is associated with the higher heating 
value of hydrogen at standard conditions. 
One observation which may be made already is the desira- 
bility of heating the reactants, if they are stored ina 
cryogenic state, prior to their coming into contact with the 


combustion chamber. 


Expander total to static efficiency is defined as follows: 


t4,0 75,st 


ak -1 ; 
Ors pls, st 


Ne t-s (Tv=7) 


Miere the subscript, st, denotes static conditions. 


Overall cycle thermal efficiency is defined as follows: 


_ net work _ (mitm) (24 outs st) Me Qui ) 
py heat added ae; 
Cc C,0,/159EF 


fo 


(Iv-8) 


For purposes of the thermodynamic assessment, one may neglect 


pump work, giving 


tA) 
uy 





(m +m_) (i, .-i ) 
a eee 476. (Opot af 
Noy = See ee = (Iv-9) 
c CpO nt 
Of course, in the final analysis all parasitic loads must be 
considered when evaluating the overall thermal efficiency. 
The specific reactant consumption is a somewhat better 


means of comparison of competitive plants. Specific reactant 


consumption is defined as follows: 


reactant mass flow per hour 


a HeE WOrk Per Nour 
= 2545 Btu/hr HP + 1.341 HP/kw 
Cc e e . 5 
+ a 
(m_. Mi) (4 4 is st) 
m 3412.8 
(m+m_) a, Gu57Se 


One may now show the effect of various plant parameters 
On overall thermal efficiency, bearing in mind that the 
efficiencies reported do not reflect parasites losses. First, 
one may consider the effects of varying expander inlet temper- 
ature and pressure without regeneration for fixed condenser 
pressure of 1.0 lbf/in? and expander efficiency of 75%. This 


is shown in Figure IV-3. 
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The beneficial effect of increasing the expander inlet 
temperature and pressure is evident. A considerable increase 
in overall thermal efficiency may be attained by increasing 
expander inlet temperatures above the 900-1000°F level, normal 
for most superheated steam plants, to 2000°F. Of course, the 
peak temperature attainable will be dependent upon materials 
and the particular expander geometry and operating conditions 
selected. The expander inlet pressure should also be as high 
as possible to attain the best overall thermal efficiency. 
This, too, must be tempered by constraints on the expander. 

Figure IV-4 shows the effect of expander inlet pressure on 
expander exit temperature for constant expander inlet tempera- 
ture. The result is to increase expander exit temperature with 
decreasing expander inlet pressure, thus making the plant a 
good candidate for regenerator, particularly at poor expander 
efficiencies. | 

Figure IV-5 shows the improvement of overall thermal 
efficiency with a constant regenerator effectiveness of 0.70 
at two expander inlet temperatures. Here a constant condenser 
pressure of 1.0 lbf/in? and regenerator steam-side pressure 
drop of 1.0 lbf/in? have been assumed. Again, expander 
efficiency has been maintained constant at 75%. The improve- 
ment in overall thermal efficiency at a relatively poor 
expander efficiency is evident. At expander inlet pressures 
above 750 lbf/in* with an expander inlet temperature of 1500°F, 
regeneration is not effective because the expander exit temper- 


ature is too close to the saturation temverature. The same 
38 
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effect is noted for higher expander efficiencies. This occurs 
at about 80% expander efficiency for an expander inlet temper- 
eeure of 1750°F and pressure of 600 ibf/in-. 

The effect of increased regenerator effectiveness is shown 
in Figure IV-6 for an expander efficiency of 66%, a condenser 
pressure of 2.8 lbf/in* and a constant regenerator pressure 
@repeort 0.7 lbf/in*. In this case, expander inlet temperature 
and pressure were kept constant at 1750°F and 600 Ibf/in? 
respectively. The relatively small increase in overall thermal 
efficiency with increasing regenerator effectiveness would 
suggest that the regenerator is not a critical component in 
the semi-closed Rankine cycle propulsion plant. 

Figure IV-7 illustrates the effect on overall thermal 
efficiency of decreased condenser pressure. Here expander 
efficiency has been kept constant at 75% with no regeneration. 
This effect is significant and points to the desirability of 
achieving the lowest possible condenser pressure consistent 
with reasonable condenser size. 

The effect of increased expander efficiency is also very 
significant; this is shown in Figure IV-8. Here a regenerator 
Setectiveness of 0.70 and pressure drop of 0.14 lbf/in*, a 
condenser pressure of 1.0 lbf/in”?, and expander inlet pressure 
of 600 lbf/in* have been assumed at two expander inlet tempera- 
tures, 1750 and 1400°F. From Figure IV-8 one may compare the 
two major expander candidates, turbine and reciprocating steam 
engine. If a 75% turbine efficiency is achievable at a 17/50°F 
turbine inlet temperature, a reciprocating steam engine 
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operating at an inlet temperature of 1400°F would need to have 
an efficiency of about 82.5% to be competitive. The figure of 
1400°F for a reciprocating steam engine is probably the highest 
attainable at present, as is discussed in Chapter VII. 

Thus it would appear that the following parameters will 
yield the highest overall thermal efficiencies for the semi- 
closed Rankine Cycle propulsion plant: 

a) the highest attainable expander inlet temperature 

and pressure 

b) the lowest possible condenser pressure 

c) the lowest possible regenerator steam side pressure 

drop 

ad) the highest attainable expander efficiency 

Each of these parameters must be tempered by size and 


Material limitations of the various plant components. 
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LIST OF SYMBOLS 


Specific heat, Btu/lbm°R 

internal energy, Btu 

enthalpy, Btu/lbm 

mass flow rate, lbm/hr 

rate of heat addition - Btu/hr 
specific reactant consumption lbm/kwh 
Temperature, = °R 

rate of shaft work, Btu/hr 
-turbine total to static efficiency 
overall plant thermal efficiency 


regenerator effectiveness 


SUBSCRIPTS 


refers to combustion products (or reactants) 

refers to reactants at cryogenic temperatures (1 atmos- 
phere, saturated liquid) 

refers to total conditions 

refers to recirculation flow 


refers to standard conditions 


General 


CHAPTER V 


Combustion Chamber 


A combustion chamber burning hydrogen and oxygen in 


stoichiometric mixtures must meet the following criteria: 


a) 


b) 


Cc) 


d) 


e) 


f) 


g) 


Must be safe with respect to explosions or possible 
uncontrolled oscillations 

Must operate at acceptable heat fluxes to combustion 
chamber walls and internal components 

Must exhibit very good mixing of combustion products 
and diluent to avoid “hot gas streaking” or diluent 
droplet carryover which might prove disastrous to 
turbine blading. 

Must have a reasonably long time between overhauls 
(4000-5000 hours is the design goal) 

Must operate at reasonably high pressure to achieve 
high overall plant thermal efficiencies 

Must be reasonably compact 

Must operate at very high combustion efficiency (~99%) 
using a stoichiometric mixture of hydrogen and oxygen 


with diluent water 


Several desirable features for such a combustion chamber 


may also be considered. 


a) 


Capability to operate with varying chamber pressure 
to improve expander (turbine) performance at off-design 


conditions 
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b) Capability to operate on gaseous reactants to eliminate 
the need to heat cryogenic reactants in the combustion 
chamber, thereby reducing the available enthalpy of 
the reactants 

The pressure at which a combustion chamber must operate 

in order to achieve high overall thermal efficiency is limited 
by the type of chamber selected and the reactant feed pressures 
available. The expenditure of power necessary to pump hydrogen 
and oxygen to high feed pressures with small pumps or the use 
Of gaseous reactants,stored at high pressures, limits to some 
degree the pressure attainable in the combustion chamber. 
Gaseous storage of reactants is not possible if the combustion 
chamber selected requires liquid or supercritical reactants at 
cryogenic temperatures. 

Three combustion chambers were investigated for application 

in the semi-closed Rankine cycle propulsion plant: 

a) Liquid propellant rocket-type combustion chamber 

b) Porous plug combustion chamber 

c) Catalytic chamber 

A 600 lbf/in? chamber pressure and exit temperature of 1750°F 
were selected because of the associated, high overall plant 
thermal efficiency attainable with these values and because 


they appeared to be achieveable. 


Liquid Propellant Rocket Combustion Chamber 


The liquid propeliant rocket combustion chamber is dis- 
cussed and analyzed in detail in Anpendix A. The chamber, 
shown schematically in Figure V-l, burns nydrogen and oxygen 
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in a manner similar to that of a liquid propellant rocket. 
Supercritical hydrogen is admitted through a porous liner in 
the injector head. Liquid oxygen, pumped to supercritical 
pressures upstream of the head, is sprayed in through a number 
of swirl type injectors. The oxygen enters the chamber as a 
two phase mixture since at the design combustion chamber 
pressure selected, 600 lbf/in*, oxygen is a liquid. The swirl 
injector is characterized by its ease of fabrication for small 
oxygen flow rates and its capability of being throttled to 50% 
of its design flow rate. Injectors would have to be stepped in 
and out for throttling down to power levels less than 50% of 
the maximum chamber fiow rate. This corresponds to 62.53 of 
design power since the chamber is designed to operate continu- 
ously at 125% of the design plant power level. At the 50 kw 
design power level, laser fabrication of at least the smallest 
injector hole is necessary to permit a pilot injector to 
operate continuously as injector stepping is performed. Smooth 
meepping Of injectors without flameout or pressure or tempera— 
ture excursions must be shown to be satisfactory for this 
chamber. In all the chambers, the capability of throttling to 
10% of the design mass flow was considered a requirement, even 
though the assumed power profile does not stipulate this. 

Operation of the chamber with gaseous hydrogen and oxygen 
is possible but would reduce significantly the available 
cooling to the injector head (by the cryogenic reactants). 
Such cooling might be provided by spraying,in a limited amount, 
@= diluent steam into the combustion zone, but care must be 
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taken to avoid quenching. Operation with gaseous hydrogen 
probably would result in a shorter lifetime for the combustion 
chamber and perhaps somewhat lower combustion efficiency than 


for liquid reactants, but this would have to be shown by 


experimental work. 
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Liquid Propellant Rocket Combustion Chamber 
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The heat fluxes associated with the liquid reactant 
chamber are very high. Adequate cooling of the walls is pro- 
vided by a combination of film and regenerative cooling by the 
diluent feed water. Film cooling is provided in the high 
temperature region, where both radiative and convective heat 
flux are very high. Mixing of the evaporated film coolant and 
the bulk combustion products is enhanced by turbulence pro- 
moters at the wall, themselves protected by film cooiniven 
After the point of wall film coolant dry out, regenerative- 
type coils are provided in the combustion chamber walls to 
evaporate the remaining diluent water which is then sprayed 
into the chamber from an annular plenum. Thorough mixing is 
necessary to eliminate hot gas streaking and must be evaluated 
experimentally. 

Uncontrolled oscillations are a danger in liquid propell- 
ant rocket chambers and there is no indication that the 
possibility for such oscillations does not exist in the proposed 
application. For this reason, the types of reactant and diluent 
(feed) pumps employed in conjunction with the chamber should 
have a Bei min of fluctuation in pressure and the frequency of 
any fluctuation should not be the same or a harmonic of the 
oscillation frequency of the chamber. Similarly, the exit nozzle 
of the chamber should remain choked at all operating conditions 
to prevent oscillations in the expander (pressure fluctuations 
due to turbine blade passage) from reaching the combustion 


chamber. 
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Startup of such a chamber would require a short time of 
hydrogen rich operation until ignition is achieved and proper 
film coolant flow is established. During this time, the small 
but quite hot flow would be directed to the condenser, bypass- 
ing the turbine. A small dilution chamber might be necessary 
to spray in diluent to reduce this gas température to one 
acceptable for the condenser. Shutdown would also require 
hydrogen rich operation, cutting out oxygen first followed by 
hydrogen. The only reason for hydrogen rich operation during 
shutdown, however, is the possibility of carryover of water 
into the turbine should the diluent water not be cut out 
early enough. 

Operation of the combustion chamber at higher pressures 
than 600 lbf/in*, the selected design pressure, is possible 
but may not be practical. At the 50 kw design power level, it 
is seen in Chapter VI that turbine blade height and arc of 
admission in the first stage arequite small. At higher design 
power levels, this may be of somewhat lesser importance. The 
other consideration is pumping power expended to pump liquid 
oxygen and hydrogen to suitable injection pressures. High 
injector pressure drops (100-250 lbf/in*) require quite high 
oxygen injection pressures with even a 600 lbf/in? combustion 
chamber design pressure. Operation at variable combustion 
chamber pressure might be expected to improve overall plant 
efficiency by keeping expander (turbine) efficiency nearly 
constant over a wide range of power levels, With the exception 


Of Reynolds number effects, variable pressure operation with a 


GJ 
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constant speed turbine could produce such a constant efficiency 
Bharacteristic. If 600 lbf/in* is the design pressure, this 
would mean off-design operation at lesser pressures. Table V-l 
gives predicted overall plant thermal efficiencies under this 
mode of operation for 60% of the design power level. It appears 
that the efficiency gained with variable chamber pressure 

below 600 lbf/in*® is not worth the increased complexity of the 
system. Variable pressure operation would require diluent 
(feed) and reactant pumps with a capability of wide discharge 


pressure. 


Porous Plug Combustion Chamber 

The porous plug combustion chamber uses a principle 
developed in laboratory studies of laminar flames. Siegler 
and Moore (1) have developed this device, explained in detail 
in Appendix B, for recombining radiolytically decomposed 
hydrogen and oxygen (off-gas) from the condenser of a boiling 
water reactor power plant. This device has proved to be 
Satisfactory and reliable in a 1700 hour test under virtually 
Operational conditions at a pressure of 1 atmosphere. For 
the proposed semi-closed Rankine cycle, the burner would 
Operate at high pressures. Figure V-2 is a schematic diagram 


of a high pressure combustion chamber. 
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Schematic Diagram of a High Pressure Porous Plug Combustion 


Chamber 
Parameter Design 60% power- 60% power- 
Point Variable Constant 
Combustion Combustion 
Chamber Chamber 
Pressure Pressure 
Combustion chamber pressure, 600 400 600 
i5f£/in? 
Turbine inlet temperature, °F 750 PFS 1750 
Turbine efficiency, % T3695 Voeo 68.0 
Gendenser pressure, lbf/in? 10 On7 0.7 
Overall Plant thermal 
efficiency!, $% 36.0 B52 34.8 
‘Based on turbine power, H., and O, at 75°F, and neglecting 


bearing and parasitic proptilsion Plant losses. 


TABLE V-1 
Variable Pressure Operation of the Liquid Propellant Combustion 
Seo ee tek” Plant 






The high pressure application incurs many penalties. 
First is the necessity to dilute the stoichiometric mixture 
(pre-mixed) of hydrogen and oxygen with the only inert available, 
steam, in order to achieve acceptably low burning velocities 
and heat flux to the porous plug. This requires preheating of 
the stoichiometric mixture by addition of the diluent steam to 
temperatures in the range of 400-500°F in order to prevent con- 
densation of the steam prior to entry into the porous plug. 
Evaporation of the diluent (feed) water is performed in the 
porous plug, where the diluent water is the coolant. Super- 
heating is performedinaseparate coil located downstream of the 
plug in the hot gas steam or embedded in the chamber wall. A 
portion of the superheated diluent is then added to the hydro- 
gen and oxygen and the remainder is sprayed into the chamber 
downstream of the Rote plug and superheater. 

Detailed investigation of the premixing requirement has 
shown that the temperatures involved do not constitute a 
hazard as far as spontaneous ignition is concerned. This, of 
course, considers complete and uniform mixing of the steam, 
hydrogen and oxygen. The actual method of mixing and the 
Order of mixing of the fluids as well as the velocities at 
which they are mixed areimportant. Care must obviously be 
taken to prevent the build-up of static charges from mixing 
nozzles. 

The most stringent requirement on the system appears to 
result from the required pore size of the porous plug material 


when a chamber is operated at high pressures. Reduction of 
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pore size to prevent flashback is necessary because of reduc- 
tion of quenching distance with pressure. Quenching distance, 
one finds, is nearly inversely proportional to pressure. 
Addition of the steam diluent acts to increase quenching 
distance and hence increase the pore size, but at a pressure 
fmeee00 lbf/in* this still results in a predicted quenching 
distance of about 37y. If a safety factor of 10 is used to 
insure against flashback, this results in a pore size of about 
3-4. 

A minimum flame velocity of about 2-3 cm/sec is imposed 
by flammability limits. At 600 lbf/in* chamber pressure, it 
is shown in Appendix B that flame velocities greater than 4.0 
cm/sec result in excessive heat flux to the porous plug, the 
effect of which is to heat the downstream side of the plug to 
temperatures in excess of that which can be tolerated for a 
Sintered nickel material. 

Reduction of pore size is detrimental, however, to gas 
Side pressure drop across the plug if adequate effective thermal 
Seneauctivity of ae porous material is to be maintained to 
transfer the heat. Several correlations were used to estimate 
gas side pressure drop at the design conditions (600 lb£/in? 
chamber exit pressure, + inch thick porous plug), and resulted 
in a range of predictions from 30 to 1400 lbf/in?. Experience 
with porous filtration materials of the sintered variety would 
indicate that the lower eeorse, 30 Ibf/in- ws most likely 
accurate, based upon experimental data in the range of interest 


f)- if reduction of pressure drop is necessary, it can be 
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achieved to some degree by decreasing plug thickness within the 
limitations of coolant tube diameter and required thickness 
downstream of the coolant tubes to achieve a uniform flow 
velocity. Such a design would be desirable in any case to 
reduce the plug exit temperature to the lowest achievable value. 
This is important with respect to oxidation of the nickel which 
could close up the pores after some period of operation. It is 
doubtful that a material with an oxidation resistance superior 
to nickel could be found which would meet the other requirements 
of sufficient strength and adequate thermal conductivity. 
Another problem associated with the porous plug burner 
is the necessity to divide the burner into segments in order 
to throttle the chamber. This aca consequence of the 
previously described limitations on flame velocity of 3-4 cm/sec. 
This implies throttling of the chamber is not possible below 
753 of its design flow rate unless the burner is split into 
Sections and sections are "stepped out" as throttling continues. 
Such an arrangement would probably require some cooling steam 
to pass through the unused sections of the plug in order to 
Meeain plus exit temperature in these sections considerably 
below the sintering temperature. This imposes an added 
requirement for mixing downstream of the plug of the coolant 
steam and the combustion gases. 
In the proposed ee oeationt slightly improved overall 
thermal efficiency could be obtained in this device over that 
Of the liquid propellant rocket combustion chamber. Cryogenic 
reactants would be heated by sea water (through an intermediate 
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heat transfer loop) and would enter the chamber as gases. 


Satalytic Combustion Chamber 


A catalytic combustion chamber was investigated as a 
possible alternative to the two combustion chambers previously 
described. This type of chamber was developed for use in the 
Dyna-Soar project and employed excess hydrogen to achieve 
acceptable temperatures. This chamber is discussed by Bailey 
(3). In the method of application which is of interest with 
respect to the semi-closed Rankine cycle, hydrogen and oxygen 
were mixed as gases at 540°R upstream of the catalyst bed by 
an injector arrangement similar to that for a liquid propellant 
rocket combustion chamber. A concentric injector arrangement 
was tested, as well as one employing impinging jets, in order 
to yield adequate mixing upstream of the catalyst bed in as 
small a volume as possible. The impinging jet arrangement 
produced the most satisfactory results. Flashback occurred to 
the injector head at high mixture ratios of oxygen to hydrogen 
but was non-destructive. 

The 3" diameter prototype chamber developed at Sundstrand 
Aviation - Denver, operated at 275 lbf/in* and a discharge 
temperature of 1960°R with a mixture ratio of 0.8 pound oxygen 
per pound hydrogen. A packed catalyst bed 4 inches in length 
was used. Measured combustion efficiencies were 9911, .9730, 
-9455 and .9080. No explanation was given for the variation 


an efficiency. 
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Bailey (3) reported that the palladium coated aluminum 
oxide pellets used as the catalyst had a melt*ng point of 3289 
and 4180°R for the palladium and substrate respectively. He 
also reported that catalytic action was not impeded by the 
presence of a liquid water film on the catalyst. In another 
combustion application, water condensed on the catalyst bed 
and did not impede startup. 

The literature was searched for reports of work conducted 
on this combustion chamber subsequent to the writing of 
reference 3. The most recent report found was one on the Dyna 
Soar APU System, apparently the last report prior to cancell- 
ation of the Dyna Soar project (4). 

The catalyst described in reference 3 is a sintered 
aluminum oxide pellet, coated with palladium. The required 
life of the pellets for the Dyna Soar application was 250 hours. 
In reference 4 is described an evaluation program of several 
substrates for the palladium catalyst. Life tests which were 
performed invariably resulted in overtemperature of the 
catalyst bed, from either human error or component failure. 
Nevertheless, the results obtained were not encouraging. Only 
one catalyst met the 250 hour lifetime requirement of the Dyna 
Soar APU system, and its performance was somewhat degraded at 
the end of testing. It was, however, puptecred to several 
temperature excursions. In one test, water was inadvertently 
allowed to flow into the combustor from the exhaust system. 
Even after vacuum drying, ignition was not possible with this 
Catalyst bed with a stoichiometric mixture of hydrogen and oxygen. 
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While the catalytic combustion chamber may prove to be 
successful for a short lifetime (250 hours) with hydrogen rich 
operation, it is very doubtful that it could be utilized in 
the semi-closed Rankine cycle of this thesis. It is questionable 
if sufficiently complete combustion could be achieved in a 
catalyst bed of reasonable size with the amount of diluent 
steam necessary in this application to reduce the temperature 
of the catalyst bed to acceptable levels. The loss of catalytic 
action when the catalyst bed was exposed to water in significant 
amounts is also not encouraging for this application. Further- 
more, because of the relatively short lifetime of the tested 
Betalyst configurations, it is doubtful that a catalyst could 
be developed having a lifetime sufficiently long to permit an 


extended time between overhauls for the propulsion plant. 


Combustion Chamber Selection 

mie porous plug burner has the greatest promise LOG 
extended times between overhauls and reliable operation of the 
three Penbustion chambers. However, the combination of possible 
Oxidation problems in the porous plug, questionable porous plug 
pressure drop, and low effective thermal conductivity of the 
best available porous material would require testing to deter- 
Mine if these problems can be overcome. Estimated quenching 
distance and laminar flame velocities must also be verified 


Since they are critical to the successful design of the burner. 





The liguid rocket combustion chamber, though probably less 
efficient a combustor than the porous plug burner, appears to 
be the best short term solution to the combustion chamber 
problem. It is expected, however, that a design lifetime of 
4000-5000 hours probably cannot be met with this chamber and 
that a lifetime of about 1000 hours is more realistic. This 
would require opening of an access to the propulsion plant 
pressure vessel to replace this component, which is undesirable, 
but possible. 

A listing of weights and volumes for a 50 kw and 500 kw 
design power level liquid propellant rocket combustion chamber 


are given in Table V-2. 


Parameter Design Power Level 

50 kw 500 kw 
Inside Length, in 1520 30.0 
Inside diameter, in 4.0 8.0 


Outside diameter, (including 


Masulation) in ORIG ar 0 
Overall length Tee0 36.0 
Wolume, ft? 2-62 3.2 
Weight, lbm 50.0 300 

TABLE V-2 


Liquid Rocket Combustion Chamber Dimensions, 


Weights and Volumes 





REFERENCES 


Siegler, M. and G.E. Moore, "Flame Recombination of Oxygen 
and Hydrogen", Chemical Engineering Progress Symposcium 
Series, Nuclear Engineering - Part XXI, Vol. 66, 1970, 
No. 104. 


Seevens, D., Pall Trinity Micro Corporation, Cortland, 
New York, personal communication, 11 May 1972. 


Bailey, R.N., "Development of Catalytic Hydrogen-Oxygen 
Reaction Chambers for Space Power Systems", Space Power 
Systems, Vol. 4 of Progress tn Adtnonautics and Rocketry, 
N.W. Snyder, ed. New York, Academic Press, pp 183-210, 
1916 2 « 


Design Analysis Report, Dyna-Soanr APU Thermodynamic System, 
Sunstrand Aviation, Denver, Report No. 24-DER-63, 29 
April 1963 (AD435892). 





CHAPTER VI 


Condenser and Related Equipment 


Because of its anticipated large size in relation to other 


propulsion plant components, considerable effort was expended 


in the development of an acceptable condenser and condensing 


system. 
warrant 


1) 


2) 


3) 


4) 


5) 


There are several vital aspects of the problem which 
consideration: 
The size and weight of the condenser based upon an 
estimated amount of non-condensables 
Selection of the best method of heat rejection from 
the condenser 
Elimination of non-condensables 
Sensing of the amount and type of non-condensables 


Storage of product water 


Condenser Type 


The selection of condenser type was predicated upon the 


requirement that the device be as compact as possible. Several 


types of condenser arrangements were considered. 


1) A standard horizontal tube and shell condenser 


2) 
3) 


4) 


5) 


A vertical tube and sheli condenser 

Condensation on the inside of vertical tubes, located 
outside the propulsion plant pressure vessel wall 
Condensation on the inside wall of the propulsion plant 
pressure vessel 


Spray type condensing 


Condensation directly on the inner surface of the propulsion 
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plant pressure wall seems attractive at first glance, particul- 
arly when it is compared with heat transfer through the same 
wall with a standard tube and shell condenser in conjunction 
with an intermediate heat transfer loop. Non-vertical walls 
may promote improved condensation on the upper portion of the 
inner surface of the sphere to a certain extent. However, if 
the steam flow is vertically downward in the lower portion of 
the sphere, one may anticipate degraded performance due to 
condensate film thickening. Some novel ways of removing the 
condensate might be devised for the lower portion of the 

Sphere or some novel arrangement of steam flow might be used 
over the entire surface. A set of parallel spiral passages 
might be used, for example, to remove the condensate and at 

the same time permit utilization of the major portion of the 
internal surface area of the sphere. Figure VI-1 shows a 
possible arrangement of such a pressure vessel wall condenser. 
Corrosion on the inner surface of the sphere, however, could 
not be adequately controlled in the presence of pure steam. 
Corrosion inhibitors could not be used here but they could be 
used in an intermediate heat transfer loop. Furthermore, direct 
contact of superheated steam with the pressure vessel wall could 
produce excessive thermal stresses in the vessel wall. Such 
Stresses would have to be carefully analyzed to prove the 


feasibility of such an arrangement. 
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FIGURE VI-1 
Possible Arrangement for Condensation Directly on the 


Inner Surface of the Pressure Vessel Wall 


Condensation on the inside of many vertically oriented 
tubes located outside the pressure vessel aiso has certain 
advantages. The primary advantage of this arrangement, shown 
in Figure VI-2, is the reduction of required propulsion plant 
pressure vessel volume by that amount necessary for an inter- 
nally mounted condenser. A disadvantage of this system is the 
Size (or number) of penetrations to the pressure vessel which 
become necessary to duct the exhaust steam from the turbine or 
regenerator to the external condenser. The effect of these 
penetrations, of course, is to increase pressure vessel weight, 
probably offsetting, to some degree, the Significant advantage 


Gained by elimination of the space requirement of an internally 


mounted condenser. 
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External Condenser Arrangement 


An additional disadvantage of the external condenser is 
the generally larger tubes or larger number of tubes which must 
be employed (to accept the high volumetric flow rate of the 
steam inside the tubes) than for a comparable internally 
mounted horizontal tube and shell condenser. Unless even 
larger tube area is used, forced convection of seawater must 
mememoloyed, requiring some sort of external shell to direct 
forced circulation across the vertical tube bank. Because of 
the compressive load to which the condenser tubes would be 
@eegected, their thickness would be greater than for tubes of 
equal internal diameter subject to internal sea water pressure. 
Each of these items adds weight to the system. The greatest 
disadvantage to this system, however, is the proximity of 
Superheated steam with the vressure vessel wall at the steam 
Penetration. While insulation could be provided at this point, 


67 





it would have to be very carefully engineered to prevent 
excessive thermal stresses in the pressure vessel. 

A spray-type condenser was considered primarily for de- 
superheating. No standard methods for analyzing this type of 
condenser were found. For the desuperheating process, however, 
it was anticipated that a droplet vaporization model similar to 
those described in Appendix A (for droplet vaporization in 
liquid propellant rockets) might be used. Subsequent analysis 
of the tube and shell condenser (see Appendix D) showed that 
desuperheating was not a significant problem. Furthermore, 
analysis of droplet condensation in the presence of a signifi- 
cant amount of non-condensable gas appeared to be particularly 
difficult to analyze. The spray type condenser would use a 
Spray of subcooled condensate to produce condensation, thus 
adding another piece of heat transfer equipment, a condensate 
cooler. 

The standard horizontal or vertical tube and shell 
condenser, mounted inside the propulsion plant pressure vessel 
is probably the easiest condenser to analyze. The standard 
method of analysis for this senaeh couanes proved to be fairly 
accurate in predicting condenser performance. The advantages 
of the horizontal tube and shell condenser over the vertical 
tube and shell condenser are probably not too significant. 
Powell (1) has shown that the size of the two types of con- 
densers for identical conditions of the entering steam is 
Comparable. Kern (2) indicates that the vertical tube and 
Shell condenser is advantageous when condensate subcooling is 
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desired but that the horizontal tube and shell condenser is 
probably somewhat smaller for equivalent entering steam 
conditions. 

The choice of condenser type, then, is not clear cut. 
Several of the alternatives require detailed analysis of the 
effect of the condenser type on other components. For instance, 
while the external vertical tube condenser appears promising, 
the problems of the pressure vessel wall penetration for steam 
would require detailed stress analysis. Furthermore, the 
reduction of pressure vessel weight due to elimination of 
required internal volume for a condenser is offset by gains in 
tube weight and pressure vessel weight at the penetrations. 
Condensation on the internal wall of the pressure vessel appears 
@engerous in that it might result in catastrophic thermal 
stresses in the wall. 

- For the reasons outlined above, the standard tube and 
shell condenser was selected for the purposes of this thesis. 
While some other choices appeared possible and perhaps even 
Slightly better choices as far as overall propulsion plant and 
encapsulation weight is concerned, this could only be proved 


after much detailed analysis. 


meat Rejection 


The selection of a condenser located within the pressure 
vessel hull imposes a requirement for a means of rejection of 
heat to the sea water. Two alternatives were considered: 


1) a skin cooler (heat rejection through the hull of the 
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pressure vessel) similar to that used on the U.S. 
Navy's deep diving submarine Dolphin 
2) a conventional, but heavy, standard forced sea water 
arrangement (requiring penetration of the pressure 
vessel hull of sea water piping) 
The first alternative, that of the skin cooler, appeared 
to be most attractive since it required no pressure vessel 
penetrations. Since the propulsion plant pressure vessel would 
be mounted inside the envelope of the streamlined body of the 
submersible, natural convection heat transfer for the outer 
surface could not be employed as in Dolphin unless unacceptably ° 
high temperatures of the intermediate neat transfer loop were 
permitted. Hence, forced circulation on both inside and outside 
of the pressure vessel wall was considered. Figure VI-3 is a 
schematic diagram of such an arrangement. The intermediate 
heat: transfer loop circulates fresh water, containing a 
corrosion inhibitor (and perhaps ethylene glycol anti-freeze 
if desired), through a standard tube and shell condenser and 
then rejects heat by contact with the cool pressure vessel 
walls. Heat is then removed by the external forced circulation 
loop employing sea water as a fluid. A submerged sea water 
Circulating pump, located outside the propulsion plant pressure 
vessel, provides the pumping power for the forced circulation. 
It was estimated, for purposes of a design study, that as 
much as 80% of the pressure vessel wall area could be used for 
the skin cooler. The remaining 26% of the area would be 
required for access, penetrations, etc. To determine feasibility 
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of such an arrangement several assumptions were made, some of 


them slightly optimistic: 


a) Outside sea water temperature was assumed constant 


b) 


Cc) 


d) 


through the external loop at 40°F, the sea water 
temperature at extreme depths. This assumption would 
result in slightly degraded plant efficiency at the 
surface. 

The intermediate loop temperature was assumed to be 
Sonstant throughout the @loop at 75°F. his as a 
reasonable assumption for condenser size as will be 
shown later. 

A coolant velocity of 7 ft/sec was assumed to produce 
an acceptable film heat transfer coefficient on both 
the inside and outside of the pressure vessel walls. 
Higher velocities, while possible, would result in 
Significant increases in pumping power required. 

Three possible cooling passage arrangements (A, B, and 
C-in Figure VI-3 were considered.) Arrangement A is a 
1" x 4" passage, long dimension parallel to the 
pressure vessel wall. Arrangement B is a modification 
of the geometry of arrangement A and includes 4" fins 
on both sides. 90-10 Copper Nickel was selected for 
the salt water side and pure copper for the fresh water 
side. The employment of cathodic protection could solve 
the corrosion problem which would exist in this 
arrangement. Arrangement C utilizes passages of the 


same size as Arrangement A, but with the long axis 
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e) 


perpendicular to the pressure vessel wall, and has 

fins of the same material as the pressure vessel wall. 
For the finned arrangements, it was assumed that the 
method of attachment of the fins would not affect the 
fin efficiency. This, of course, is untrue. Welding 
of ferrous to non-ferrous materials (copper base 
Materials to steel, for instance) is a possible method 
of attachment. Such a bond as this might not affect 

the fin efficiency. Certainly welding of fins of the 
Same material as the pressure hull would be most 
reliable in not degrading fin efficiency. Welding to 
the pressure vessel walls, however, is not desirable, 

if it can be avoided. Certain thermally conductive 
epoxy materials have been successfully used for 
attachment of heat rejection devices to deep submergence 
pressure walls (3). Two such materials are STYCAST 2850 
and DEVCON. The thermal conductivities of these materials 
mee not reported by the Manufacturer aaa must be deter— 
Mined. Furthermore, the thermal contact resistances 
associated with the application of fin attachment, are 
unknown. An additional requirement for this application 
is that the adhesive be capable of withstanding the 
cyclic flexing of the pressure vessel hull in compression 
without failure of the bond. 

A heat load of 3 x 10° Btu/hr was selected for the 
feasibility analysis. This is representative of a 


40 kw propulsion plant operating at 31.5% overali 
eermlalwerriciency. 
13 





=) 


g) 


h) 


i) 


3) 


A 3-4 foot I.D. propulsion plant pressure vessel is 
capable of housing a propulsion plant corresponding to 
the power plant ae above. 

For ail hull materials, no allowance was made for 
protective coatings. Aluminum, known to be subject to 
stress corrosion, normally employs a protective coating 
of several layers of plastic, and steel hulls would be 
painted, at least on the outside. 

Assumed scale coefficients: 

= SOOM Guy Wie hyn 


h 
sc,Sw 
h = 2000 Btu/hr ft? °R 
sc,fw 
Corrosion inhibitors can be employed on the fresh water 
side, which will permit the maintenance of a high scale 


eoerficient on the inside surface. 


Pressure vessel hull thermal conductivities: 


Steel (HY-130 and HY-180) - k 21 Btu/hr £e°R (AY =s0ndaca) 


me HY 130 = k 4.2 Btu/hr ft°R (based 


upon data in reference 4 for similar alloys) 

eo «6d 7079 - T6 k = 70.2 Btu/hr £t°R (based 
on data in reference 4 for similar alloys) 

Actually, one would anticipate a slightly lower thermal 
conductivity for HY-130 and HY-180 steels than that of 
fe-oO Steel. 

The total and usable pressure vessel areas are given 


below, based upon internal sphere diameter: 
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Sphere I.D. Total Surface Area Usable Surface 


(£t) (fee) Area 

Gite) 
aye 8) 2622 22.0 
4.0 50.4 40.3 


With the above assumptions, calculations of required 
surface area were performed for cooling passage configurations 
A, B, and C of Figure VI-3. These results are given in Tables 
feet through VI-3. 

From Table VI-1, the case without fins, it is evident that 
none of the pressure hull materials would permit transfer of 
the assumed heat load with a 35°R terminal temperature 
difference. Aluminum exhibited the best performance, but this 
did not include any allowance for a protective coating. Craven 
(5) discusses the problem of stress corrosion cracking in alum- 
inum pressure vessels in an underseas environment, which has 
precluded welding of heavy pressure vessel castings. Aluminaut, 
which utilized a bolted 7079-T6 aluminum construction, also 
used several layers of a plastic coating on the exterior of the 
hull to guard against stress corrosion cracking. 

In subsequent calculations only steel was considered for 
the pressure hull material. The poor thermal conductivity of 
titanium eliminates its use in the application of a skin cooler. 
The poor thermal conductivities of plastic coatings would 
eliminate aluminum. 

Table VI-2 is an analysis of the second configuration, 


arrangement B of Figure VI-3. Improved verformance is noted, 
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Pestgn e&ressure Pressure Pcs ote 19 yp re Uveralrl Reece 
Deoth Vessel Vessel Vessel saad tepcb luda! Heat Trans- Surface 
Ste Internal Material Thickness Btu/hr ft7°R Btu/hr £t*°R fer Coeffi- Area 
Diameter oor event fee 
ft Btu/hr £t7°R 
8090 320 Steel ums 1260 is 30 169 50.7 
FY = 10S Ore te 
HY = E30. 
T2 HY=120 Oi, S622 152-5 
A1l-7079-T6 0.9 216 Bo 7 
4.0 Steel Oud 149 57S5 
HY-130 & 
ny poe 
Ti HY-120 0.9 45.7 188.0 
A1l-7079-T6 ee 7, ¥ 186 46.0 
20,0100 BO Steel Bae Ie 1260 Tiss) desz 64.9 
HY-130 l 
Steel 8 an 60.7 
HY-180 
T1-HY-120 ilies 420.5 2A? 
. Al-7079-T6 20 ESs Bid alt 
ao (@ Steel inee3 fs ee) i ogres: 
HY-130 
Steel er ALAA, y glee 
HY-180 : 
Ti~-HY-120 eS | 29753 293 
Al-7079-T6 Pla eet | 142 G0) 42 
Heat Load-3 x 10° Btu/hr Passage Size-1" x 4" 
RECS Watkeierd emperatire—7 5° FR no fins 
Salt Water Temperature - 40°F Pressure vessel» thickness (safety factor 


of 1.5) from reference 6 


TABLE VI-1 


Skin Cooler Analysis for Coolant Passage Configuration A of Figure VI-3 





z+ 
eoly 
asodejyans 
perzatnbey 


€-IA eanbty 


JO g UOoTReANHTzZUOD sbhesseg juejlooD zrA0F stsATeuy AeTOOD utyxs 


C-IA ATE 


u¥/T X wIT/E SZTS ebesseg 
»9T/T ssouyoTtyL uta 


epTsS 797eM YSerzd -ND ,F 
eptS 797FeM ATES TN-ND ,F SuTd 


8ST aE 
97T Ccae 
G6T 08° 
yere LT6° 
OTZ Le 
7S2 OSST O8ET S° 
YoztF AU/nIig 
RZUuSeTS “th 
-TFJFSOD AOF YozAF AU/NAG Yoz3F AU/NIq ssouyoTUL 
-SUBPA]L, 2eOH Tesson 
TTP 2A9A0 Ag UT Ty AS SULT eAansselidg 


d1,o07 eaAnjeroduey Ae}eEM 3TeS 


AoGL sAanjzeradusay zsjeM Ysery 


zu/njg ,OT X € = peoTy }eeH 


O8T-AH 
Tee4sS - 
O€T-AH oy 
Tee4sS 
O8T-AH 000° 0¢ 
Tee4sS : 
O€T-AH De 
Tee 4s 
O8T-AH 
20€T-xulL 
Tee3s¢r O°? 
eee 0008 
80€1-AH 
Feance Oe 
ala 
Asajowe tq 
TetzejeW TeureAuL 27 
Tessoan Tessa yjdaq 


eansserig oanssexzg uUubtsedqd 


a7 





¢-IA eanbTy FO D uoTZeANbtjuoy sbhesseg JueTOOD AO;J stsATeuy AeTOOD uTxXS 


zt 
Poly 
aoejans 
peaztnbay 


nl * av/l 22tS ebessed 


u8/T SseuyoTtyy uta 


€-IA ATaVh, 


1,07 etnqZerZSeduay zr904eM ATeS 


JIoGL eaznjerodusy, reAZeM YserzT 


Teeys ,T = sutd ay/nig >OT X € PROT 3eAH 

O8T-AH 

O° TVT Olt Te94S 

: Ve l-An 

O° OCT Sy Aaa Te9e3s 0°? 

O08 1-20 

O° real 08°0 T8e4S 

“ OC LT—-AH 
O° 2ST Evo 0 T8935) Oe 000‘02 

OS T—-An 

; 7 Ole Le 

Geert OL C T?94S OT 
O23 TA 

S Ce l-—An 

O ele O€ST 09ZT 0S°0 T2938 O£ 0008 
Yoz4z AU/NIE ait 
ZUSTIOTFFSOO no i § ZoRzoueTtTG 
ASFSueAL Yoz35 ay/nig Words azyu/njag sseuyotyuL Tetz9zeW TeurS RUT alas 

Zea i F Tesssa Tesssa Tesson yjzded 
TT e429A0 meer UL TS y ee Ut TF yy oansSeag eansserzg eAnssexrq uUbtsed 


78 





Pat it is still not sufficiently ancreased to permit transfer 
of the assumed heat load at the assumed tempevature difference. 

Table VI-3 shows the results of using 1 inch steel fins. 
Such fins could be welded to the pressure vessel or might be 
machined from the casting. Machining would prove to be a 
difficult task, however, when one considers the effort expended 
to achieve near perfect sphericity in such a pressure vessel. 
Performance in this configuration is better than for arrangement 
A but inferior to arrangement B. 

Forced circulation of the intermediate coolant (fresh 
water) and sea water requires the expenditure of considerable 
electrical power for pumping. If a 70% pump efficiency and 
60% motor efficiency are assumed for both the salt water and 
fresh water pumps, the estimated electrical power to pump 
coolant through the cooling passages and associated piping for 
eet £.D. sphere is aS listed in Table VI-4. Power required to 
pump coolant through the condenser is not included here, since 
it is common to the alternate arrangement, forced circulation 
of sea water through the condenser tubes. Condenser pressure 
drop is small, however, in comparison to the pressure drop in 
the internal panels. 

For any detailed analysis of the skin cooler, of course, 
pumping power required and its effect cn plant thermal efficiency 
must be evaluated. A 1.3 kw drain on power produced would not 


represent a serious increase in reactant load for a 1000 kwh 


mission. 





Arrangement Electrical Power Required 
for Pumping (kw) 


A UES: 

B 2.5 

C 335) 
TABLE VI-4 


Electrical Power Required for Skin Cooler Pumps 


(neglecting condenser pressure drop) - 3 ft I.D. sphere 


The skin cooler could, of course, be employed easily ie G| 
sufficiently high terminal temperature difference is provided. 
For a condenser designed for a 75°F coolant temperature, this 
means a higher condenser pressure and reduced overall thermal 
efficiency. Hever ees it can be seen that at the low condenser 
pressures shown in Chapter IV to give high overall thermal 
efficiencies, large increases in condenser volume will be 
necessary to provide a sufficiently high terminal temperature 
difference to reject the assumed heat load. 

An estimate of the additional weight incurred by the 
inclusion of a skin cooler without fins (Arrangement A) is 
given in Table VI-5. Without a knowledge of the condenser 
geometry and required flow rates of intermediate coolant loop 
fresh water, one cannot yet estimate the weight and total 
power requirements of a fresh water circulating pump, or the 


weight of the intermediate coolant. 
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Sphere Internal Diameter Weight 


(£2) (1lbm) 
oO 0 Panels @ ribs” 290 
Piping 40 
External SW 24 


pump & motor?” (308GPM) 


Total 354 

4.0 Panels & ribs 515 
Piping ae: 

External SW 43 


pump & motor (410 GPM) 
Botal Gia 
Mee” steel panels & ribs 


* 400 cycle geared, sea water compensated (oil filled) 
Pnduction motor 


TABLES VY LS 


Additional Weight Incurred for a Skin Cooler Without Fins 


The panels and ribs which form the cooling passages should 
be removable for cleaning and preservation of the pressure 
vessel surfaces. While some other lighter material than steel 
Might be found for this application, steel was selected for a 
first approximation of cooler weight. The large volumetric 
flow rate of the sea water cooling pump results in a minimal 
rise in temperature of the sea water. Thus the assumption of 
a constant temperature of 40°F for heat rejection is reasonable 


for extreme devths. 
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Non-Condensable Gases 

The amount of non-condensable gases which the condenser 
must accommodate is dependent upon combustion efficiency, the 
achievable accuracy of reactant proportioning, leakage into the 
condenser by the normal routes of turbine glands and condenser 
gaskets, and the purity of the reactants. Information obtained 
from oxygen suppliers indicates that the commercial grade of 
oxygen has a Significant amount of contaminants. Table VI-6 


gives a summary of commercial grade oxygen contaminants (7). 


Contaminant PPM 
Ny 150 
Ar 3500 
Hydrocarbons <20 
co. 5 
Moisture 6 
Total Non-Condensables 8675 
TABLE VI-6 


Commercial Oxygen Contaminants-Typical Analysis 


By far, argon is the most serious oxygen contaminant to be 
considered. Moisture, of course, does not present a problem in 
the condenser. Since argon is not condensable at normal temp- 
eratures and pressures and is non-reacting, it will represent 
a Significant disposal problem. Only research grade oxygen 


Gives any significant improvement in the amount of non-conden- 


82 





sables and this is at an unacceptable increase in cost. 


Hydrogen contaminants are given in Table VI-7 


Contaminant PPM 
co, ik 

Hydrocarbons 3 
N, 2) 
0, zZ 
CO al 

Moisture 1 

Total non-condensables 10 

TABLE VI-7 


Commercial Hydrogen Contaminants-Typical Analysis 


Combustion efficiency is difficult to assess for a 
combustor which has never been built. However, the porous plug 
combustion chamber,described in Appendix B, is expected to have 
a very Bich combustion efficiency, on the order of 99-100% for 
Steady operation. For this reason a 99% combustion efficiency 
has been selected for condenser design purposes. The accuracy 
of proportioning equipment at the small reactant flow rates for 
the 50 kw propulsion plant is unknown, but it is estimated that 
control of the hydrogen and oxygen flow rates without feedback 
could be as accurate as 1% (8). Leakage of non-condensables 
into the condenser through normal means is estimated to be 


meet .05%3 of the steam mass flow rate (9), (10). 
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The selection of a design figure which takes into account 
all of the possible sources of contaminants is somewhat arbitrary. 
Clearly the effects of combustion inefficiency and non-stoichio- 
metric proportioning are interrelated. The total effect is not 
simply additive. It was assumed, therefore, that unreacted 
hydrogen and oxygen flow rates of 2% of the design flow rate of 
each reactant would reach the condenser, plus an additional 
oxygen flow rate of .05% of the steam flow rate. This was 
considered to be a very conservative estimate of the non- 
condensables from all sources. 

Comparison of the effect of hydrogen and an equal mole 
fraction of oxygen on condensation rates showed that condensa- 
tion rates were somewhat greater in the presence of hydrogen. 
This might be expected from the higher diffusion coefficient 
for hydrogen-steam mixtures than for comparable oxygen-steam 
mixtures. 

Computer costs for the analysis of condensation in the 
presence of more than one non-condensable gas (see Appendix D) 
were significantly higher than for just one non-condensable 
gas. For this reason computer analysis of condenser size and 

weight was conducted with oxygen Only as Ehe non-condensable. 
To simulate the presence of hydrogen, an equal number of moles 
of oxygen were substituted. An allowance was also made for an 
additional amount of oxygen, .05% of the steam mass flow rate 
by weight. Substitution of an equal molal flow rate of oxygen 
for hydrogen resulted in an approximate 8% increase in condenser 


murrace area for the steam conditions tested. Thus this 
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assumption is conservative. The ratio of reactant to steam 
flow rates were based upon a representative design at 31.5% 
overall thermal efficiency (11). The resultant design figure 
for non-condensables was an oxygen flow rate of 1.39% of the 


steam mass flow rate. 


Analysis Of the Horizontal Tube and Shell Condenser 


The horizontal tube and shell condenser was analyzed for 
the assumed non-condensable flow rate at various inlet pressures 
and temperatures of interest and for a wide range of mass flow 
rates of entering steam. Condenser weights were determined for 
the following cases: 

1) a standard condenser to be used in conjunction with a 

skin cooler 

2) an 8000 foot depth "deep-sea" condenser circulating 

sea water 

3) a 20,000 foot depth "deep-sea" condenser circulating 

sea water 

The computer program developed and described in detail in 
Appendix D was used to determine condenser surface area and 
dimensions. 

Table VI-8 gives a summary of the condenser geometry 


assumed for the family of condensers studied in this thesis. 
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Tube Material (deep sea condensers only) Ineo Alloy CA7 i? 


muibe O.D.7- 11. 2250 

mube [.D.. in. (20,000 £t depth) ESO 

Tube Spacing, in. (see Appendix D) ~6605 

Tube Arrangement Staggered Rows, 
Square Pitch 
Rotated 

Coolant (sea water) velocity, ft/sec ea) 

Coolant (sea water) inlet temperature, °F 75.0 

TAB |Vl—6 


Assumed Data for Condenser Studies 


Tube wall thickness for the 8000 ft condensers and the 
Standard condenser were proportionately less than the .030 inch 
figure used for the 20,000 ft depth condenser. 

Figures VI-4 and VI-5 give condenser surface areas required 
for a range of condenser inlet BecearneS Gon 1.0) to 1020 beni. 
Two curves are presented for each pressure, one for saturated 
steam and a second for superheated steam at 300° superheat. 
Studies of the effect of superheat indicated that the additional 
required surface area for superheated steam was relatively 
small. Calculations at 100, 200, and 300°F superheat showed no 
more than a 5-10% increase in required surface area over that 
for saturated steam at the same pressure and mass flow rate. 

The standard condenser, for use in conjunction with the 
Skin cooler was designed using the specifications of references 


12 and 13 as a guide in determining shell and header weights, 
86 
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hot well volume, etc. The weights presented in Figure VI-6 
represent those of a single pass condenser and include a margin 
of 25%. A double pass condenser might prove to be less heavy 
in the tradeoff of condenser weight and volume with intermediate 
coolant weight, but this was not investigated. Figure VI-7 is 
a diagram of the standard condenser. 

The deep sea condenser uses a header (waterbox) design 
Similar to that of header or steam drum design for a boiler. 
The cylindrical headers are kept small, usually 6" or less in 
internal diameter to minimize wall thickness and the weight. 

No detailed stress analysis was performed on the headers but a 
safety factor of 2 based on yield stress was applied because of 
the large number of penetrations. The tubes and headers are of 
the same material and can be joined by MIG welding techniques 
(14). The use of non-ferrous materials for the condenser tubes 
is desirable from the aspect of scale buildup and fouling. If 
the tube and headers materials are identical, this would then 
entail joining ferrous to non-ferrous materials in order to 
connect the piping to the pressure vessel wall. Such a joining 
process could be accomolished by MIG welding techniques and 
performed under laboratory conditions if necessary (14). If 
Such joining proves to be unsuccessful, ferrous materials 
could be used throughout the condenser with a consequent in- 
Crease in condenser size due to steam side and waterside 
corrosion. It would be wise to minimize the size of the hull 
Piping penetration because of the effect of the size of the 


penetration on pressure vessel weight. Too small a penetration, 
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however, would result in the excessive expenditure of pumping 
power. Figure VI-8 is a schematic diagram of the deep-sea 
condenser. 

A detailed stress analysis must be performed on this 
design to determine the combined effect on piping penetrations 
and header joints by thermal stresses resulting from the slight 
heating of the condenser tubes and by stresses due to the com- 
pression of the pressure vessel. An alternative design, and 
probably a much simpler design to fabricate, is shown in 
Figure VI-8. This design is basically a two pass condenser 
which employs U-bend shaped tubes to provide the return pass. 
It is believed that this design might be somewhat better inso- 
far as stresses are concerned than the single pass condenser. 
Impingement of superheated steam directly on the tubes may 
require detailed stress analysis for startup conditions, and 
insulation of the headers from direct contact with the steam 
may be advisable. Under normal operating conditions the tubes 
would be covered by a film of condensate, and the outside wall 
would reach a temperature no greater than the saturation 
temperature of the entering steam. | 

In the computation of overall condenser weights for the 
deep-sea application, an additional margin of 50% has been 
applied to the overall condenser weight for piping and founda- 
tions and additional. pressure hull vessel weight as a result of 

‘ ra) 
the piping penetrations. Thus the condenser weights presented 
On Figure VI-9 and VI-10 represent what are believed to be very 


conservative preliminary design figures. Condenser volumes for 
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the standard condenser and the deep-sea condenser are nearly 
equal. These are presented in Figure VI-1l. 

Pumping power and pump weight for the deep-sea condensers 
areanticipated to be relatively small but must be traded off 
with pressure vessel penetration weights. An estimated 2.6 kw 
of electrical power is required for the 40 kw propulsion plant 
previously mentioned. This figure represents a crude design 
with 2-tees, 2-90° elbows and 2 - 45 degree elbows and 2" 
Piping throughout. Detailed weight tradeoffs might show an 
advantage in using a two or three pass condenser to reduce the 
Size of the piping penetration in the hull. Such an arrange- 
Ment could also reduce header weights since smaller headers 
could be employed, but this would be at the expense of some 
increase in condenser size with an attendant increase in 
pressure vessel size and weight. Table VI-9 gives a summary 
for the pump for the 40 kw plant. It is estimated that with 


careful piping design, the power requirement might be reduced 


as much as 25 - 30%. 
Pump pressure rise, lbf/in? 16.4 
Electrical power required, kw 2.6 
Pump & motor efficiency, % - 42 
Pump weight, lbm 60 
Mass flow rate, lbm/sec 24 
Power 400~AC 
Condenser Single Pass 
Piping 2” ~Enroughout 


TABLE VI-9 
Submerged Condenser Circulating Pump Requirements 


for the Deep Sea Condenser (40 .kw plant) 
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Figure VI-12 is a composite graph of condenser weight and 
volume as a function of condensing surface area. By the use of 
mugure VI-12, one may compare the weights and volumes of the 
two types of condensers. Assuming now a design depth of 20,000 
ft, condenser pressure of 1.0 lbf/in*, and a 50 kw plant using 
the turbine presented in Chapter VII, one may predict the 
intermediate heat transfer loop minimum temperature necessary 
to reject the waste heat from the plant. Table VI-10 is a 
summary of these weights. 

While the assumed sea water temperature for the deep sea 
condenser is 75°F, the condenser size and weight could clearly 
be reduced if the same sea water temperature were selected as 
for the skin cooler. It should be recognized that the pressure 
vessel size and weight may also be reduced with the smaller 
deep sea condenser. 

While the deep-sea condenser is clearly a better choice in 
Overall weight and volume than a comparable skin cooler, it 
must be recognized that the deep sea condenser will be consid- 
erably more expensive to build and is more susceptible to leaks, 
which would be disastrous at extreme depths. 

The single greatest unknown in the skin cooler is the 
scale heat transfer coefficient and the effect of any preser- 
vative which might be applied to the outside surface. Experi- 
mental work will be necessary to prove the feasibility of the 


Skin cooler for this application. 
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ASSUMED CONDITIONS 


Sea water depth, ft 20,000 
Plant nominal power level, kw io On) 
Turbine efficiency, % T2320 
Overall plant efficiency, % 36.0 
Assumed pressure vessel diameter, ft 4.0 


Available pressure vessel surface for heat transfer, ft* 40.3 


Sea water temperature, °F 40.0 
Heat rejection rate, Btu/hr eye eee tN ia 
Condenser pressure, lbf/in? | 120 
Condenser steam inlet temperature, °F 243.5 
Regenerator effectiveness, 2% 70. 
Regenerator pressure drop, lbf/in?” Gea? 
Condensation rate, lbm/hr 29727 
Oxygen flow rate, lbm/hr 4.13 


Skin Cooler - 4" Cu-Ni and +" Cu Fins - Arrangement B of 


Figure VI-3 
Pressure vessel material HY-130 Steel 


Intermediate heat transfer loop minimum temperature, °F 86.2. 


Condenser Surface area, ft? 70.5 
Condenser volume, ft? 3.9 
Condenser weight (including foundations), lbm 2680.0 


Skin cooler weight (including piping & fw and sw pumps), lbm 888.0 
Total weight (less pressure vessel), lbm* TAGE 0 
Pumping power required, kw 235 


*Weight of fresh water coolant not included 


Deep Sea Condenser 


_Assumed design coolant temperature, °F 75.0 
Condenser volume, ft? D4 
Condenser weight (including piping & foundations), lbm 460.0 
Pump weight, lbm 66.0 
Total weight (less pressure vessel), lbm 526.0 
Pumping power required, kw 32 


TABLE VI-10 
Comparison of Skin Cooler with Standard Condenser and Deep- 


Sea Condenser Designs at 20,000 ft 
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Removal and Disposal of Non-Condensables 


The non-condensables may be removed from the condenser by 
a standard vacuum pump with the expenditure of a small amount 
of power. For a 50 kw plant it is estimated that the steam 
flow rate which is uncondensed for the deep-sea condenser of 
Table VI-10 is 2.28 lbm/hr. With a 4.13 lbm/hr oxygen flow 
rate (representing the total non-condensable flow rate), a 
vacuum pump capable of removing 1.69 standard cubic feet per 
minute would be sufficient. The startup condition of large 
amounts of excess hydrogen must also be met. Here one must 
assume a time required for hydrogen rich operation until 
sufficient diluent steam can be admitted into the combustion 
chamber. If a ten second startup time is assumed, and if the 
pilot section of the porous plug burner (see Appendix B) is 
used in startup, a period of about 3.5 minutes would be required 
to clear the condenser of hydrogen, assuming the presence of no 
other non-condensables. To permit somewhat longer startup 
times, it is recommended that two vacuum pumps be installed, 
Sach rated for 2.0 SCFM. A single pump withdrawing 1.69 SCFM 
from a oon REER at 1.0 lbf/in’* pressure and discharging to 
16.4 lbf/in? would use approximately .625 kw at a combined 
motor and pump efficiency of 50%. The estimated weight for 
the two pumps and motors is 20 lbm. A single pump would be 
Capable of the hogging operation. 

The sensing of the amount of non-condensables in the 
Condenser would best be accomplished at the vacuum pump dis- 
Ehearge. Instruments for detecting oxygen using the paramagnetic 
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Principle could be used for this application (15). Equipment 
designed originally for flue gas analysis (and acting on the 
Gatalytic principle) could be used to detect the hydrogen 
concentration. This would first, however, require dilution 
with a predetermined amount of oxygen, Since the equipment is 
limited to a volumetric concentration of combustible gas of 
about 15%. With this known amount of diluent oxygen and the 
oxygen concentration sensed from the paramagnetic detector, one 
could then determine the hydrogen concentration. Effluent from 
the hydrogen and oxygen detectors could then be redirected to 
the condenser. The two gas measurement systems are available 
in a single commercial instrument approximately 18" x 12" x 18" 
Which operates on 115 VAC, 60-hertz, drawing 75 watts of power. 
This instrument would, however, require some modification for 
the proposed application. Its time response is about 30 seconds. 
Elimination of the non-condensables is accomplished in two 
steps. The first step is a porous plug recombiner very similar 
to that described by Moore (16) for the recombination of radio- 
lytically decomposed hydrogen and oxygen from the steam condenser 
of a Beiiing water reactor power plant. The device has proved 
_to be reliable for this application. The principles of operation 
of the porous plug burner are discussed in Appendix B. If the 
plant is run very slightly oxygen rich, any excess hydrogen may 
be recombined at the porous plug burner. A second porous plug, 
downstream of the burner plug, may be used to condense water 
vapor from the burned mixture, which may then be routed to the 


Condenser via a steam trap. A receiver downstream of the 


OZ 


condensing plug may act as a reservoir for uncondensed gases, 
chiefly argon with some oxygen and uncondensed water vapor. 

The content of this receiver might also be sensed for the 
mresence of hydrogen and the contents could be dumped to the 
condenser if the hydrogen concentration becomes excessive for 
compression of the mixture to 5000 lbf/in?, the storage pressure 
selected for on board storage of non-condensables. For the 
hydrogen rich startup condition, a suitable amount of oxygen 
could be bled into the line downstream of the vacuum pump and 
upstream of the porous plug burner. 

While overboard discharge of the non-condensables is a 
possible alternative, on board storage at high pressure has 
distinct advantages. The possible discharge of hydrogen at 
the surface could be hazardous. The amount of the non-conden- 
sables is small. It is estimated that a 7.5 inch ID sphere of 
HY-130 could house all the argon from a 1000 kwh mission at 
36% overall thermal efficiency if the argon were compressed to 
5000 1lbf/in?. Such a sphere would weigh about 7.5 pounds. Two 
such Spheres would probably be sufficient to house all the non- 
condensables. The electrical power to compress the non~-conden- 
. Sables is estimated to be about .2kw, assuming an overall 
compressor and motor efficiency of 20%. Intercooling would have 
to be provided, necessitating some sort of auxiliary cooling 
water system, which could easily be provided for this component 
and others. A sea water/auxiliary cooling water (fresh water) 


heat exchanger would be necessary. 
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Table VI-11 is a tabulation of required equipment and 
estimated weights for the non-condensable removal and storage 


system. Figure VI-13 is a schematic diagram of the system. 


Component Estimated Wt. Estimated Vol. 
lbm fe 
Vacuum pumps (2) and motors 3.0 0.08 


Porous plug burner, condenser 

and non-condensable gas 

receiver Lay 40 0230 
mydrogen-Oxygen Detector Ise, PAS) 


Non-condensable gas compressor 
(allowance for space and weight) 20.0 330 


Non-condensable gas storage 
Spheres (2) 14.8 «30 


Total 67.8 a23 
TABLE Vi-il 


Non-Condensable Removal and Storage System 


The condensate pump and feed water pumps have not been 
designed in detail, but the best solution appears to be a 
combined three stage condensate and feed pump. The first 
stage of the pump would serve a dual purpose. It would pump 
the product water to the product water tank as well as supply 
the second stage of the pump with sufficient net positive 
Suction head. The product water mass flow rate, of course, 
would be nearly identical to the reactant flow rate. For 
the very small volumetric flow rate, 0.172 gpm discharge at 


design power level, a gear pump would probably prove to be 
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most satisfactory for the proposed application. Table ViI-12 
lists the requirements and estimated weight and volume for such 


a pump for the 50 kw power plant described in Table VI-10. 


Type Gear -3 stage 
Discharge pressure, lbf/in? 800. 

Design power flow rate (discharge) gpm O21 72 

125% design power flow rate (discharge), gpm 0.214 

Drive Variable speed AC motor 


Electrical power required at design flow rate, kw G25 
Weight, lbm 2020 


Volume, ft? 0.3 


TABLE VI-12 


Condensate/Feed Pump 


A schematic diagram of the condensate/product water/feed 


system is shown below. 


TO REGENERATOR 
AND COMBUSTION |PRODUCT WATER 
CHAMBER STORAGE 


SAS a 


PRODUCT ILL VALVE 
WATER ¢ 

CONTROLT HOTWELL 
VALVE 


QQ _1Q) conpensate 


BONDENSATE/ |O © QO 
FEED PUMP 


FIGURE Vi-14 


Condensate/Product Water/Feed System 
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The opening of the product water regulating valve could be 
coupled to the flow rates of one of the reactants, arena 
hydrogen, and closed if the hotwell level drops below a 
specified level. 

An additional system is required for cooling of electronic 
components, motors, and other heat loads within the propulsion 
plant pressure vessel. A list of such equipment is given in 


Table VI-13. 


Lube oil cooler (turbine) 
Cycloconverter (if utilized for AC power transmission) 


Inverters (for startup and emergency battery to AC power 
conversion) 


Porous plug recombiner 
Propulsion plant pressure vessel atmosphere cooling 
Non-condensable compressor intercooling 


Hydrogen and oxygen heaters (if provided) 


TABLE VI-13 


Heat Loads for Auxiliary Cooling Water System 


No detailed study of the auxiliary cooling water system 
Was conducted. Rather a weight/space/power allowance 1S pro— 
vided. A weight of 150 lbm for all the heat transfer equipment 
and pump, including a sea water/fresh water heat exchanger,is 
believed to be conservative. Similarly, a volume allowance of 
1 ft*® and a power requirement of .2 kw for the system seems 
adequate. The atmosphere inside the pressure vessel could be 
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hydrogen saturated with water vapor or oxygen saturated with 
water vapor. Electrical components must be designed to operate 
in the environment selected. If a high speed generator is 
selected for the transmission system, which appears most pro- 


mising, hydrogen would be preferable in reducing windage losses. 


Condenser Off-Design Performance 

@hrough use of Figure VI-4 or VI-5 one may construct a 
graph of the estimated off-design performance for a condenser 
of a specific design condensing rate and pressure. Such a 
graph is shown in Figure VI-15 for the design conditions of 
Table VI-10. While extrapolation of Figure VI-4 or Figure VI-5 
will predict a zero condenser pressure at a zero condensation 
rate, a better estimate of the existing condenser pressure at 
zero steam flow is given in reference 10, estimating a pressure 
of slightly eee than 0.5 Ibf/in- for “standard air veqection 


equipment" and a coolant inlet temperature of 75°F. 
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CHAPTER VII 


Expander Selection and Design 


There are two basic choices which can be considered for 
the expansion of high pressure high temperature steam. These 
are the reciprocating engine and the turbine. Until recently 
very little work has been done in perfecting the reciprocating 
steam engine, particularly at the low power levels of interest 
in this thesis. With the current interest in low emission 
engines for automotive purposes, however, reciprocating steam 
engines have been reevaluated in the light of new technology. 
Steam engines in the past have operated on saturated steam, 
mach in addition to a liberal supply of lubricating oil, pro- 
vided sufficient lubrication. At the high superheated steam 
temperatures which, as shown in Chapter IV, give high overall 
thermal efficiencies, a different type of lubrication must be 
used. Syniuta (1) describes the development of an automotive 
reciprocating steam engine, which with the proper dry lubricant, 
permits Pee rat ion at stéam temperatunesyas high as 1400°F. The 
dry lubricant serves the purpose of replacing the condensed 
steam, which provides a large percentage of the lubrication in 
a conventional reciprocating steam engine, and eliminates the 
contaminating nature of lubricating oil which would foul conden- 
ser surfaces and reduce heat transfer at low condenser pressures. 
Detailed information on these engines is not available, but 


they present a promising alternative to the turbine. 
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A second type of reciprocating engine is one which burns 
the hydrogen and oxygen reactants in the cylinder. Such an 
engine has been built (2) which operates on excess hydrogen 
to reduce combustion temperatures to acceptable levels. Two 
significant problems were noted with this engine. The oxygen 
injector experienced rapid deterioration, and at the low ex- 
haust pressures of a space environment, there was a significant 
loss of lubricating oil into the engine exhaust. While steam 
Might replace the excess hydrogen as the diluent, the loss of 
lubricating oil could probably be rectified only through sub- 
stitution of a solid lubricant. The low condenser pressures 
shown in Chapter IV to give high overall thermal efficiencies 
would again result in the collection of oil on the condenser 
heat transfer surfaces. It is doubtful that the problem of the 
oxygen injectors could be easily overcome for an engine with a 
design lifetime of about 5000 hours. 

The turbine appeared, then, to be the best selection at 
present for the role of the expander. Baljé (3), in a classical 
paper on the optimization of turbine design, showed that at low 
Specific speeds, the single stage axial impulse turbine is 
Superior to the Terry turbine. He also showed that at even 
lower specific speeds, the partial admission axial impulse 
turbine and single disk, multiple-stage axial impulse re-entry 
turbine were superior to the full admission turbine. It is in 
the range of specific speeds where the re-entry turbine is most 
efficient, that the expander of this thesis must operate. This 


1S particularly true of the 50 kw power level. Of course, if 
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large enthalpy drops are to be taken across the turbine, a 
single disk re-entry turbine, even though it may have as many 
as four stages, will not be as efficient as a multidisk 
machine. In normal steam turbine practice, multistage 
machines are the rule and partial admission is common only in 
the first stage. In the proposed application, however, 
simplicity of design requires a small machine of not more than 
three or four disks. Radial inflow turbines are a possible 
alternative to the single stage axial impulse turbine but 
could not be easily staged. 

It is fortunate that the re-entry turbine design permits 
the utilization of a single disk for more than one-stage. Thus 
a machine of three disks could incorporate as many as twelve 
Stages, at four stages per disk. The repeated routing of 
steam through the blading of a single wheel, however, requires 
the addition of re-entry ducting. Thus the simplicity of the 
rotor design is negated to a certain extent by the addition of 
meme Cumbersome re-entry ducts. The re-entry ducts may become 
quite large in some configurations, effectively doubling the 
diameter of the machine, compared to a normal machine of 
equivalent rotor diameter (4). 

Figure VII-l is a schematic diagram of a two stage re- 


entry turbine with the crossover duct geometry. 
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FPEGURE ViIL=1 
Two Stage Re-Entry Turbine - Schematic Diagram 


(Crossover Duct Geometry) 


A primary limitation of the axial impulse re-entry turbine 
is interstage and tip leakage. The leakage phenomenon is 
described in detail in Appendix C. When the ratio of axial and 
tip clearance to blade height is small, the machines are quite 
efficient, but the efficiency deteriorates rapidly with in- 
creasing clearance ratios. A clearance of .004-.005" with a 
tip to tip rotor diameter of 8" gives very good efficiencies, 
but increasing the clearances bv a factor of 3 may result ina 


drop in efficiency by as much as 15 efficiency points at the 
igs Be 





same pressure ratio and entering steam conditions. With the 
high inlet steam temperatures, thermal growth of first stage 
rotor blading will require careful design of component clear- 
ances. It is doubtful then that clearances less than .005" 
can be maintained. This becomes a problem when one investi- 
gates the required blade heights for a re-entry turbine to be 
used as the first disk at the 50 kw power level. If a four 
stage re-entry turbine is used for the application, one must 
sacrifice efficiency over a considerable range of the available 
total pressure ratio. Here multidisk staging is advantageous. 
By using, for the first disk, only a. two stage re-entry turbine 
with small blade heights and large relative clearances, one 
may limit the low efficiency portion of the turbine to just a 
small portion of the total enthalpy drop available. 

Another problem associated with the proposed application 
is material limitations. Of the many high temperature, high- 
strength superalloys developed, few have sufficient strength 
at 1600-1750°F to withstand rotative speeds in excess of 80,000 
mae «6THis, of course, is coupled with disk diameter, which 
practically could not be much less than about three inches, and 
blade height, which could not be less than 0.1 inches. Ideally, 
One might set some limitation such as 0.25" on blade height, 
but even for re-entry turbines this is too large for the first 
Stage of a 50kw turbine operating with high pressure, high 
temperature steam. 

Another advantage of the re-entry turbine is also related 


tO Materials. In the normal partial admission turbine, the 


ms) 


blade temperature reached in the first stage is nearly the 
stagnation temperature corresponding to the relative inlet 
conditions to the blades. In the re-entry turbine, the blading 
experiences this same stagnation temperature in the first stage 
arc of admission but the same blading experiences lower stag- 
nation temperatures in succeeding stages on the same wheel, 
which cover a much wider sector of the arc of admission than 
the first stage. Even a two stage re-entry turbine operating 
at moderate stage pressure ratios is able to reduce the average 
stagnation temperature experienced by the rotor blades, allowing 
a higher rotative speed than for a comparable turbine with a 
simple partial admission first stage. 

The speed to be selected for a small turbine is not only 
a function of the best speed for highest turbine efficiency, 
but also the best speed for the transmission system selected. 
It is here that the electrical transmission system is advanta- 
geous since the weight of an electrical generator decreases 
with rotative speed for the same power produced. With solid 
rotor machines, rotative speeds of 24,000, 36,000 (5) and even 
93,500 RPM (6) may be reached. Of course, at these high 
rotative speeds and for an extended time between overhauls of 
4000-5000 hours, bearing reliability may become significant. 
For this application, ball bearings would probably result in 
the lowest friction and hence the least power loss. Syniuta 
(7) reports bearing de ratangs Gf as highwas 1— 2.5 x 10° “may 
now be utilized for high speed ball bearing aainleae TOnEP where 


qs = (bore diameter in mm) x (RPM). 
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Of course bearing life is of vital importance, and the 
more conservative the bearing qd rating, the longer the ex- 
pected lifetime. Calculations of a 30 mm bore diameter, 

36,000 RPM (d_ = 1.08 x 10°) design for the closed Brayton 
cycle engine (8) estimated a system time between overhaul of 
3100 hours. This system used two angular contact (combined 
journal and thrust) ball bearings with gas cooling and oil-mist 
lubrication. Total bearing loss for this machine was about .4 
kw. The operating conditions for these bearings are similar to 
that for the proposed cycle of this thesis with the exception 
of the presence of water vapor and the greater thrust which 
Must be absorbed from a steam turbine. Subsequent models of 
the Brayton unit have used gas bearings with satisfactory re- 
sults (9). Gas bearings using high temperature superheated 
steam would not be suitable for the semi-closed Rankine cycle. 

While the desired transmission system speed is of import- 
ance in selecting the turbine speed, one finds in practice that 
in order to achieve acceptable turbine efficiencies at low 
power levels, (50 kw) the turbine speed is of primary importance. 

It is here that the electrical ep enenseeon system is most 
compatible, in that it is also adaptable to high rotative speeds. 
Geared drives for propulsion purposes at the 50 kw power level, 
however, must provide gear ratios in excess of 100:1 if at all 
acceptable turbine efficiencies are to be realized. This high 
reduction ratio is at considerable expense in gear transmission 
efficiency. The efficiency is still high, however, but other 


factors, discussed in Chapter VIII, make the gear transmission 


unattractive. 
is 


It would then seem best to operate the turbine in its 
mest efficient speed range, consonant with material limitations 
and at a speed compatible with alternator frequencies which are 
most useful. A speed of 24,000 RPM corresponds to a 400-hertz 
2-pole machine, and iS a desirable speed for low power (50 kw) 
applications. Other possible rotational speeds are 36,000 RPM 


(4-pole, 1200-hertz) (5) and 48,000 RPM (4-pole,1600-hertz). 


faepine Optimization 

The methods described by Baljé (3) and amplified by 
Linhardt and Silvern (10) and Linhardt (11), were used to 
develop turbine optimization computer programs for full admiss- 
ion, partial admission, and multistage re-entry axial impulse 
turbines. Baljé and Binsley (12), describe the turbine optimi- 
mieeaa, technique used in this thesis, the pattern search. The 
particular method used is that of Hooke and Jeeves (13). The 
equations for the re-entry turbine are developed in detail in 
Appendix C. “it was aeeseere to modify the leakage analysis 
of Linhardt (11) considerably. lLinhardt neglected leakage on 
the exit side of the rotor both in the radial and tangential 
directions. The possible leakage paths in a re-entry turbine 


ere Shown in Figure VII-2. 
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PIGURE Vii-—2 
Leakage Paths in a Multistage Axial Impulse Re-Entry Turbine 


(3-stage) (crossover duct geometry) 


The leakage paths are described below: 
- tangential leakage in the direction of rotor 


bean AL 
rotation from one stage to the succeeding stage, 


upstream of the rotor 


eZU 


tan 2 


tan 3 


tan 4 


tip 1 


eip 2 


fcip 3 


Diad,u 


L 


rad,d 7 


tangential leakage in the direction opposite to 


rotor rotation from the first stage to the last 


stage, upstream at the rotor 


tangential leakage downstream of the rotor in 
the direction of rotor rotation from the exhaust 
of one stage to the exhaust of the succeeding 
stage 

tangential leakage in the direction opposite 
rotor rotation from the exhaust of the first 
stage to the exhaust of the last stage 

normal tip leakage, over the end of the blade 
from the pressure to the suction side of the 
blade 

leakage across the tips of the biades in the 
direction of rotor rotation from one stage to 
the exhaust of the succeeding stage 

leakage across the tips of the biades in the 
direction opposite rotor rotation from the first 
Stage to the exhaust of the last stage 

radial leakage into the rotor cavity upstream 

of the rotor 

radial leakage into the rotor cavity downstream 


of the rotor 


Tip leakage for the cases of interest from paths 2 and 3 


was found to be small in comparison to radial and tangential 


leakage and was not included in the computer optimization. An 


analysis of this leakage, however, is developed in Appendix C. 
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In addition to the leakage'paths described above, there 
exists another possible leakage path. This leakage results 
from blade movement into the exhaust duct region of the 
succeeding stage while in the blade passage there still remains 
working fluid. This effectively bypasses the nozzle of the 
succeeding stage. This type of leakage is termed dynamic 
leakage (10). In the computer analysis it was assumed that no 
work was performed by this leakage. Baljé (14) has assumed 
that only half of this leakage does no work. An interesting 
aspect of the dynamic leakage is the fact that it is speed 
dependent, i.e., it increases with speed. 

The modifications to the leakage analysis of Linhardt are 
presented in detail in Appendix C. Basically they include 
tangential and radial leakage downstream of the rotor, omitted 
by Linhardt, and assess leakage per stage as a percentage of 
the total steam flow rate for the multistage disk. lLinhardt 
assessed the leakage fraction on the basis of the nozzle flow 
for a particular stage, which overestimates the effect of 
leakage. While Linhardt's assessment on the basis of stage 
nozzle flow rates may be adequate for small clearance ratios, 
it is not satisfactory for the applications of this thesis. 

Reference 15 reports the experimental testing of a two 
Stage re-entry turbine operating at a pressure ratio of 300:1. 
Observed efficiencies for this turbine fell short of the 
efficiencies predicted by methods developed by Baljé (16). 

The principal causes of this discrepancy were deemed to be 


leakage by path L in excess of that anticipated and 


tan, L 


Ze 


dynamic leakage, unknown prior to that time. This reportedly 
was experimentally substantiated by a "static" test in which 
the rotor was held stationary and the upstream axial leakage 
gap reduced to zero. The remaining leakage was nearly equiva- 
lent to the value predicted for tip leakage. The difference 

in the leakage observed in the static test and the total leak- 
age which would produce the observed efficiency was evaluated 
as dynamic leakage, which agreed with a theoretical assessment 
of dynamic leakage. This assumed that no work was done by the 
dynamic leakage. The conclusions of this experiment are 
questionable, in that the conditions in the turbine change when 
no work is extracted from the rotor. Nevertheless, the success 
of this turbine is indicative of the promise of re-entry tur- 
bines for low power output and high enthalpy drop. It also 
indicated that the method of analysis developed by Baljé, 
Linhardt and Silvern can give reasonable approximations of 


expected turbine efficiency after further refinement. 


Presentation of Data 
The computer program employing the equations set forth 
in Appendix C was used to evaluate turbine efficiency over a 
Wide range of conditions. The results of this study are 
shown in Figures VII-3, VII-4, VII-5 and VII-6. The conditions 


assumed are listed below: 


16.0°. 


Q 
II 


Pr.= 1.97, constant for all stages 
Re. = .003 
23 





SSINISUNL TVIXV AYLNSSY GNV NOISSIWGY 

TVILYVd ‘NOISSINGY T1N4 ¥O4S WVYOVIG 

AONSZIOISS5 °N SJLISOMWOD SATA SYNIIS 
~N-Q43dS 91419 3ds , 


















001 10 
UT TET 
TALE EEE A I 

a €09'°O2 4Yy 

BPuLS cu - 
eee ee: 
wees : —— WILad Of yy = 
a Ne nA 

SA ae 1o0-G | OV 7) 
WLLL ee TLL | fee 
, aa og 4 
|) ee 

WA aoe 

ij p% AMLNDSRY S9vlS + a 

d ZA on ~< 








MA) 
cage Bee 20s 
BOWLS -S NVHLINAIDISIS = 
set BYOW Dow1S-+ _1§000'9 = O/€s 2 A/t%s aS/'s 


- S00 °9 =T/,76 
AMLNA-3Y BOWS € Py dee 
Lo‘l-4a POWLsS 
Wild AQANB-3a 











2 
DILVLS OL 7 


424) 





S 3 ONG 1 Sere GN 

SACV 1g 30 SYSEWNN TIVWS YOsS WVHS VIC 

AON SID ae aaa SNigant Rae eeice 7 eee ail} | 
a °N - Gl eieve,S JlalOads 











| ai Lard 26046 
esau 2 Cjrg =A/ig OF [TN] 
se St i a 
il 
Nba ean AL 
Sm 
YA. Re 
= Ze IO 
We MI ATT LV 
aA Las 

HEL Ze BOVIS 2 ttt 


mee vis OL VIO! Bai) eves I 


125 





SSONVdV4a 10 VIG Gye er sae Oe Sea iN 


SLVIGSNYSLNI HLIM SSNISYNL AYMLNSSY 
HOS WVYH9VIG AONZIOISSS—-N a 


ONE oases YlabJads 
To I TH 


HET Eo EP] pp | 
Lb'l=%q avDvwASsS 
lik SLO000'O 2A/€s = d/ts = O/'¢ Coos 
| S1090°0 =A7y y 
li Pe 4 
AdIN3IJ3d 39VLS vb —-— Gil 
i AWING3Y 39VLSe ----- Bai 
lr AWULN3dY 2 2 
Ita i ae 


HELL — aed Tl WEEE 


OO| 


” yA ONSIOISd3 
BaeeviS Ol IVWLOl aNident 








S3qdv1d JO YSEWNN TIVWS GNV SSONVYVST10 
S9YV1 ATSLVYSCON-SSNISUYNL AYLNISY 
HOIWVYOVIG AONSIOISSS°N 9TTA SYHNO!IS 


>N - Gaard S Sao aes 
00! Ol | 














| 
Eee LT II 
9100'0=AUeEe aQ jes aig O | 
ERs Pitta 
Ba Nasi BON 1S 2 aa ae OZ 
Lifsansay govse | | OTT | Dm 
INfeeeirabkaerh MLC: 
Lr A (NTT 8 
| ve Fo 
ee 
@) 
= Oe —< 
UL wea | | 
5 al Oz 
HBR | 
Lal eee 
teeth calllpoleba “ 


Sil Ope Sebel Ul 


= 


SL vrs ‘OL 


UAE | 





~~ 
It 
© 
© 
[= 
> 


13 
ke = 1.25 
vy = 0.96 
where a, = nozzle angle measured from the tangential direction 


Pr.= the stage pressure ratio 


Ky = the disk friction coefficient (see Appendix C) 
Kp = pumping loss coefficient (see Appendix C) 
k = ratio of specific heats 


by = nozzle velocity coefficient 

In the evaluation of the re-entry turbine the ratio of 
mecor blade throat width, a*, to rotor tip to tip diameter, D, 
was varied as well as the ratios of tip (s3) and axial clear- 
ance (sj, So) to rotor diameter. The significant variation of 
turbine total to static efficiency at constant specific speed 
as a*/D, s,/D, S5/D, and s,/D are increased is apparent. Tur- 
bine total to static efficiency has been plotted against 


specific speed. 


1 
N Q% 
—=- —— (VII-1) 
Ss q 
H : 
ad 
= Specific speed in the notation of Baljé (3) 
Q, = stage exit volumetric flow rate, ft*/sec 
Hoa isentropic head, total to static, ft lbf/lbm 
N = RPM 


Baljé (14) in unpublished lecture notes, reports similar 
findings of rapid deterioration of efficiency as clearances 


and throat width are increased for a fixed rotcr diameter. 
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The effect of increasing the ratio of throat width to rotor 
diameter is to decrease the number of blades, thus increasing 
the filling and emptying losses. Filling and emptying losses 
are common to partial admission turbines and were first 
identified and evaluated by Stenning (17). 

Because of the emphasis on high pressure ratio turbines 
in references 10, 11, and 15, the effect of pressure ratio on 
turbine efficiency at constant specific speed was evaluated. 
Large pressure ratios invariably resulted in significantly 
lower efficiencies. 

For the case of stage pressure ratios greater than the 
eritical value, the general first stage leakage factor, %, may 


be represented by the following relationship: 





: Piste 
_ ™e,1 _ YRTi st 
=U UDCA VEL Se) 
, my To throat 
YRT 
1,0 
where 
My 174 general leakage mass flow rate for the first 
7 
stage, lbm/hr 
m, = the first stage nozzle mass flow rate, lbm/hr 
st. the static pressure in the first stage cavity, 
7 


assumed in this case equivalent to the total 
pressure upstream of the leakage area, lbf/ft?* 
Pio 7 the total pressure upstream of the first stage 
7 


mozzle, Lotro 


; 
| 





A, = the leakage area, it 
= 2 
S hroat = the nozzle throat area, ft 
T) pet the total temperature upstream of the leakage 
? 


area, assumed in this case to be equal with the 
exit static temperature of the first stage rotor, 
oR 
T1,0 = the total temperature upstream of the first 
stage nozzle, °R 
R = the gas constant for the gas in question 
It may be seen that the leakage factor is inversely 
proportionai to the stage pressure ratio. Hence, for increased 
pressure ratios at constant Ne one might expect the leakage 
factor to decrease. 
For subcritical pressure ratios, in a turbine of equal 
pressure ratios per stage, the same relationship holds, but 
Bnly because the pressure ratios for the first and second 


Stages are equal. 


K+1 
Pi,st "2 /2k [pucss, nae K 
‘ /RE ed Pi,st PB 


Bai de’ 1 at ee 


™  Pi,o*throat /2k /Pr sty ?/*_ Pi,st, 
RT, > ar 15,0 Pi ,0 
? 


Post Peer Piist 
l,st P1,0 


ei, , 1 


where the ratios are equal. IS 
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The results of the computer study show, however, that 
turbine efficiency, for a constant No decreases with pressure 
ratio. This effect is small at moderate increases in pressure 
ratio but becomes significant at stage pressure ratios as high 
as 4.5. 

In practice one finds that many factors combine to produce 
the degraded efficiency at constant specific speed as pressure 
ratio increases. The major effect, however, is that of reduc- 
tion of the arc of admission in the first stage. This can be 
understood if one considers an incompressible fluid in a four 
Stage turbine, for example. In this case the arcs of admission 
for all four stages are equal. With a compressible fluid, the 
first stage arc of admission becomes progressively smaller with 
respect to the other stages, increasing the filling and emptying 
losses in the first stage, as pressure ratio increases. Even 
as the total degree of admission is increased, as seen from 
Table VII-1, efficiency decreases. Because the effect of fill- 
ing and emptying losses becomes so significant at small arcs of 
admission, the optimum blade height becomes smaller, resulting 
in a compromise between degraded efficiency due to reduced 
nozzle arc length and increased leakage due to increased rela- 
tive clearances with respect to blade height. 

The values given in Table VII-1l are computed at a constant 
ratio of blade throat width to rotor diameter, a*/D and axial 
and tip clearance ratios s,/), s,/D, and s,/D. From Table VII-1, 
One may see the effect of decreased first stage nozzle arc 


length and reduced blade height, which combine to increase 
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Data Common to Both Pressure Ratios 


4 Stage Re-Entry Turbine 


k es 2) S,/D = .00125 Ky = .003 
> i wes Gli a*/D = .015 
Ss 

X = 0.6 

Si mS ile Vy = 0.96 
mop, s./D = .00125 
1 2 
Pressure Ratio 4054 80:1 
Ls 4.04 S204 
B, 23.46 22. 5.2 
h/D OS sill .0276 
aee/c,)° .078 .084 
Ue 8143 .7924 
u/c, SL620 SS 
O13 510. Gia 46/2° 
Nps BS) olSdl - 4938 
A/A* 1.0604 Leek 270 
M,/M._» 1.228/0.849 1.34/0.959 
Soc -016 ORs 
SD 002 0.0016 
Sp 03:0 01-0 
a,/D O2 128 .1019 
% admission Sou I! 99.5 
t/2a, 0.1487 0.1923 
Cr. tip 0.0668 9.0790 
Gr, Pad 0.0738 0.0795 
Cr. een 0.0418 0.0495 
Sr, dyn 0.0424 0.0462 


TABLE VII-1 


Effect of Pressure Ratio on Turbine Efficiency at Constant 


Specific Speed 
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leakage as pressure ratio is increased. Further increases in 
pressure ratios would result in additional losses when the 
relative Mach number of the fluid entering rotor, M2! exceeds 
unity. It is interesting to note the increase in the ratio of 
blade pitch to twice the first stage nozzle are length, t/2a,, 
as pressure ratio increases. This is the parameter which 
increases filling and emptying losses. 

It is noteworthy that, with increased leakage, one cannot 
express analytically the function for the optimum ratio of 
blade height to rotor diameter, h/D, which is necessary in the 
computer optimization, and is of interest in the effect of 
pressure ratio on leakage. Baljé (3) and Linhardt (11) obtained 
approximate relations for h/D by partial differentiation of the 
expression for total to static efficiency. One finds for large 
clearance ratios, however, that this requires the solution of a 
fifth order algebraic equation and that reduction of this 
equation to a lesser order is not possible by simplifying 
assumptions. For this reason, the quantity h/D was selected as 
the optimization variable in the pattern search optimization. 
Thus the dependence of the blade height on pressure ratio 
cannot be seen directly. 

It should be noted here that the efficiencies presented in 
Figures VII-3 through VII-6 are based on an equal pressure 
ratio split between stages. This assumption is based upon the 
investigation of Linhardt and Silvern (10), (15), which com- 
pared the cases of equal pressure ratio and equal isentropic 


head pec stage for a two stage re-entry turbine with an overall 


JESS) 





pressure ratio of 300:1. It was shown that over a wide range 
of specific speeds, equal pressure ratio staging was clearly 
preferable to equal isentropic heads in each HEE) in the 
upper range of specific speed the advantage was as much as ten 
efficiency points. This does not imply that the equal pressure 
ratio split is the best choice in all combinations of number 

of stages and overall pressure ratio, and it probably is not. 
It does, however, result in an estimate of achievable turbine 
efficiency. An advantage of the assumption of equal pressure 
ratios is the relative ease with which the turbine performance 
may be evaluated. The equations for equal pressure ratio are 
themselves quite complicated (see Appendix C). The inclusion 
of the element of varying stage pressure ratios renders evalua- 
tion considerably more difficult, especially when considering 

a family of turbines over a wide range of specific speeds. 

Wong (4), however, used the concept of unequal pressure 
ratios in each stage, probably to reduce the leakage from the 
first stage rotor cavity to the second stage and thereby to 
Paice overall leakage. This may be seen from the leakage 
equations derived in Appendix C. Leakage from the first stage 
rotor cavity to the second stage in a turbine with stage 
pressure ratios greater than the critical is truly determined 
by the total pressure ahead of the leakage gap and not the 
interstage pressure ratio, as might be inferred from equation 
VII-2. Wong's turbine is of considerable interest in that it 
is the only four stage re-entry found in the literature to have 


actually been built and tested. The svecifications for Wong's 


turbine are given in Table VII-2. 
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Number of Stages 


Meror tip to tip diameter, D, in 


Blade height, h, in 


Average nozzle angle, Oo 


Blade angle, Bo 


Overall pressure ratio 


Design speed, RPM 


Flow rate (air), 


lb/sec 


Blade pitch, t, in 


Stage pressure ratios: 


Stage l 
Stage 2 
Stage 3 
Stage 4 


Inlet conditions, 


Tip and axial clearance, in 


N 


> 


H 


ad’ 


Pressure, lbf£/in* 


Temperature, 


et Ibf/lbm 


mee £t°/sec 


Predicted total to static efficiency (Wong) 


Observed total to static efficiency 


rN 


TABLE ViiI—=2 
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42125 
OA Jers) 
224° 
Se.” 
55.56 
24,000 
0.0044 


0.10 


4.332 
2.045 
2505 
20726 


On ZZ 
Ono 19 


0.432 


Data on a 4 Stage Re-Entry Turbine with Unequal Pressure Ratios 





Using the computer program developed for re-entry turbine 
analysis and the overall pressure ratio, clearance ratios, and 
nozzle angle for Wong's turbine, a total to static efficiency 
of .341 was predicted for a turbine with equal pressure ratios. 
The observed total to static efficiency of Wong's turbine was 
0.432. Several conclusions may be drawn from this analysis. 
The first is that the computer program, based on the design 
method of Baljé, Linhardt and Silvern and modified by the leak- 
age equations in Appendix C, is conservative. A possible 
second conclusion is that a re-entry turbine with unequal 
pressure ratios may very well be superior to one with equal 
pressure ratios. Since Reynolds number effects were not 
included in the equations of reference C, it is concluded that 
Reynolds number effects in Wong's turbine, for which he allowed 
in his analysis, were overshadowed by leakage effects. Although 
Wong showed an appreciation of interstage leakage, it is 
apparent that he did not consider dynamic leakage and that he 
was probably unaware, at the time, of the parallel work of 
Baljé. The small blade height, 0.125" and blade pitch, 0.010", 
used in Wong's turbine are important in that they show that a 
turbine with very small blades and blade passages can, in fact, 


be built and give reasonable efficiencies. 


Reynolds Number Effects 

The effect of Reynolds numbers on turbine efficiency has 
been investigated by Baljé (18), Bullock (19), Holeski and 
Stewart (20) and many others. A rapid deterioration of turbine 


efficiency occurs with separation from the back of the rotor 
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blades. Before separation occurs there are increasing viscous 
losses below some critical Reynolds number, that of the transi- 
tion from turbulent to laminar flow. The difficulty is to 
predict the flow conditions at which this will occur. One may 
use a common method described by Horlock (23) involving loss 
factors to predict the effect of Reynolds numbers. If one 
considers that all the losses are Reynolds number dependent, 
the loss factor, &, may be replaced by 1-np 


S° 


mee lost brs 


—,)* = — (VII-4) 
Re 1I-NnLs 


where the primed values indicate the condition of interest. 
Here it is implied that at Reynolds numbers above some value, 
10°, based upon the hydraulic diameter of the blade passage 
and blade inlet relative velocity, viscous losses are constant. 
The exponent, #, in equation VII-4 varies, depending upon the 
investigator (19). Another common value for the exponent is #. 
Baljé based his rotor Reynolds number on blade chord 
‘length and relative inlet velocity to the rotor blade and dis- 
tinguished between rotor and stator hoseeen He also established 
a machine Reynolds number based on blade tip speed and tip to 
tip rotor diameter for correlating Reynolds number effects at 
Various specific speeds and pressure ratios. Baljé found that 
the critical Reynolds number, based upon chord length, at which 
the transition from turbulent to laminar flow occurs, varies 
with specific speed. The eritical Reynolds number did fall, 


Mewever, in the range of 10° to 10°. 
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Bullock (19) discussed the various attempts at analyzing 
Reynolds number effects on turbomachinery ene en@ted ste: wide 
variance of opinion which exists in predicting Reynolds number 
losses. Bullock also discusses the causes of the transition 
from turbulent to laminar flow peculiar to turbomachines. Some 
of these are the thickness of the trailing edge of the blade, 
the smoothness of the blade finish, and the existence of time 
unsteady effects which promote turbulence and delay transition 
and separation. 

Holeski and Stewart (20) studied a large sample of full 
admission axial flow turbines for Reynolds numbers which varied 
from 10* to 2 x 10®. Again the Reynolds number was based on a 
different characteristic length, in this case blade height. 

The velocity used by Holeski and Stewart is different from 
others, in that they used the axial velocity, assuming that the 
variance in nozzle angles and blade angles is small. Viscosity 


mee based on stage inlet conditions. 


Repoleski and Stewart ux 


where = rotor radius at mean blade height, ft. 

Holeski and Stewart attempted many different methods to 
correlate Reynolds number losses and reported that a plot of 
Stator throat area against loss parameter was most successful. 

The loss parameter defined by Holeski and Stewart is quite 
complicated and is given in reference 21. It is evident, 


Nowever, that Holeski and Stewart's evaluation does not cover 


the case of off-design performance. 
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Stewart does, however, give data on the variation of loss 
factor with Reynolds number. The data for Reynolds numbers in 
the region of interest for the 50 kw application of this 
thesis are that of Wong and Nusbaum (22). Wong and Nusbaum 
evaluated the performance of a full admission turbine using 
the rotor of the 4-stage re-entry turbine of reference 4. 

Again a different definition of Reynolds number was used, based 
on blade height, relative Velocity to the rotor sand jctarer 
blades, and a viscosity based on mean static conditions. 

In the interest of comparison of the various reported loss 
correlations with the data of Wong and Nusbaum, Figure VII-3 
was plotted against the loss factor ratio of equation VII-4. 
Equation VII-4 was also plotted as well as a similar relation 
with an exponent of 3. It is of interest that in the Reynolds 
mumber range from 10* to 10°, (in the notation of Holeski and 
Stewart), Wong's data consists of three curves. At low 
Reynolds numbers, Wong and Nusbaum's turbine operated with 
three distinct operating curves, attributed to varying degrees 
Of boundary layer separation in the rotor. Furthermore, the 
data for Wong's turbine was based on a blade to jet speed ratio 
(u/c.) Ge .225, Lar from the normal ratio of abcout.0..5 for 
axial impulse turbine. In the lower range of Reynolds numbers, 
efficiencies at higher ratios of blade to jet speed were signi- 
Brcantly poorer. 

It is now possible to define a comparable Reynolds number 


for partial admission using Holeski and Stewart's notation. 
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m 
Re partial _ a/2m u 
admission 


(VII-6) 


where a= the nozzle arc length of the are stage. 

Mere it is noted that Holeski and Stewart"s notation, in 
addition to ignoring the effect of nozzle and blade angle, also 
omits a constant factor of 27 in the definition of Reynolds 
number (21). 

BO base the turbine design O@& this thesis on thesdacaper 
Wong and Nusbaum is probably somewhat overoptimistic, primarily 
because of the degree of reaction on which their data is based. 
Wong and Nusbaum's data does, however, give an indication of 
Reynolds number effects in the region of interest for this 
thesis. It was therefore decided to utilize a compromise of 
the + and + power curves with a critical Reynolds number of 
mexelO°. With proper design, it is evident from reference 20 
that turbine efficiency is not significantly affected above 


this Reynolds number. 


Turbine Design 
Using the design curve for Reynolds number from Figure 


fta-7, and the design curves of Figures VII~3 through VII-6, 


several 50 kw turbine designs were attempted. The number of 
possible designs using full admission, partial admission, and 
re-entry turbines of various pressure ratio and number of 
Stages is very large and not easily subjected to computerized 
Optimization. The thermodynamic studies, discussed in Chapter 


iV, determined that as high a turbine inlet temoerature and 
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pressure as possible would produce the highest overall thermal 
efficiency. Limitations on the porous plug burner restricted 
turbine inlet pressure to about 600 lbf/in*. Based upon studies 
of available turbine blade materials, a turbine inlet tempera- 
ture of 1750°F was selected. It was also shown in the thermo- 
dynamic studies that as low a condenser pressure as possible 
meatd be best. With 1.0 lbf/in*? as a practical lower limit of 
condenser pressure, disregarding for the time its effect on 
condenser size and volume, the turbine exit pressure was 
Meeablished as 1.14 lbf/in* allowing for a 0.14 lbf/in? 
pressure drop across a regenerator with an effectiveness of 

me 70. 

For the design of the 50 kw Aes eee a re-entry turbine 
was selected for the first disk because of the effect on blade 
stagnation temperature described earlier in this chapter. Two 
and four stage re~entry turbines were compared at a stage 
pressure ratio of 1.97,slightly greater than the critical value. 
fee tour stage re-entry turbine yielded poor efficiency over a 
considerable range of the available pressure ratio of the tur- 
bine, and for the limiting blade height of 0.10", resulted in 
a rotative speed in excess of 90,000 RPM. A two stage re-entry 
turbine was therefore selected for the first disk, at an over- 
all pressure ratio of 3.88, an RPM of 24,000, and a total to 
Seacic efficiency of 0.364. 

The second disk stage also used a stage pressure ratio of 
1.97, but because now larger disks could be used, this disk was 


designed as a 4 stage re-entry turbine. The final disk, a two 
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stage re-entry turbine, was designed for an overall pressure 

of 9.05. This higher pressure ratio in the last two stages 
(third disk) was selected on the basis of a desire to minimize 
Reynolds number losses in this disk due to the high specific 
volume of steam at low pressures. As will be shown later, this 
mgorce was fortuitous in that it permitted the turbine to 
operate with all nozzles choked down to power levels below 40% 
of the design power level. In each of the stages it was assumed 
iat the entire dynamic head was lost in the re-entry ducts and 
from disk to disk. It was found that a full admission stage or 
stages in place of the third disk would have to operate at 
Significantly higher rotative speeds and a consequently smaller 
eeeameter and blade height to be at all efficient. 

An overload nozzle for the 125% design power level was 
not included in this design but could be easily added. If the 
turbine should be used with a geared drive, it would be 
necessary to provide an additional re-entry disk of at least 
two stages for backing purposes. The windage and blade pumping 
losses of this disk would detract slightly from the efficiency 
predicted under ahead operation. 

Figure VII-8 is a diagram of the rotor of the proposed 
turbine, and Table VII-3 lists the detailed characteristics of 
the 3 disks. The predicted overall turbine efficiency, neglect- 
ing bearing losses and labyrinth leakage, for this turbine is 
Sraghtly greater than 0.73. This value is believed to be 
conservative, particularly in view of the Reynolds number 


POrrelation used. 
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STI 


Overall Turbine and eycre Data 


Turbine inlet pressure, lbf/in* 600. Turbine exit temperature, Slee) - 
Turbine inlet temperature, °F 50. Regenerator exit temperature (steam) 243.5°F 
ferbine exit pressure, Lb£f/in- Tel4 Nozzle angle, Qos (all stages) Ie) 
Condenser pressure, Ibf/in* La (60) Total number of stages 8 
Turbine mass flow rate, lbm/hr IRS) Total number of disks 3 
Cverall turbine pressure ratio 526.0 RPM 24,0010 
(Disk and Blade Material- 

Overall total to static Ng? bo ne 

efficiency (neglecting 0734 

labyrinth seal and bearing losses) Nozzle velocity coefficient, Wx 0.96 


Actual turbine power developed-kw-53.5 (neglecting bearing and labyrinth seal losses) 


Data on Individual Disks 


——_—_ 


* Type of disk 


Number of stages 

Total pressure ratio, Pr, 

Stage pressure ratio, Pr; 

Rotor diameter (tip to tip), D, in 
Blade height, h, in 

Axial clearance, S1/Soe in 

Tip clearance, Sor in 
Stage admission, % 


Total admission, % 


Recor Olade pitch, t, 2 


Rotor blade chord length, C, in 


Rotor blade angle, Ba 


Stage Reynolds number, a 
(after Holeski and a 
Stewart) 

ate, ee 9 


Nozzle area ratio, A/A* 


Length of first stage are of 
admission, aye Th}9. 
Ipig 


Maximum blade stress, lbf/in? 


Maximum disk stress, (tapered disk) 


Lbf/in?¢ 


50 kw, 
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TABLE VII-3 


3 Disk - 8 - Stage Re-Entry Turbine 


Disk 3 
Re-entry 
2 
905 
5.02 
O72 
YAO S 
.0054/.0054 
-0054 

S57 9.2 
54.7 
Onions. 
Orne? 
28.4° 


Geox 0 7 AosenOe 


1.345/0.774 
16 LA 


Drea 


0.764 
15300 
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The design of the turbine caSing and crossover (re-entry) 
ducts was not attempted in this thesis and will require some 
innovation. The re-entry ducts should be as small as possible 
to reduce casing weight and size, yet a desire to recover the 
leaving loss from a previous stage would suggest large ducts. 
Clearly, some compromise solution considering these tradeoffs 
and the three disk geometry must be reached. In the estimation 
of turbine weight and volume, a conservative estimate of casing 
weight and size has been used. 

’ The 500 kw turbine was designed in much the same SV EIETE as 
the 50 kw turbine. In the case of the 500 kw turbine, however, 
a four stage re-entry turbine could be used on the first disk. 
The same rotative speed, 24,000 RPM was selected in order to 
keep the size of the machine small and because of its compati- 
bility with a 400 hertz, 2-pole alternator. It is of interest 
that aside from having one additional disk, the 500 kw turbine 
Size is nearly identical to that of the 50 kw turbine. The 
Second disk is a two stage re-entry turbine and the third and 
fourth disks are full admission disks. Dynamic head between 
all Beces, including full admission stages, was considered 
lost. The overall predicted turbine efficiency was 84.3%, 
neglecting diaphragm labyrinth seal and bearing losses. It is 
estimated that a 1 kw bearing loss for this machine and a .66 kw 
bearing loss for the 50 kw machine is reasonable for two 
combined thrust-radial ball bearings based on reference 8. It 
Must be shown, however, that such bearings are capable of main- 
taining the small axial clearances required by this machine. 
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Figures VII-9 and Table VII-4 give a detailed description of 
the 500 kw turbine. A regenerator was not utilized with the 
500 kw turbine, since its high turbine efficiency results in 

a reasonably low turbine exit temperature. Estimates of weight 
and volume and dimensions for both the 50 and 500 kw turbine 
are given in Table VII-5. 

The small blade throat widths associated with the re-entry 
portions of the turbines are a consequence of the optimum blade 
angle and improved performance gained from maximizing the 
number of blades. This reduces the filling and emptying losses 
of Stenning (17), common to all partial admission turbines. 
While these throat widths are quite small, Wong's turbine (4) 
is an indication that they can be fabricated. The small blade 
and blade passages will require very close attention to such 
items as trailing edge thickness and surface finish. The disks 
have been designed with a central hole and would be keyed to 
the shaft. A taper of 2:1 from hub to disk periphery has been 
assumed in the disk stress calculations. Blade and disk 
stresses are well within the present state-of-the-art Of high 
temperature alloys. The blades and disk would be manufactured 
from a single casting. Electrochemical machining methods would 
probably be best for forming the blades. Critical speed 
calculations were not performed on either turbine. Rather, a 
shaft size of approximately 1.5" was assumed for disk stress 


calculations and bearing loss approximations. 
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Qverall Turbine and Cycle Data 
Turbine inlet pressure, lbf/in* 
Turbine inlet temperature, °F 
Turbine exit pressure, lbf/in? 
Condenser pressure, lbf/in* 
Turbine mass flow rate iDeyuy Aegis 
Overall turbine pressure ratio 

No regenerator 

Overall total to static efficiency 
(neglecting labyrinth seal and 
bearing losses) 

Data_on Individual Disks 

Type of Disk 

Number of stages 

Total pressure ratio, Pr, 

Stage pressure ratio, Pr. 

Rotor -atancter (E1p to tip), D, in 
Blade height, h, in 

Axial clearance, S1/85 

Tip clearance, s 


3 


Stage admission % 


Total admission % 

Rotor blade pitch, t, in 

Rotor blade chord length, C, in 
Rotor blade angle, B 5 


Nozzle angle, «a 


2 
Stage Reynolds number 
20m a 
Ay u Tm" 


(after Holeski and Stewart) 
Mo/M 9 
NOZZLG arca FratLro 


Length of first stage arc of 
admission, in 
ps 


Maximum blade stress, lbf/in? 


Maximum disk stress (tapered 
disk), lb£/in? 


600. 
L730. 
Ihe 

1. 
2550. 


600. 


0 


0 


Turbine exit temperature, °F 360. 
Total number of disks 4 
Total number of stages 8 
RPM | 24,000 


Actual turbine power developed 
(neglecting bearing and labyrinth 
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ies, 

9.95 
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005 
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TABLE VII-4 


500 kw - 4 Disk Turbine 


seal losses) 3500 
Disk and blade material IN713LC 
Disk 2 Disk 3 Disk 4 
Re-entry Pula ude 
2 it a 
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50 kw turbine 


Casing diameter without return ducts, in 20 
Maximum diameter, in NSO) 
Overall length, in 18.0 
Casing length, in 13.0 
Turbine weight, lbm 105.0 
Foundation weight, lbm 30.0 
Total weight, lbm 135 0 
Volume, ft? 1.8 


500 kw turbine 


Casing diameter without return ducts, in 14.0 
Maximum diameter, in 20.0 
| Overall length, in 20.0 
Casing length, in L520 
Turbine weight, lbm 186.5 
Foundation weight, lbm B70 
Total weight, lbm | 223.5 
Volume, ft? | Da 


TABLE Vit=5 : 


Turbine Weight, Volume and Dimensions 


" P50 





Control, Method of Operation and Offdesign Performance 
With the small length of the arcs of adnission of the 


first stage of the first disk in both the 50 and 500 kw turbines, 
throttle control with perhaps a single nozzle control valve 

seems the best method of speed regulation. The nozzle control 
valve would be sized to close SERS portion Of thesarc of 
admission at the cruise condition, in this case assumed to be 

40% of design power. In both the 50 and 500 kw turbine, 
sufficient unadmitted arc exists in the first disk for the 
inclusion of an overload nozzle to achieve the 1.25% design 

power specification. 

Both constant speed and propeller load off-design charact- 
eristics have been considered. The constant speed application 
would require an electrical transmission system and a turbine 
governor for speed regulation. The propeller load application 
would require a suitable geared transmission for speed 
reduction to propeller RPM. 

The ‘off-design study was conducted only at 41% of the 
design power level and for the 50 kw turbine. At this point, 

a 41% Pere tnrorciea condition and a condition with the single 
nozzle control valve being fully shut and the throttle valve 
fully open were considered. The results are shown in Figure 
VII-10. The propeller load curves were estimated from the 
constant speed curves modified by data on propeller load 


turbine characteristics at optimum nozzle openings (24). 
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In the calculation of off-design performance it was noted 
that the last stage nozzle was still choked at the 41% power 
condition, which implies that all nozzles remained choked. 
Data from Lee (25) were used to estimate losses due to under 
and over expanded nozzles. In the 41% power nozzle condition, 
the reduction of the length of the arc of admission in the 
first stage negated, to a great extent, gains made by the 
elimination of throttling. At lower power levels, Reynolds 
number losses were also found to increase in importance. A 
graph of stage exit pressure for the design condition, the 
41% power throttled condition, and the 41% power condition 
with the single nozzle control valve fully shut and throttle 


valve fully open are shown in Figure VII-1l. 
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List OF SsYMBOLS 


nozzle arc length 

blade throat width 

nozzle area ratio 

absolute velocity 

absolute velocity leaving the stator 
absolute velocity leaving the rotor 
spouting velocity, V2g. Hog 

blade chord length 


BetLOt= tip tO tip dlamerer, Et or in 


Sig 
D Hoa 


¥Q 





specific diameter- 


Conversion factor 32.2 ft lbm/1lbf sec? 


isentropic head, ft lbf/lbm 

ratio of specific heats, assumed to be 1.25 for steam 
dusk friction coefficient 

blade pumping loss coefficient 

Mach number of fluid leaving the stator 

relative approach Mach number of fluid entering the 
rotor blade | 

mass flow rate, lbm/hr 

RPM 

specific speed (Baljé) NYO/H,* 

pressure, #bfL7im-morelbe7 ce 

pressure ratio 

exit volumetric flow rate from a disk, ft*/sec 
Reynolds number 
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r - mean rotor blade radius 


m 
Sy - axial clearance upstream at the rotor 
sm - axial clearance downstream of the rotor 
S, - tip clearance 
c ae Dlade pitch 
a. - nozzle angle measured from the tangential direction, ° 
B. - blade angle measured from the tangential direction, ° 
Cr, - leakage factor 
Np-g ~ Total to static turbine efficiency 
u - dynamic viscosity, lbm/hr ft 
3 - loss factor 

- blade density factor 

- velocity coefficient 

SUBSCRIPTS 
dyn - refers to dynamic leakage 
il - refers to stage conditions 
L - refers to leakage : 
m - refers to a general stage of a multistage re-entry 
turbine 

N - refers to nozzle 
O - refers to total conditions 
rad - refers to radial leakage 
SC - refers to scavenging losses 
st - refers to static conditions 
mo -— refers to tangential leakage 
mip - refers to tip leakage 
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CHAPTER VIEL 


Auxiliary/Emergency Power and Transmigsion System 


The auxiliary/emergency power systems and the transmission 
system for a deep submersible power system are closely inter- 
related. While it is not the intention in this thesis to 
optimize such a system, the selection of a satisfactory system 
is necessary to prove feasibility of the propulsion plant. 

Such a transmission system must meet acceptable standards of 
reliability, efficiency, weight and volume. 

The geared transmission and shaft seal, from the aspect of 
high efficiency, appears to be very promising despite the lower 
turbine efficiency of propeller load turbine operation compared 
with constant RPM turbine operation. For example, consider a 
240 RPM propeller speed, which would entail a reduction of 
100:1 in speed from a design turbine speed of 24,000 RPM. 
According to reference 1, a spur gear set with this reduction 
ratio would exhibit about 94% transmission efficiency. This 
figure would probably be somewhat reduced by the high rotational 
speed of the high speed pinion. If one assumes a 90% trans- 
mission efficiency for such a spur gear set for the 24,000 RPM 
application and neglects, for the moment, thrust bearing and 
Sealing losses, a combined turbine/transmission efficiency of 
66.1% is predicted for a 50 kw (turbine shaft) plant at design 
power. At the 50% design power level, this would result ina 
combined turbine/transmission efficiency of 53.5%. The propul- 


sor efficiency here has not been considered. 
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Electrical transmission systems are penalized in power 
losses in the generator, the speed regulation device and the 
motor. Typical electrical transmission efficiencies for a 
50 kw generator power supply are presented in Table VIII-1 (2) 
and compared with a possible reduction gear drive. The geared 
drive does not include thrust bearing or shaft seal losses. 

The selection of a transmission system is not quite so 
simple as the results of Table VIII-1 would suggest. The 
reduction gear system, while seemingly attractive, incurs many 
penalties, the least of which are seal leakage and required 
memping power to control it. Consider, for example, a l 3" 
shaft penetration at a 20,000 foot depth. The thrust load in 
such a penetration would be about 16,000 1bf. For the 8,000 
foot depth, a penetration of the same size would result ina 
Barust load of about 6,300 ibf. A Kingsbury thrust bearing to 
carry this load would have a diameter of about 1 foot at the 
20,000 foot depth and eight inches at the 8,000 Loot depth 
Considerable structural support would have to be supplied at 
the thrust bearing, increasing overall weight. Additional 
Seructural eiahe would also be necessary at the shaft penetra- 
eion. Practically, the shaft should be as small as possible to 
reduce the thrust and buckling load and still support the 
torque. A proper evaluation of such a system would require 
determination of shaft length, diameter, torque and buckling 
load to minimize the thrust load. Shaft seals have not yet 
been developed for either the 8,000 foot depth or the 20,000 
foot depth. While a shaft seal may be possible at the 8,000 
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foot depth, it is doubtful that one can be designed for 
20,000 feet with an acceptable leakage. 

Another problem associated with shaft sealing is the 
flexing of the hull of the propulsion plant pressure vessel 
under compression and its effect on misalignment of the 
reduction gear. A second, somewhat less serious problem is the 
difficulty of internal arrangement posed by a relatively large 
condenser, a turbine, and a reduction gear and thrust bearing. 

Another major problem associated with the reduction gear 
transmission is the need to supply auxiliary electrical power. 
If the turbine does not operate at constant speed, some other 
power conversion device must be provided for electrical power. 
Possible candidates are batteries, fuel cells, or a small tur- 
bine. In addition to added weight and volume, this would 
increase system complexity considerably. Table VIII-2 gives 
a summary of auxiliary electrical loads. 

A significant and perhaps crucial disadvantage of the 
geared transmission system is the need for more than one 
Meerpatsion device for the typical deep submersible. To meet 
the needs of hovering, reasonably exact positioning and manuver- 
ability, there is a requirement for more than one propulsor. 

It is here that the electrical transmission system is most 
versatile. 

The case of the reduction gear transmission system has not 
been studied in detail. It is apparent, however, that associated 
with it are technological problems which are considerably be- 
yond the state-of-the-art. At the 8,000 foot depth, the 
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problems for a 50 kw plant could probably be overcome but limits 
on shaft size may rule out higher power levels. At the 20,000 
foot depth, the use of such a transmission system is considered 
very unlikely. For this reason, an electrical transmission 


system has been selected for purposes of this thesis. 


Electrical Transmission 

At the high rotational speeds compatible with the turbine, 
solid rotor, brushless AC generators are the most likely 
candidates. Two such generators are the Lundell-Rice alter- 
nator and the homopolar inductor alternator. 

The Lundell-Rice alternator has been proposed for use with 
the Brayton cycle for a 50 kw submersible power plant (3) with 
an efficiency of 94% (neglecting bearing losses) at a rotational 
speed of 24,000 RPM. Homopolar alternators have been built at 
rated speeds of 93,000 RPM and 10 kw power output (5). A 
description of a 4-pole, 400 hertz- 12 kw homopolar alternator 
Operating at 12,000 RPM is given by Corcoran (6). The 
Mietred efficiency of this machine was 91.7%, neglecting 
bearing losses. It appears that the efficiencies, weights and 
volumes and associated electrical components (voltage regulator 
and exciter) for the two types of machines are comparable. 

The Lundell-Rice rotor has been tested for windage losses 
(7) with rotors as large as 12" in diameter for 24,000 RPM. 
This rotor is about the size required for a 500 kw alternator. 
Thus the Rice alternator may be used at a rotational speed of 
24,000 RPM for all rated power levels from 50-500 kw. The 
rotor diameter (6") for a 50 kw machine is given by Balukjian 
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and Rackley (3). The characteristics of the generator and its 
voltage regulator and exciter system are discussed by Corcoran 
(6) and given in detail by Dunn (8). Table VIII-3 gives 
estimated weights and volumes of two such alternators, rated 
at 50 and 500 kw. The use of hydrogen in the high speed 
alternators may reduce the windage losses and slightly improve 
the estimated efficiency of 94%. 

The problems which have been encountered with the compon- 
ents of the electrical transmission systems listed in Table 
VIII-1 are discussed by Bloomquist (9). Based upon the 
experience of the Naval Ship Research and Development Labora- 
tory, one would conclude that AC submerged motors and encapsul- 
ated power conversion systems would be the best approach at 
present. It is apparent that placement of inverters and other 
such power conditioning equipment in the sea water environment 
will require a considerable development effort before reliable 
equipment becomes available. 

While the results of Table VIII-1 indicate that the AC 
hydraulic propulsion system is the best choice with respect to 
Mieteiency, it is apparent that the AC Srellowanie ver system 
is also very promising. The AC cycloconverter system is 
slightly lighter but dynamic braking would be required for 
this system while it would not be necessary for the electric- 
hydraulic system. The AC hydraulic submersible system operates 
at constant voltage and therefore achieves Porat ively Gocd ort 


design performance. 
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Alternator Power Level 


50 kw 500 kw 
Lundell-Rice Alternator! 
Rotor diameter, in 6 12 
Machine diameter (frame), in D2 ay 24 
Machine frame length, in 27 ir 
Shaft diameter, in 135 2.0 
Estimated volume, ft? 0.93 Bede) 
Estimated weight, including foundations 230.0 720.00 


meltage Regulator-Exciter and Speed Regulator 
Weight and space allowance lbm/ft? 207 0. > 6071.5 


Inverters, Rectifiers and Startup/Emergency Batteries 


DC-AC inverters (2), (Mounted internally) 


power rating kw 5 50 
volume (est) 0.3 2.0 
weight (est) 10 67 
Startup/femergency batteries (AgZn), kwh 5 50 
Weight (in air @ 47 wh/lbm)?, lbm 106 530 
Volume (@ 3.5 wh/cu in)*®, ft? 0.83 diye 


Power conditioning equipment weight 
and volume allowance, lbm/ft? 10/.5 30/1.5 


‘Weights and Volumes estimated from drawing in reference 12 
and reference 7. 


*Reference 13. 


3estimated from reference 14. 


TABLE VIII-3 


Alternator and Related Components and Extraneous 


Electrical Equipment 
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The cycloconverter system is discussed by Thomas (10) 
and Barnes (11). The primary advantage of tre cycloconverter 
system is the significantly improved off-design performance 
of frequency control of induction motor speed over voltage 
control. With the assumed power profile of a 40% design power 
level for nearly 85% of the total mission duration, the cyclo- 
converter appears to be an interesting and promising alternative 
to conventional AC induction motor speed control. 

Koegel (15) reports that a cycloconverter system has been 
used with a 10 KVA Lundell-Rice generator operating at 100,000 
RPM with frequency reduction to 400 hertz. The application is 
a portable power supply for the U. S. Army driven by a gas 
turbine. The equipment was supplied by Lear-Siegler and 
General Electric. It is apparent, however, that little 
development work has been performed in the application of 
cycloconverters to propulsion systems for deep submersibles. 

The two systems, AC hydraulic and AC cycloconverter, 
cannot be adequately compared without an in-depth study involving 
cost, present and potential availability for deep submergence 
applications, weight and volume, and reliability. Such a 
study is beyond the scope of this thesis. 

The AC submersible hydraulic system was selected because 
of its current availability and high efficiency. The slight 
weight advantage of the cycloconverter system (about 50 lbm) 
is considered insignificant. The AC submersible hydraulic 
System is listed in Table VIII-4. It includes two 400-hertz 
AC submerged induction motors and hydraulic pumps and motors, 
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Equipment Wt (lbm) BEficiency 100% 
power/50% power 


Submersible Motors (2-30 HP) 
16-24 poles 
36, 400-hz 800 S07. 63 


Hydraulic Pumps 
(8) IL SiS) -88/. 86 


Hydraulic Motors 


2-15 HP 
6-10 HP 135 .86/.86 
Total 1250 .68/.61 


Data from reference 2 


TABLE VIII-4 


AC Propulsion Transmission System 


169 





which would be coupled to the propulsor(s). 

Other power equipment may be necessary to provide 
specialized power to various vehicle equipment. Startup and 
emergency batteries, located outside the pressure hull, will 


require inhull inverters for conversion of DC to AC power. 
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CHAPTER IX 


Reactant Supply and Storage System 


The reactant supply system and storage system is the 


largest and heaviest system comprising the semi-closed Rankine 


cycle propulsion plant. In all configurations this system 


must be encapsulated in hard pressure vessels, with considerable 


penalty in weight. Hydrogen-oxygen reactant supply and storage 


systems for deep submergence applications have been analyzed in 


detail by various contractors in studies for the U.S. Navy. 


Reference 1 reports the results of such a study which was 


conducted in conjunction with an analysis of a closed Brayton 


Cycle propulsion plant for deep submersibles. 


Several reactant supply systems are possible for the semi- 


closed Rankine cycle propulsion plant. Four possible configura- 


mons are listed below: 


(a) 


(b) 


(c) 


storage of hydrogen and oxygen as gases at high 

pressure with injection into a combustion chamber 
utilizing gaseous reactants 

storage of hydrogen and oxygen as subcritical cryogenic 
liquids, pumping them to pressures sufficiently in 
excess of combustion chamber pressure, and injecting 
them into a combustion chamber at cryogenic temperatures 
storage of hydrogen and oxygen as subcritical cryogenic 
liquids, pumping them to pressures sufficiently in 
excess of combustion chamber pressure, and heating 


them in a heat excnanger of an auxiliary cooling water 


~~ 
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system prior to eee into the combustion chamber 
as gases 

(d) storage of cryogenic hydrogen and oxygen at super- 

critical pressures, sufficiently in excess of combustion 
chamber pressure, and injection with or without prior 
heating. 

The first alternative has been found to be considerably 
heavier than the others at a design operating depth of 20,000 
feet. Reference 1 reports about a 7000 pound increase in the 
overall plant weight of the closed Brayton cycle propulsion 
plant with the high pressure gas storage system over that using 
a subcritical cryogenic liguid system. For lesser depths than 
20,000 feet, one can expect cryogenic storage to remain 
preferable to high pressure gas storage (2). Two advantages of 
the high pressure gas storage scheme, however, are its simpli- 
city and the elimination of the product water storage pressure 
vessel, in that product water can be pumped into either the 
hydrogen'‘or oxygen tank or both. This, however, must be 
accomplished by pumping against the existing pressure in the 
tank, Ortien’ represents a parasitic load. A considerable amount 
of unused reactant remains in the tanks when the tank pressure 
drops to the feed pressure for the combustion chamber. This is 
common as well to the other alternatives, pParcrcularly to the 
Supercritical system. 

The second alternative, subcritical storage, requires 
adequate insulation to reduce boil~-off losses and is consider- 
ably more complex than the high pressure gas system. This is 
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true for all the cryogenic storage systems. To achieve the 
relatively high combustion chamber pressures necessary for 
good overall plant thermal efficiency requires pumping the 
reactants as liquids. Hydrogen and oxygen become supercritical 
fluids when they are pumped above about 200 and 750 lbf/in’, 
respectively. This pumping is another parasitic load and 
requires additional reactant to accomplish it for the same 
mission energy requirement (e.g. 1000 kwh). Cryogenic pumps 
require cooldown prior to operation,and sufficient suction 
head to operate satisfactorily. This second alternative is 
particularly applicable to the liquid-propellant rocket-type 
combustion chamber. 

The third system is identical to the second system except 
for heat exchangers downstream of the cryogenic pumps, which 
heat the cold fluid to normal temperatures. If the selected 
combustion chamber utilizes gaseous reactants, increased 
efficiency can be realized by heating the cryogenic reactants 
by sea water (through a suitably designed intermediate heat 
transfer loop) prior to their coming in contact with the com- 
bustion tide. This increases the available enthalpy of the 
Beactants. 

The fourth alternative allows heat to be admitted to the 
cryogenic storage tanks, either by an electrical heating 
element or by returning cryogens heated by sea water to "heating 
coils" in the tanks. An inner pressure vessel is required in 
this application to contain the supercritical fluids at storage 


Pressures from 100-250 lbf/in? above combustion chamber pressures, 
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the pressure depending upon the type of combustion chamber and 
chamber pressure selected. The supercritical storage method 
then incurs a weight penalty due to the inner pressure vessel 
and requires slightly less insulation than the subcritical 
‘systems. Its advantage is the elimination of the cryogenic 
pumps. 

Using the electrical loads of Chapter VIII for a nominal 
50 kw, 1000 kwh endurance plant, based upon generator electrical 
load over and above all parasitic loads, a comparison was made 
of alternatives (b), (c) and (d). First, the electrical load 
was used to establish a reactant load required for the assumed 
mission profile for the 50 kw, 1000 kwh plant. These loads are 
listed in Table IX-l. A 94% generator efficiency and a total 
of 1 kw in bearing losses for both generator and turbine (3 
bearings) were assumed in computing the turbine power required 
for each of three conditions: full power (5 hours), 60% power 
(1.7 hours) and 40% power (35 hours). The period of operation 
at power levels above full power is included in the full power 
operation. The power levels apply here to all loads other than 
the parasitic loads, deo@ sa; ae eeailsien power, manned sphere 
power requirements, and external auxiliary equipment. The 
parasitic load includes reactant pumps, which of course are not 
necessary for the supercritical storage system. This is 
corrected in the reactant load necessary for the supercritical 
configuration. Parasitic loads actually would be slightly 
increased over the values in Chapter VIII, since the turbine 


Power required is ahout 10% greater than that used for the 50 kw 
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Loads Power Level 


Parasitic Loads 


kw 
Condensate/feed pump 0.55 
SW circ pump 3.20 
Vacuum pumos 0.625 


Auxiliary cooling water pump 0.20 
Non-condensable compressor 0220 
Reactant supply pumps 1.14 


H,-0, detection equipment 


and propulsion/reactant 





supply controls 02575 
Bearing losses 1.00 
Total parasitic load 7.49 
Propulsion Load 44.28 
Manned sphere power 2.44 


External Auxiliary Equipment 3.28 


Full Power 











Total Electrical Load 5. 49 
Non-parasitic Load 50.00 
Turbine Load = Gly. L0 
(b) 
Subcritical 
»torage—Pump 
H/05 


Reactant required, lbm* 2217/7 Ess 


*assumes a 98% combustion efficiency by weight 


tabi x= 


60% Power 40% Power 
kw kw 
0.33 O22 
ao. 20 3.20 
02625 0.625 
0.20 0.20 
0.20 C220 
0.684 0.456 
02575 02575 
1.00 La00 
Ge cl 6.48 
24.23 14.28 
2.44 244 
3e28 B29 
362.8 1 26.48 
30200 20.00 
39 150 28.20 
Alternative 
ic) (d) 
Suberitical Super- 
Storage-Pump & Critical 
Heat H,/O Storage 
He 
216/1726 2117 boo. 


Reactant Required for a 50 KW, 1000 KWH Propulsion Plant 


LL 





turbine design and used for computing the parasitic loads in 
Chapter VIII. This was ignored for the purposes of this 
analysis. 

For alternatives (b), (c), and (d), the amount of unused 
reactants at the end of the mission was determined by assuming 
meroo lbf/in* gas pressure at the end of the mission for 
@iternatives (b) and (c), and a 400°R, 700 lbf/in*® residual gas 
condition for alternative (d). It was also assumed that fil ding 
the cryogenic storage tanks to greater than a 95% capacity 
would not be possible since some space must be left to permit 
boil-off without inordinate pressure increases. Table IX-2 
gives a summary of the required reactant volume, pressure 
vessel sphere inner diameters, and pressure vessel weights for 
the cryogenic containment spheres. 

It was not considered necessary to calculate detailed 
weights for containment spheres. Balukjian (3) gives curves 
of pressure vessel thickness and weight for spherical pressure 
vessels of various materials for 8000 and 20,000 foot depths, 
but these are found to be considerably overconservative when 
Seared with detailed pressure vessel design figures in 
reference 4, a more thorough investigation of the Brayton cycle 
propulsion plant for deep submersibles. In the computation of 
Table IX-2, HY-180 Steel was used for purposes of comparison of 
the three reactant storage alternatives, despite the fact that 
it has not yet been proved acceptable for welding in heavy 
Sections (5). The pressure vessel weights, then, were derived 
Brom data in reference 4 for various spheres of similar size 


i 
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Total reactant required 


(b) 
Subcritical 


Storage-Pump 
BO a, 


264/1829 


(including unused reactant), 


lbm 


Required storage volume 
meeeost fill", ft° 


Pressure Vessel ID, in 
Vessel thickness, in 
Pressure vessel OD, in 
Pressure vessel weight 
Inner tank thickness, in 
Inner tank weight 


Insulation, various 
internal component 
weight (est.), lbm 


Cryogenic pumps and motors 


(est. ) 


Weights 
Subtotal H,/O,, lbm 


Total system weight, 
(including reactants) 
lbm 


M@ensity of H., at fill - 


2 


Memsity of O. at fill - 


Z 


63.07 27.0 


63.0/48.5? 
370 
65.6/50.5 
6000/2200 
O70. 1 
310/177 
150/130 


5/15 


6739/4351 


11090 


ApAQ ibm f° 


7le2 lbm/ete 


Alternative 


(ce) 


Subcritical 
Storage-Pump 
& Heat 

H5/0, 


252/1740 


60.0/25.6 


62.0748), 2- 
VA (6 
64.6/50.2 
5600/2180 
Ona 7/.0ral 
301/174 
150/130 


ube 3 es) 


6318/4239 


10547 


includes 2" thick multi-layer insulation 


3includes vacuum annulus and vapor shield 


PABEE Si) 


(d) 


Supercritical 


Storage 
Boy oh 


229/1848 


BAe O67 27.2 


59.0/48.0° 
1622 AO 
61.44/50.0 
4500/2160 
0.368/0.292 
LO Fey 537 
85/70 


5884/4615 


10499 


Cryogenic Storage Weights and Volumes for Various Storage 


Meeernatives for 20,000 Ft. Depth-HY 180 Steel Pressure Vessels 
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fa are only approximate. A safety factor of 1.5 based upon 
buckling has been used in the analysis. 

The supercritical storage system may utilize a vacuum 
annulus and vapor cooled shield for interruption of a consider- 
able amount of heat leakage through the pressure vessel walls. 
The subcritical system, storing of liquid in a subcooled state, 
will require multilayer insulation. Insulation techniques have 
been developed which permit storage of liquid hydrogen for as 
long as 24 hours in a Space environment without venting. This 
can be improved by initially filling the tanks with subcooled 
hydrogen and oxygen. 

It may be seen from Table IX-2 that the weights of the 
three systems are comparable. One may conclude then that 
other factors than weight will govern the particular system 
selected. At the 8000 foot depth, the subcritical systems are 
about 500 pounds lighter than the supercritical storage system. 

Because of the choice of the liquid propellant rocket type 
combustion chamber as the best candidate for immediate develop- 
Ment, reactant storage alternative (b) was selected. This 
system requires the development of very small reactant pumps 
Which can operate -with a minimum of fluctuation of discharge 
pressure. Large fluctuations in discharge pressure would not 
be acceptable for the combustion chamber. 

Small reactant pumps of the reciprocating or diaphragm 
variety could probably be devised. Tyree (6) reports that a 
particular reciprocating pump he devised could operate on as 


much as 50% gas by volume. These two tynves of pumps would not, 
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however, fulfill the requirement of minimum discharge pressure 
fluctuation for the liquid propellant rocket combustion chamber. 

A very small variable speed gear-type pump is proposed. 
Such a pump could be built in several stages (3-4) to very 
exacting specifications. Although the pumps would experience 
considerable leakage, this could be reduced by very small 
clearances. Pump cooldown could be conducted by bleeding the 
cryogens through the pump and burning the gases in the recom- 
biner of the non-condensable removal system (see Chapter VI). 
A diagram of this pump is shown in Figure IX-l. Table IxX-3 


gives a description of the pumps. 





FIGURE Ix-1 


Reactant Pump Stage 


iba 


eacaneter Oxygen 
Assumed reactant flow rate, lbm/hr 66.5 
Reactant density, lbm/ft? 71.25 
Volumetric flow rate, ft’/min 0.0156 
Gear OD, in O25 
Tooth length, in 0. 125 
mooth width, in OG 125 
Gear length, in 0.50 
rump, rpm 308 
Number of stages 3-4 
Assumed pump efficiency 0220 
Assumed motor efficiency 0.60 
mischarge pressure, 1bf/in’ 885. 
Electrical power required, kw 0.565 
Estimated pump and motor weight, lbm 15 
Estimated pump volume, ft? O)52 
TABLE IX=3 


Reactant Pumps for a Nominal 50 KW Plant 
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0.0314 
0.5 
Oe 5 


O25 


It would be necessary to provide these pumps with sufficient 
net positive suction head to prevent cavitation at the inlet of 
the first stage. This could be accomplished by pressurizing 
the cryogenic storage tank to a sufficient pressure, perhaps 
100 lbf/in*, thus subcooling the reactant. An arrangement to 
provide such subcooling is shown in Figure IX-2. A small 
volume of the cryogen is kept in a very small but separate 
dewar located inside the main cryogen storage tank. A small 
amount of gas boils off when an electrical heating coil is 
energized in the pressurizing dewar. This gas is sufficient 
to raise the pressure of the gas above the main cryogen liquid 
storage while the bulk cryogen temperature remains essentially 


constant. PRESSURIZER 


NORMAL FILL CONNECTION 





PRESSURIZING 


HYDROGEN OR 
OXYGEN STORAGE 
SPHERE 

_ _PRESSURIZER 588 
INITIAL FILL CONNECTION 


FIGURE” IX-2 


REACTANT 
PUMP 


Pressurizing Arrangement 
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The power drain of the pressurizer heaters was determined 
to be less than 2.5 kwh in a 1000 kwh mission. It may be 
necessary to devise some means of emergency fill of the 
pressurizer dewar should all the cryogen evaporate from the 
pressurizer dewar because of delayed plant startup. In this 
regard, it may be necessary to locate the pressurizing dewar 


deep in the tank, rather than at the top as shown in Figure IxX-2. 


Assessment 

The problems of cryogenic storage of hydrogen and oxygen 
for deep submergence applications are by no means completely 
solved. It appears, however, that with some developmental 
work in the area of transfer couplings and other areas, the 
System can be made safe and reliable. No detailed analysis of 
pump leakage rates for minimum anticipated clearances has been 
performed, but the assumed pump efficiencies are very low. If 
the reactant pumos for some reason should not prove satisfactory, 
it will be necessary to use the supercritical storage system 
and burn effectively gaseous reactants in the rocket type com- 
bustion chamber. The possible operation of a rocket type 


combustion chamber on gaseous reactants is discussed in Chapter V. 
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CHAPTERS 


Regenerator 


The regenerator is a relatively small piece of equipment 
which may be designed using standard techniques. The methods 
of Kays and London (1) have been used for the regenerator 
design. Regenerator effectiveness has been shown in Chapter IV 
to have only a small increase on overall plant efficiency 
compared to increases in turbine efficiency. Nevertheless, it 
was decided that if regenerator size could be kept within 
reason, perhaps less than one cubic foot in volume for the 50 kw 
plant, the maximum regenerator effectiveness attainable should 
be utilized. Regenerator effectiveness, according to Kays and 


London, is defined as 


Np = ~hot Cp hot ‘Th, in’ th,out? 2 ne Spc cout 7c, in) 
Mo da Ee asiae (ml 5/min*'h,in -¢,in 
(x51) 
where 
Tot = mass flow rate of the fluid on the hot side of the 
heat exchanger, lbm/hr 
mo = mass flow rate of the fluid on the cold side of 
the heat exchanger, lbm/hr 
fae) min = minimum product of mass flow rate and specific 
heat, hot or cold side 
e = specific heat, Btu/lbm°R 
Th = hot side temperature, °R 
we = cold side temperature, °R 
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Using a finned tube heat exchanger(core geometry figure 100 
of reference 1), a regenerator was designed for approximately 
half the pressure drop of 0.14 lbf/in* between turbine exhaust 
and condenser at a regenerator effectiveness of 0.90. This 
resulted in a regenerator of dimensions 8" x 8" frontal dimen- 
Sions, one foot in length. 

The regenerator could be placed in the exhaust ducting of 
the turbine, between the turbine and condenser or directly 
above the condenser. It has been assumed that about 14 feet 
of 8" diameter exhaust ducting would be adequate. 

A drawing of the regenerator is shown in Figure X-1l below. 


Specifications for the regenerator are given in Table X-l. 


STEAM IN 


a ve 






LEE 


iw” 





FEED OUT 
0.402" 
| FEED IN 
STEAM OUT 
FIGURE X-1l 


Regenerator Design 


Head) 





mass flow rate, lbm/hr 
Frontal area, ft? 
miboing, OD, in 

Mubing, ID, in 

Tube material 

Fin material 

Fin thickness, in 

Fin spacing (# per inch) 
Core length, ft 


Regenerator steam, lbf/in?* 
Side pressure drop 


Regenerator Exhaust Pressure, iof/in- 
Steam inlet temperature, °F 

Steam exit temperature, °F 

Feed inlet temperature, °F 

Feed exit temperature, oe 
Regenerator volume, ft? 

Turbine exhaust ducting volume, £t? 


Regenerator weight, lbm 


Exhaust ducting weight, lbm 


TABLE X-l1 


Regenerator Data 
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351 
0.450 
0.402 
0.382 
IN CASS 
cu 


O70 re 
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CHAPTER 


Encapsulation and Weight and Volume Summary 


The preceding chapters of this thesis have established 
information on weight and volume of the components of the 
proposed propulsion plant. The propulsor (propeller, Kort 
nozzle, ducted thruster, pump jet, etc) has not been specified 
Since the propulsor is dependent on the detailed mission and 
body form of the submersible. 

One may then, at this point, select power plant conditions 
and determine overall weights and volumes. These figures may 
then be compared with comparable power plants for the deep 
submergence mission. From the weight and volume data developed, 
one may add or subtract components for a more direct comparison 
with other propulsion plants. 

One item of importance is the selection of a suitable 


encapsulation material. Several materials are possible 


candidates: aluminum, titanium, and ultra high yield strength 
steels. It is not considered of great importance in this thesis 
to justify a selection of an encapsulation material. Rather 


it was considered necessary only to select a good candidate 
which would facilitate comparison with other proposed propulsion 
plants. The same is true, to some degree, with respect to the 
selection of a transmission system. Each proposed propulsion 
plant, whether it be fuel cells, Brayton cycle, Stirling cycle 
or the cycle proposed in this thesis has certain equipment in 
common. Certainly all encapsulated plants would utilize the 


0 


best material available, or which will probably become available 
in the near future. For this thesis, spheres have been selected 
for the shape of the encapsulating vessels. HY-180 steel was 
selected for the encapsulation material. It is realized that 
HY-180 steel could not be used for this application at present 
because of fabrication difficulties (1). As in Chapter X, a 
Barety factor of 1.5, based upon buckling load, has been 
assumed. 

The power plant conditions selected were those shown to 
effect the highest overall thermal efficiencies attainable. 
These conditions are shown in Table XI-l. A 98% combustion 
efficiency has been assumed with a 75°F sea water temperature. 
The efficiency of 27.1% for full power is based upon a turbine 
power of 62.4 kw to carry 50 kw electrical load other than 
propulsion related auxiliaries. The mission is as specified 
in Chapter III, 1000 kwh. 

The weights and volumes for a 500 kw plant have not been 
included Be caliee of the lack of a specified mission, without 
which one cannot determine reactant welenes. From previous 
chapters, one may, if desired, estimate weights and volumes 
of components for a 500 kw system or some intermediate power 
level between 50 and 500 kw. 

Tables 4l-2 and X13 give a summary of weights and volumes 
for two power plants. The first plant utilizes subcritical 
storage of the reactants and reactant pumps to deliver them 
at pressures sufficiently in excess of the 600 lbf/in? combus- 


tion chamber pressure for the liquid propellant rocket combustion 


oe 


SuUbGmiEtical SUpereritacal 


Reactant Storage Reactant 
System with Liquid Storage System 
Propellant Roc! et with Porous 


Combustion Chamber Plug Combus- 
tion Chamber 


miebine, rpm 24,000 
Turbine inlet temperature, °F 1, 720 
Turbine inlet pressure, lbf/in? 600 
Turbine efficiency, % 74.0 
Generator efficiency, 3% 94.0 
Condenser pressure, lbf/in? 10 
Condensate/feed pump discharge 

meesscure, lbf/in* 800 
Regenerator effectiveness, % 90.0 
Turbine exhaust temperature, °F 449.0 
Regenerator steam exit temperature, °F 145.6 
Regenerator and exhaust ducting 

pressure drop, lbf/in? Or4 
Regenerator feed inlet temperature, °F Eons? 
Regenerator feed exit temperature, °F 365.5 360.5 
Turbine mass flow rate @ full power 

lbm/hr 351 338 
Reactant mass flow rate @ full 

power, lbm/hr 9025 82.7 
Turbine load @ full power, kw 62.4 61.0 
Reactant required for 1000 kwh 

mission’, lbm 2150 1990 
Specific reactant consumption’, 

lbm/kwh 1.63 eS 
Overall thermal efficiency”, $% Pag fs 2071 


11000 kwh mission is over and above parasitic loads; assumes 
98% combustion efficiency by weight. 


“Based upon HHV of hydrogen @ standard conditions for 50 kw 
useful electrical power output. 


TABLE XI-1 


Final Plant Conditions 50 kw - 1000 kwh Plant 
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EG 


Component Volume Dimensions Weight 


fee L/W/H, £t (1bm) 
Condensing System 20,7000 £7 S000 - re 
Condenser foundations and piping 205 e972 A] eee 467 Z95 
SW circ pump and motor* —- Ou3 - 70 
Condensate feed pump and motor Case = Zao 
Non-condensable removal and storage 
system Se ~ 80 
Auxiliary cooling water system ee) ~ TSO 
Poe ene and "Regenerator 
feeb inewsand foundations dae oe de 5/1. 5 135 
Regenerator O25 206070 .667 1.0 85 
Exhaust ducting 0252 - 2 
Combustion Chamber 5 570.037 0205 50 
Picctrilcalmoys tem 
Lun@ell-Rice generator 0.93 1.06/1.06/1.06 210 
Inverters and rectifiers Ors = nko 
Power conditioning equipment Os = 10 
Propulsion plant controls 0:5 = . Ths 
Start-up/emergency batteries?’ 0.83 - 106 
Submersible motors and hydraulic 
pumps and motors? 5-0 test) ~ 1250 
Miscellaneous | 
Reactant pumps aS) = 30 
Piping Owd - 50 
Fluids (including sea water condenser 
coolant) - = 60 
Product water tankage (empty) EO =.9 ~ 190 
207 OC OMe S000 sre 
Total volume of Internal components 252 25 Ge 
Total volume of external components 4.13 4.13 
Total weight of internal components spel y) 1425 
Total weight of external components 1426 1426 
Total buoyancy of external components 
@ 64 lbm/ft*, lbf£ 264 264 
Propulsion plant pressure 
Vessel sphere, ID, in 48.0 48.0 
Sphere internal volume, Ft? 33.5 33.5 
Sohere packing factor (volumetric) O75 0.275 
Propulsion plant sphere weight, lbm 2,400 250 
Propulsion plant sphere buoyancy 
@ 64 lbm/ft’, lbf 2425 2265 
‘includes hotwell 
mounted externally 
TABLE XI-2 


Propulsion Plant Pressure Vessel and External Equipment - Weights and Volumes 


vet 


Subcritical Reactant Storage 
System with Liquid Propellant 
Rocket Combustion Chamber 


Reactant Storage system 


YD: 
Reactant 
usable, lbm 239/1910 
unusable, lbm 39/20 
ftoeal, Lom 278/1930 
Inner tank weight, 
11,/0. tanks, lbm 83207 192 
Miscellaneous weight 
insulation, etc., ibm £507130 
Inner tank diameter 
lea7O2, in 60/45.4 
2D 
Pressure vessel, ID 
14/051 Thy 64./47.4 
20,000 (£e78,7,00072E 
Pressure vessel thickness, H,/O5, in ie 377 LOe 0.65/0.48 
Pressure vessel weight, H5/O., lbm 6800/2220 2800/900 
Total storage vessel weights 
including cryogens, H/O5, lbm 7548/4472 3548/3152 
Buoyancy @ 64 lbf/ft?* 5780/2340 5500/2220 


Subcritical Reactant 
Storage System with 
Liquid Propellant 
Rocket Combustion 


Chamber 
Product Water Storage 
edie water volume, total, £t° oes 
Product water stored in propulsion 1s, 
plant pressure vessel, ft’ 
Product water stored in separate 2376 
product water sphere, ft? 
“roduct water pressure vessel ID, in 42.6 


20,000 ft/s, 0C0Ggrt 


Product water pressure vessel 


thickness, in 029270.43 
Product water pressure vessel, 

weight, lbm 1700/650 
Buoyancy at 64 lbf/ft?, lbf 1725/1590 


‘Operating supercritical pressure-700 lbf/in? 


TABLE Xl—=3 


Supercritical Reactant 
Storage System!’ with 
Porous Plug Combus- 
tion Chamber 


H,/0. 


221/769 
18/165 
239719 34 


1080/544 
85/70 


57/45.4 


60.6/49.0 
20,000 £t/8,000ft 
07 1. OS Ol C07 Ooo 
520072450. 22507 L000 


6604/4998 3654/3548 
4930/2620 4590/2390 


Supercritical Reactant 
Storage System with 
Porous Plug Combustion 
Chamber 


ogo 
OR 2, 


LALA 

41.0 

20,000 £t/8,000 £t 
0.88/0.41 


1550/600 
1510/1415 


Reactant Storage and Product Water System Weights and Volumes 


chamber. The second plant utilizes supercritical reactant 
storage at 700 lbf/in*? and the porous plug combustion chamber. 
Both plants were analyzed for depths of 8000 and 20,000 feet. 
A diagram of the internal arrangement in the propulsion plant 
pressure vessel is shown in Figure XIi-l. Lowe (2) states that 
volumetric packing factors of greater than 0.60 to 0.75 are 
Sceficult to achieve. A packing factor of 0.75, including 
product water tankage, has been selected for the prcomenon 
plant pressure vessel. Approximately 1/3 of the product water 
storage volume is provided within the propulsion plant pressure 
vessel sphere in otherwise unused space. The remainder of the 
product water is stored in a Separate product water sphere. 

Thus four spheres comprise the propulsion plant: a 
propulsion plant pressure vessel, two reactant pressure vessels 
{one for hydrogen and one for oxygen) containing cryogenic 
reactants, and a product water pressure vessel. Electric 
propulsion and associated hydraulic pumps and motors, startup/ 
emergency batteries, and a salt water circulating pump for 
propulsion plant pressure vessel cooling are mounted externally. 
meewarrangement of the spheres within the hull form of the sub- 
mersible should be such as to facilitate transfer of cryogens 
from their storage tanks to the propulsion plant pressure 
vessel sphere. 

The overall weights and volumes for the propulsion plant 
are given in Table XI-4as well as the required weight and 
volume of flotation material to achieve meReAL buoyancy. The 
assumed density for the flotation material is 34.5 isn ete 
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Subcritical Reactant Supercritical 
Storage System with Reactant Storage 
Liquid Propellant Rocket System with Porous 
Combustion Chamber Plug Combustion 


Lot 


Chamber 
Item 207 000 Vie S000 Tit 2 0r0 00) aes o, 000 set 

Weights 
Propulsion plant pressure vessel and 

internal components, weight, lbm 3997 2375 3997 2375 
Hydrogen storage pressure vessel weight 

ineiuding fuel, i1bm 7548 3548 6604 3654 
Oxygen storage pressure vessel weight 

mncluding exidant, lbom 4472 2 ia le 2 4998 3548 
Product water storage tank 1700 650 1550 600 
External components, lbm 1426 1426 1426 1426 
Total weight without flotation material, | 

Lbm 17 143 eal: S575 Tes 3 Oe 
Buoyancy 
Propulsion Plant pressure vessel, 1bf£ 2425 2265 2425 22605 
Hydrogen pressure vessel, Il1bf 5780 5500 4930 4590 
Oxygen pressure vessel, lbf 2340 VAY 2620 2390 
Product water pressure vessel, lbf Lg pas. 1590 a5 10 1a 
External components, lbf 264 264 264 264 
Total buoyancy without flotation material, 

ib bag 12,634 iMees so 11,749 10,924 

Net buoyancy, lbf Sep ie, +688 -6,826 S57), 
Required flotation material for 

neutral buoyancy @34.5 lbm/ft?, ft?/lbm 220.5/7610 - 231/7990 230/794 
Overall weight in air, lbm ZG D3 ly Sall 26,0200 1 Sem 

TABLE XI-4 


Overall Weight for 50 kw, 1000 kwh Propulsion Plant 





Syntactic foam is not currently available in this light density 
mer the 20,000 foot depth application. 
It is apparent that the two systems have comparable 


weights in air at both the 8000 foot and 20,000 foot depths. 
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CHAPTERS CEL 


Conclusions 


The hydrogen-oxygen fueled semi-closed Rankine cycle is a 
feasible propulsion plant for deep submersibles. With care in 
design of individual components, the plant can be designed for 
safe and reliable operation. The design goal of 4000-5000 hours 
time between overhauls probably cannot be met immediately, 
largely because of potential problems in the combustion chamber. 
Ball bearings for the turbine/generator unit may also shorten 
the time between overhauls somewhat from the design goal. 

It is difficult to assess the design lifetime of a combus- 
tion chamber. The experience in space flight with rocket 
combustion chambers has been for relatively short design life- 
times. High heat flux and a desire to keep weights to a minimum 
have resulted in combustion chamber designs which would not meet 
the goal of long-term reliability. The availability of relative- 
ly large amounts of coolant (diluent water) and, within reason, 
no Peesasity to minimize the weight of the combustion chamber 
should permit a more conservative and reliable combustion 
chamber design to he developed. The porous plug combustion 
chamber appears to have the greatest promise for essentially 
complete combustion and for reliability if problems of oxidation 
of the porous metal can be kept within limits and if the device 
behaves as predicted. This can only be proved through experi- 
mental work. The porous plug chamber would also allow the 


Beactant pumps to be eliminated. 
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The porous plug chamber permits the use of gaseous reac- 
tants at a saving of about 5% in reactant load over that for 
the plant with the liquid propellant rocket-type combustion 
chamber. The weight of the porous plug combustion chamber 
plane Could be reduced slightly if reactant pumps were used 
instead of the supercritical storage system. It is believed, 
however, that the supercritical storage system is more reliable. 
An anticipated better combustion efficiency for the porous plug 
chamber could also reduce the estimated reactant load, based 
upon a 98% combustion efficiency by weight. 

Overall thermal efficiency (or Beenately specific 
reactant consumption), reflected by overall plant weight in 
air, is not as good as was anticipated, largely because of 
parasitic loads. The largest of these loads is the condenser 
salt water circulating pump. The required pumping power of this 
pump is primarily determined by the size of the pressure hull 
penetration and internal piping. In the trade-off of pressure 
vessel weight due to buildup around penetrations with overall 
plant weight due to increased reactant load, it may be possible 
co improve the pumping power considerably. It is likely that 
other parasitic loads may also be reduced in magnitude by more 
careful component design and analysis. A conservative approach 
was taken in the determination of these loads. 

Overall plant weights in air for the proposed plant are 
heavier than the closed Brayton cycle plant proposed by 
Balukjian and Rackley (1) and outlined in detail in reference 2. 
Although the proposed plant includes certain equipment not 
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contained in the weight summary of references 1 and 2, notably 
startup/emergency batteries and the propulsion motors and 
hvcraulic system, one may conclude that the proposed plant is 
about 4300-7400 pounds heavier than the Brayton cycle for the 
20,000 foot, 50 kw, 1000 kwh mission. The variation in weight 
is apparently due to elimination of the need for a hard 
pressure vessel to contain product water in reference 2. A 
better comparison between the two plants is made by noting the 
required reactant load for the two sollamnts forgthe 20,000 £60t 
depth, 1000 kwh mission, 1711 pounds for the Brayton cycle 
plant and 2149 pounds for the proposed plant. At the 8000 
foot depth, the proposed cycle is about 1000 pounds heavier 
than the Benaeon cycle system. 

It is estimated that the proposed plant is about twice 
as heavy in air as a comparable cryogenic fueled, encapsulated 
hydrogen-oxygen fuel cell system. Variations in mission 
length and definition do not permit exact comparison. For a 
Pemeaeaole mission it is estimated that the proposed system is 
about half the weight of a silver-zinc battery system. 

Achievable turbine efficiency was found to be much better 
than anticipated for the 50 kw power level. This is a result 
of utilization of a combination of disks, each having several 
re-entry stages. This led to a relatively simple rotor design 
at the expense of a more complicated casing design with 
multiple return ducts. Good turbine off-design performance 
was found to be achievable if eOnocene speed operation is 


employed. The detailed analysis of leakage in the re-entry 
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turbine permitted more accurate assessment of expected turbine 
Pe ficiency than would have otherwise been possible. 

The detailed condenser analysis permitted reasonably 
accurate sizing of the condenser. The condenser analysis, 
although not a major portion of the analysis of the proposed 
plant, required extensive detailed work and represents a new 
contribution to condenser design. For this reason the computer 
program for condenser analysis has been included in this thesis. 

The proposed porous plug combustion chamber is the first 
known evaluation of the application of this device to high 
pressure combustion. It appears that the concept is feasible 
for hydrogen and oxygen premixed with steam and at pressures 
meta 600 lbf/in*. 

Although the proposed plant does not, in its present 
configuration, equal other proposed propulsion plants for the 
Same mission, it is considered that a more detailed design of 
parasitic loads could make it comparable with the Brayton cycle 
plant. The relatively high anticipated turbine efficiencies 
would permit this if parasitic loads can be sufficiently 
reduced. The development of reliable,very small reactant pumps 
in conjunction with development with the porous burner would 
assist in achieving this goal. The Brayton cycle plant has 
been subject to a much more detailed examination and development 


than has the proposed plant of this thesis. 
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APPENDIX A 


Liquid Propellant Rocket Combustion Chamber 


General 

The liquid propellant rocket combustion chamber is a 
proven device which, with the availability of large amounts of 
coolant in the form of diluent water, should provide a 
reliable high pressure combustion chamber for deep submergence 
applications. 

Its disadvantages are the high temperatures achieved in 
stoichiometric burning of such energetic reactants as hydrogen 
and oxygen. An additional difficulty is the problem of in- 
jection of the reactants at sufficient velocities to permit 
burning at some specified distance downstream of the injector. 
Given the fact that in most cases the reaction will be com-. 
plete within an inch or two of the injector head, a tradeoff 
between injection velocities and mean droplet diameter is 
required to prevent overheating and burnout of the head or 
Peeing of the injectors. Another difficulty, an outcome of 
small reactant flow rates and high injection velocities, is 
the extremely small injector hole size associated with low 
power applications. These problems may be overcome by the 
selection of the nroper type injector. If hole sizes become 
extremeiy small, laser "drilling" may be advantageous. While 
injection of both reactants in the gaseous state would be 
advantageous with respect to injector manufacture, it would 
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The prevention of uncontrolled combustion chamber 
oscillations is of major concern in any rocket combustion 
chamber design. While some conditions which would lead to 
uncontrolled oscillations may be possibly eliminated or delayed 
by the dececion of large amounts of diluent water spray down- 
stream of the combustion zone, this effect is not anticipated. 
Large pressure fluctuations could result from such a small 
exiting force as the passage of turbine blades downstream of 
the exit nozzle. For this reason, operation with choked exit 
flow from the chamber is desirable under all power conditions. 

The achievement of very high combustion efficiency (~99%) 
is essential to the overall efficiency of the plant and in 
Minimizing condenser size, Since any unburned reactants 
represent an increase in the amount of non-condensables in the 
condenser. Even a small fraction of non-condensables will 
Srgnstticantly decrease heat transfer. Such nom-condensables 
must also be stored at high pressure, pumped overboard, or 
feuened and condensed. Accurate control of the reactant feed 
rates is essential to the maintenance of stoichiometric 
conditions in the chamber. Any departure from stoichiometric 
conditions is also reflected as an increase in non-condensables 
in the condenser, as are impurities in the reactants. 

An additional problem posed by the selection of hydrogen 
and oxygen as reactants is their dissociation at temperatures 
above 3600°R. The dissociation reaction is temperature as 
well as pressure dependent. If, in the process of introduction 
of diluent water to the combustion gases to achieve the desired 
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turbine inlet temperature, cooldown is too rapid, the reassocia- 
tion reaction may not reach an equilibrium state within the 
combustion chamber. The reverse reaction (reassociation) 
proceeds fastest at higher temperatures and for that reason 
one would prefer introduction of a cooldown spray not in a single 
step, but rather in several stages, of which two would pro- 
bably be sufficient. 
High combustion zone temperatures, while detrimental to 
the injection head and walls if sufficient cooling is not 
employed, are more effective in achieving complete combustion. 
Quenching at the walls must be minimized to preclude large 
amounts of unreacted hydrogen and oxygen. It would appear, 
then, that the following items are desirable for high 
combustion efficiency: 
1) As high a combustion temperature as possible 
2) Minimizing of quenching at the walls by the 
use of as little wall cooling ds practicable 
. (minimizing cooled wall area) 
3) Provision for sufficient combustion chamber volume 


4) Adequate mixing in the combustion gases 


Previous Work 

Several gas generators burning hydrogen and oxygen have 
been reported in the literature. All of these have employed 
hydrogen rich mixtures to keep combustion gas temperatures 
within reasonable values. Ball (1) reported generally 


Satisfactory operation of a liquid oxygen, gaseous hydrogen 
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combustion chamber at chamber pressures from 47 to 232 lbf/in’ 
with a gas exit temperature ranging from 1260° to 2030°R. The 
total mass flow rate for this chamber varied between .18 to 
170 lbm/sec, roughly 8 times the total reactant mass flow rate 
of a 50 kw (shaft) hydrogen-oxygen semi-closed Rankine cycle 
propulsion plant at design conditions. This design utilized 
impinging hydrogen and oxygen injectors in the primary combus- 
tion zone. Dilution with excess fuel was accomplished by 
Spraying radially into the cylindrical chamber 12.5" down- 
stream from the injector head. 

Acker and Fenn (2) developed a liquid hydrogen, liquid 
oxygen combustion chamber which operated at a pressure of 
220 lbf/in?. Combustion chamber exit temperature was 1823°R. 
This chamber utilized a fuel rich mixture with an oxidant to 
fuel weight ratio of 1.0. A shower head injector was utilized 
in this design for oxygen injection while hydrogen was sprayed 
radially inward .05" downstream of the oxygen injector head 
from an annular plenum. Reported combustion efficiencies 
ranged from 93 to 96%, with the highest efficiency for a 
longer chamber. 

Sekas and Acker (3) reported the performance of an 
improved model of the device described in reference l, in 
which hot streaks of gas or oxidant had been responsible in 
PumEning Out thermocouples located in the exit nozzle of the 
chamber. Such hot gas streaks could result in disastrous tur- 
bine blading failures. The chamber developed by Sekas and 


Acker was to provide hot gas for a 1000 kw turbine. Several 
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types of injectors were compared in the development of this 
chamber. An efficiency of .992 was reported for a 13 inch 
(length), 2.0 inch (ID) cylindrical combustion chamber over 
the range of mass flow rate tested. This best efficiency was 
obtained with an annular injector design, i.e. hydrogen 
admitted through an annular nozzle surrounding each individual 
oxygen injector. This combustion chamber operated at a 
pressure of 240 lbf/in* with a hot gas temperature of 1823°R. 

Reference 4 reports the satisfactory testing of a high 
pressure liquid hydrogen-liquid oxygen fueled XLR-129 rocket 
engine developing 250,000 lbf thrust with a 2740 lbf/in’® 
chamber pressure. In the thrust chamber, fuel rich turbo 
pump exhaust gases mix with additional oxygen at a 6/1 oxidant 
fuel weight ratio, still fuel rich, but not far from stoichio- 
metric conditions. This engine is designed for 300 starts and 
a 10 hour lifetime, and is the forerunner of models for use in 
the space shuttle. 

It should be noted that all of the small gas generators, 
fPererences 1, 2 and 3, were designed for total operating life- 
times of minutes rather than the minimum 1000 hour lifetime 
deemed necessary for deep submergence applications. The 


smallest injector holes reported were .0135 inches in diameter. 


Design Requirements 
The most critical areas to be considered in the design 
of a liquid propeilant rocket type combustion chamber are as 


follows: 
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i) 


2) 


3) 


4) 


5) 


The achievement of very high combustion efficiency 
(99%) through 
a) Proper selection of injector type and 
arrangement 
b) Proper sizing of the combustion chamber for 
adequate combustion length and thorough 
mixing of reactants 
Smooth operation over a wide range of reactant 
mass flow rates through 
a) Injector throttling and stepping 
b) Control of reactant stoichiometry, dilution 
ratio and chamber exit temperature and pressure 
Adequate wall and injector head cooling to ensure 
safe, reliable operation over an extended lifetime 
(71000 hours) 
Reduction of combustion chamber pressure oscilla- 
tions and the elimination of uncontrolled oscillations 
through 
a) Proper injector selection 
b) Maintenance of choked flow in the exit nozzle 
c) Reactant and diluent feed pressures free of 
oscillations 
The ability to operate over a very wide range of 
combustion chamber pressures (~2:1 pressure ratio) 
which is desirable under constant turbine speed 


operation with varyina mass flow rates (This requires 
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reactant and feed water pumps whose speed can be 
varied continuously or in multiple svceps to meet 
the wide range of discharge pressures and flow 


rates required) 


Proposed Combustion Chamber Design 


The power levels of interest in this thesis range from 
50 kw to 500 kw. Since mission power profile of the 50 kw 
plant is known and because it is most requiring with respect 
to injector size, it was selected for a preliminary design 
analysis. The assumed conditions for the 50 kw plant are shown 
in Table A-1. For variable turbine speed operation, a constant 
combustion chamber pressure of 600 lbf/in* is assumed. 

Previous studies (5) showed that a specific reactant 
consumption (lb/kwh) of 1.68 could be expected for these 
conditions and a 1750°F chamber exit temperature. A 
regenerator effectiveness of 0.7 had been established, 
specifying the diluent water inlet temperature at 421°F. 

While these conditions do not represent the optimum cycle 
conditions for the 50 kw (shaft) power level, they are 


representative and allow use of tabulated data (7). 
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Oxygen injector sizes and stepping arrangements were 
studied considering the assumed power profile, and resulted in 
the selection of a swirl type injector. The wide range of 
throttling capability and ease of fabrication of this type of 
injector (larger size holes) led to its selection over the 
impinging and shower head types. It is estimated that a given 
injector may be reliably throttled to 50% of its design maximum 
flow rate. 

From Barrefte (6) one may estimate the Sauter mean diameter 
of the droplets for a given pressure drop. Any considerable 
amount of throttling will result in multiple stepping of injec- 
tors. This stepping must a See Laos rapidly and reliably 
to preclude any pressure or temperature excursions or extinguish- 
ing of the flame in the chamber. 

While the mission power profile assumes a cruising power 
of 40% of the design power leyel, it is realized that an ability 
to reduce power to 10% of its design level is desired when only 
‘auxiliary power is needed. For this reason injector stepping 
down to 10% of the design flow rate has been investigated. 

Since turbine efficiency would deteriorate considerably at 
these off-design power levels, the off-design flow rates would 
probably never reach less than 20% of the design maximum flow 
Bates. 

Since liquid hydrogen is a supercritical fluid at pressures 
above 200 lbf/in* it is much simpler to introduce it into the 
combustion chamber than licuid oxvaen. It can be easily 


admitted through a porous liner in the injection head, which, 
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in addition to distributing the fuel evenly, will provide vital 
cooling to the injector head. 

In order to determine the length of the combustion chamber, 
it is necessary to investigate the combustion process in some 
detail. Analysis of the combustion process in a rocket combus- 
tion chamber has been attempted by several investigators, 
emene them Spalding (8), (9), (10), (11), Williams (12), (13)) 
Priem (14), (15) and Salant and Toong (16). Each of these 
models had its advantages and disadvantages, but all represent 
certain oversimplifications to the combustion process, which is 
an extremely complex phenomenon. The assumption may be made 
that the controlling factor in the "burning" process of an 
oxygen droplet in a hydrogen atmosphere is the vaporization 
rate of the oxygen droplet. Salant and Toong (16) include 
counter diffusion of the fuel and combustion products as well. 
Using the simplified analysis of Spalding (10), it is possible, 
knowing a mean droplet diameter, to sete ne the distance from 
the injector head at which a droplet will be completely 
vaporized. This analysis includes droplet drag, using a Stokes 
Law drag term and a bulk gas velocity. Somewhat unrealistically, 
it assumes a constant injection velocity for all droplets as 
well as a single droplet diameter. Ignition delay is ignored 
as well, and the temperature exterior to the droplet is assumed 
to be everywhere the adiabatic flame temperature. The gas 
viscosity is that for the fully burned gas. Additionally, a 


full droplet spray is assumed to exist at the point of injection. 
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Despite all of the simplifications of Spalding's model, 
it does give an indication of the distance downstream of the 
injector head at which combustion is essentially complete. 
Using this as a basis, one may apply an "insurance factor" of 
2 to 4 for preliminary design purposes. Since the combustion 
chamber is small in comparison to the overall plant, it is 
reasonable to permit considerable "overdesign™ in size compared 
to rocket combustion chambers to insure reliability ame high 
combustion efficiency. For example, Spalding's analysis of 
the combustion region of the proposed chamber results in a 
combustion distance of about 2.0" for an injection velocity 
of 125 ft/sec and a Sauter mean diameter of 65u. 

One may use the method of Zucrow and Graham (17) to 
estimate convective heat flux to the walls of the combustion 
chamber if the Reynolds number (based on chamber diameter) is 
G@meater than 10° and if the mass flux of combustion gases is 
sufficiently high. This method probably may also be applied 
to transpiration cooling. Film cooling, however, would probably 
result in a more reliable distribution of liquid over the 
surface, and has therefore been selected for the proposed 
Giambper. 

Radiative heat flux was estimated by the use of charts in 
reference 18. Radiation heat transfer, at first thought to be 
relatively insignificant, was found to require a considerable 
amount of film cooling. While it was originally intended to 
produce re-association of the combustion gases by means of 
injection of a droplet spray into the chamber in the downstream 
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direction, it was found that film cooling alone in the combus- 
tion region was sufficient to reduce the bulk temperature of 
the combustion gases to 3600°R. This, of course, assumes 
complete mixing of the evaporated wall coolant and the bulk 
combustion gases, which will not be the case unless sufficient 
turbulence exists to produce it. For this reason turbulence 
promoters, to be shielded by adequate film cooling, are recomm- 
ended for placement downstream of the main combustion region, 
(about 4 inches from the injector head). Materials such as 
silicided molybdenum, which will withstand an oxidizing atmos- 
phere of 3600°R, might be useful for this application. 

If the remainder of the diluent is then introduced as a 
liquid spray, the required evaporation distance becomes signi- 
ficant, 9-10 inches for a 100u droplet spray injected at 
100 ft/sec. One also finds that so little diluent remains 
(roughly 15% of the total diluent) at the point where the bulk 
temperature reaches 3600°R that no injector pattern placed at 
‘the walls can be sized to produce a sufficiently uniform spray 
weress the chamber. This becomes especially evident when the 
chamber is throttled to 10% of its design maximum flow rate. 

For this reason the remainder of the diluent is split, 
roughly 30% of which will be used for wall film cooling down- 
stream of the combustion- re-association zone, and the remain- 
der introduced radially inward from an annular plenum as a 
steam spray. The exact proportions of this split as well as 
the distribution and sizing of the steam orifices in the 


annular steam plenum snould be the subject of experimental 
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work. Wall cooling should continue, either through film cooling 
or a regenerative cooling to a point where the bulk gas temper- 
ature is 2000°F. Turbulence promoters downstream of the point 
of discontinuence of wall cooling can be placed appropriately 

£O produce more thorough mixing prior to entry into the chamber 
exit nozzle. A centrifugal type separator may also be built 
into the chamber exit geometry to remove any remaining water 
droplets prior to entry into the turbine. A bypass valve to 

the condenser might be advisable for turbine protection during 
Seartcup. 

One finds that, for the smallest swirl injector, the mini- 
mum chamber diameter of the design is severly limited by the 
combustion distance of a droplet in the 50% throttled condition. 
Using the equations of Barrete (6) for this type of injector 
(with a .015" hole diameter at 50% flow rate) one may predict 
a pressure drop of 60 lbf/in*, corresponding to an injection 
velocity of 62.5 ft/sec. If one may use again the data of 
Kling shown in reference 6 to estimate the Sauter mean diameter 
of the droplet spray, a 210u droplet size is predicted with a 
corresponding combustion distance of 4.8". The pressure drop 
at maximum flow rate for the smallest injector would be 
220 ibf£/in* for a 65u droplet injected at 125 ft/sec and with 
combustion distance of 2.0". When accompanied by the very 
wide spray cone apex angle of 107° of this nozzle, a 
chamber diameter of at least 9" is required with injector 
placement as near to the center of the head as possible. As 
Ememinvector is throtti&td, the combustion zone will move 
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further and further from the injection head, from an average 
imo at maximum flow to 3.9" at»50% angectoneslo,. 

While the 9" chamber diameter is somewhat unrealistic 
with respect to uniform hydrogen distribution as the injector 
is throttled, both the 9" chamber and a 3" diameter chamber 
were utilized for detailed wall cooling calculations. A 
smaller chamber diameter would result in less radiative heat 
transfer per unit length and a slightly longer combustion-re- 
association length. Wall cooling calculations indicate, how- 
ever, that even for a 3-inch diameter chamber, the wall cooling 
Giluent is more than sufficient to reduce the bulk gas temper- 
ature to 3600°R within a distance of 4 times the calculated 
combustion length. In any event, a 6-8" length for combustion 
and re-association seems adequate and is not strongly dependent 
on the chamber diameter. 

The injector hole size chosen for the smallest injector, 
.015", was selected on the basis of ease of fabrication by 
conventional drilling methods. A smaller hole, which could be 
fabricated by laser techniques, might be possible for this 
application, but would be more susceptible to saemeanen The 
maximum injector flow rate for the smallest single injector is 
the governing parameter. 

It would then seem advisable to limit the threteling vot 
the combustion chamber to 20% of the design flow rate in order 
to achieve acceptable overall chamber performance. Under these 
conditions a .015" hole size injector would operate at maximum 


or nearly maximum Flow rate. Four other injectors of varying 


218 





hole diameters can therefore be sized at acceptable cone apex 
angles to.cover the remainder of the desired combustion chamber 
flow rates. The maximum throttling of any injector then is 
roughly 63.5% of the maximum injector flow rate. Figure A-1l 
shows the proposed injector arrangement, and Table A-2 gives 
details of the injector geometry and specifications. While it 
would be desirable to have the smallest injector operating 
continuously to prevent extinguishing of the flame during 
injector Peeping; this is not possible if both operating 
injectors experience the same pressure drop. When the larger 
injector is throttled, the smallest injector should remain at 
maximum flow rate, a situation which would require elaborate 
controls. An alternative, of course, 1s to use laser fabricated 
holes for the smallest injector. If the smaller holes prove to 
operate adequately without plugging, a continuously operating 
"pilot" injector would then be possible. 

Unlike a rocket combustion chamber, where axial gas 
velocity is high, carrying evaporating fuel or oxidant droplets 
rapidly downstream, the proposed rocket type combustion chamber 
exhibits very low axial velocities. Thus the oxygen is 
effectively sprayed into a stationary gas. For this reason, 
increase in the chamber diameter is not particularly effective 
in reducing chamber length. Figure A-2 shows a recommended 
design with a 4-inch diameter. The specifications for the 


chamber are given in Table A-3. 
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TABLE A-3 


50 KW Design Maximum Power Level (shaft) Combustion 
Chamber Characteristics (sized for 125% design maximum 


power level) 


Length of Combustion - Re-Association Zone, in....... 8.0 


Serculated Combustion Lengths, in.......- Sewers: odes eueraee Lider 


wordt Exit Steam Mass Flow Rate, lbm/hr........2.c+c2. 427 


ioeerie Reactant Mass Flow Rate, Lbm/hr.....cccssccssecs 106 
Ho oeeoeeee8« eee @ @ @ ee eoeeee#eeee e°e e eoee#s# oe @ eoeeeee e isis 
Oar errr cree cece reece reece crores oe @ 8 @ @ eoeeeese¢ ees 94.2 

te@eat Diluent Mass Flow Rate lbm/hr.......sssese0> « Sane oes 

unMmeMe eT JCCTCOLS « 6 secs s+ oo 6 web 5 oo 5 leo « seein e sarees 5 


Combustion - Re-Association Zone Film Coolant Evaporation 
PeeeINSIN/ 11 os oes ew os Ks Soe co © 0 slere wees + + eRe eee dG 
Cooldown Zone Film Coolant Mass Flow Rate, lbm/hr.... 48. 
meomdown Steam Mass Flow Rate, Ibm/hr......-cccesceces deg 6 
Meee tnternal Diameter, IN... scccsscccsseces cous 4 
Sitemoer Total Length, in....-scsceee ewer sorte. | LES 
Seer Exit Pressure, Lbf/in*?.......ceceesecneeces -. 600 


Chamber Exit Temperature, °F....... sie) seen enomemeteme suete oe oy 
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Ignition is best performed by a sitall gaseous hydrogen- 
gaseous oxygen igniter located adjacent to the main injector. 
head and initiated by a spark plug. Reference 1 reports 
Satisfactory operation of this device. If some film cooling 
is provided during chamber operation, burnout of the igniter 
should be preventable. 

As indicated in reference 1, startup and shutdown may be 
accomplished by admission of hydrogen followed by oxygen in 
the startup procedure with the reverse order during shutdown. 
Extremely fuel rich mixtures would be advisable during startup 
to preclude combustion chamber burnout when diluent spray has 
not yet eee introduced. 

In order that a constant combustion chamber pressure and 
exit temperature may be achieved over a range of mass flow rates, 
it will be best to select one of the reactants, either hydrogen 
Or oxygen as the controlling parameter to be varied with changes 
in chamber pressure. Selecting oxygen flow rate as the base 
variable, the relative amounts of hydrogen and oxygen may be 
sensed at the non-condensable gas exit from the condenser. 

The hydrogen flow rate may be adjusted to keep the hydrogen- 
oxygen ratio at the condenser exit at zero or at a stoichio- 
metric value, representing other than 100% combustion efficiency. 
A thermocouple in the combustion chamber exit nozzle may be used 
to detect stagnation chamber exit temperature and this value 
used to regulate diluent spray to the reassociation and cool- 
down nozzles by maintaining a specified base ratio between 


oxygen flow rate and diluent flow rate adjusted by chamber exit 
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Variable chamber pressure may be desired over a certain 
range of mass flow rates. For example, chamber pressure might 
increase from 600 lbf/in* at 40% design power level (shaft) and 
below to 1200 lbf/in* at 80% design power level and above. 
Under these conditions, sensing of steam mass flow rate from 
the chamber could be accomplished over this range and chamber 
pressure compared with a series of preset signals for various 
specified mass flow rates. Hence chamber pressure could be 
varied in discrete steps over this range. The other controls 
would remain essentially the same as before. 

Reactant pump and diluent water pump sizing and control 
are not considered in this appendix, since they are common to 


any high pressure hydrogen-oxygen combustion chamber. 


Assessment 

The rocket type combustion chamber appears to be a possible 
alternative for inclusion in the proposed propulsion plant. Its 
Major advantage is its capability for variable pressure opera- 
tion. Its main disadvantages are its extremely high tempera- 
tures and resultant heat fluxes to the chamber walls, the ever 
present possibility of catastrophic chamber pressure oscillation 
as in any rocket combustion chamber, and the difficulty of 
insuring against hot gas streaking in the exiting steam. An 
additional difficulty is the problem in achieving smooth 


operation during injector stepping. 
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APPENDIX 8 


Porous Plug Burner Combustion Chamber 


General 

A cooled porous plug burner has been used as a laboratory 
tool in the studies of laminar flames. The characteristics of 
laminar flames have been studied extensively through the use 
of a cooled porous plug burner. Kaskan (1) and Botha and 
Spalding (2) have shown that the velocity of propogation of a 
flame into an unburned gas mixture may be slowed if heat is 
extracted from the flame. Indeed this is the phenomenon which 
allows flame stabilization to occur. Prior to Kaskan's work, 
however, such stabilization was accomplished by extraction of 
heat at the edges of the flame, rather than uniformly through- 
out a flat flame front. This device has permitted the study 
of laminar flames in premixed, highly explosive gases even if 
they normally exhibit very high flame velocities. Figure B-1l 
is a schematic diagram of the porous plug burner. 

Siegler and Moore (3) have developed the porous plug 
burner from its laboratory use to a reliable engineering device 
for safely burning essentially undiluted stoichiometric mixtures 
of hydrogen and oxygen at atmospheric pressures. The particular 
industrial application for which the device was designed is the 
rTecombination of radiolytically decomposed hydrogen and oxygen 
removed from the condenser of a boiling water reactor plant. 
For this purpose, the device appears to be well suited and 
reliable. The burner, however, is not limited to burning 


hydrogen. To the contrary, undiluted stoichiometric hydrogen 
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gaa Oxygan, with their very high adiabatic flame velocity 
4000 cm/sec at 25°C and 1 atm), represent an extremeweace 
which the burner can safely and reliably meet. 

thie particular application in a boiling water reactor 
power plant requires that the device operate for long periods 
of time (essentially years) with little or no attention. 
These elements of safety and reliability, coupled with relative- 
ly low temperature combustion,make this device an interesting 
candidate for the semi~closed Rankine cycle proposed a this 


thesis. 


Principles of Operation 

The porous plug, fabricated from some porous material of 
high thermal conductivity, adequate porosity and sufficient 
strength, acts in three capacities: 

1) a flame stabilizer 

2) a flame arrestor 

3) a compact heat exchanger 

Its function as a flame arrestor is achieved by making the 
pore size sufficiently small to prevent flashback into the un- 
burned premixed gases. The quenching distance for an undiluted 
stoichiometric oxyhydrogen mixture at 25°C and atmospheric 
pressure is about 190u (4). By making the pore size a small 
fraction of this value, perhaps one-tenth, Siegler and Moore 
have satisfactorily prevented flashback from occurring. High 
values of open porosity of the plug material, while desirable 


from the standpoint of reducing pressure drop of the premixed 
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gases from the plug inlet to exit, recdwem the eifective 
thermal conductivity of the porous plug material, K.- This 
effective thermal conductivity has been investigated by 
Grootenhuis, Powell and Tye (5). 

Henee for a given heat flux to the porous plug, cooler: 
temperature, plug thickness and material, the temperature of 
the flame side of the plug will increase with increasing total 
porosity of the plug. Maintaining this temperature at a mini- 
MmmeLs important for two Begone. First, quenching distance 
decreases with increasing temperature making it possible for 
the flame to "march" into the porous plug. This would result 
in destruction of the porous nature of the plug, but probably 
not in an explosion of the premixed gases if the coolant flow 
through the plug is maintained. Secondly, if the porous nature 
of the plug is obtained by sintering, there is some temperature 
within 200-300°F of the sintering temperature at which the pores 
would begin to close up after periods of operation as short as 
several hundred hours (6). Thus a high sintering temperature 
for sintered porous materials would probably be desirable with 
as low a maximum operating temperature as possible, never 
exceeding a temperature 300 degrees F below the sintering 
temperature. 

Several porous materials of high conductivity may be used. 
Siegler and Moore (3) first used a commercial porous bronze 
plate but found that this had inadequate thermal conductivity. 
Finally several generations of sintered copper porous plugs 


were constructed from cut copper wire 180-250y in diameter. 
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In a 1700 hour "life test" of the final generation burner 
Gf reference 3, a 5,’ oxide thickness was noted on the copper 
particles upon metallographic examination. With a pore size of 
about 30yu (7), this did not adversely affect the operation of 
the burner. Only a slight increase in plug temperature and gas 
pressure drop was noted. While this did not prove to be a 
Major problem in the burner constructed by Siegler and Moore, 
it does restrict the use of sintered copper to applications 
where the downstream face of the porous plug can be kept at 
low values, somewhere in the range of 130-150°C. Porous 
Sintered nickel can be made with average pore sizes as small 
femosst with surprisingly good uniformity (6). The significantly 
higher thermal conductivity of copper and copper alloys would 
seemingly make sintered copper a more desirable choice than 
nickel, but this neglects the oxidizing atmosphere and the 
problems of reduction of pore size by oxide buildup. 

Minimizing the thickness of the plug on the flame side of 
the cooling passages is an obvious solution to keeping the 
maximum plug temperature within satisfactory limits, but this 
distance must be sufficiently great to achieve a uniform 
velocity distribution of the premixed gases at the exit (flame 
side) of the plug. Pressure drop across the plug and maximum 
plug temperature must therefore be traded off to achieve a 
Satisfactory design which will preclude flashback. 

The relationship of heat flux to the plug in conjunce en 
with flammability limits is another area in which a trade off 


is necessary. The heat flux absorbed by the plug, Q" er may 
Za2 





Bemcrelated to the mass flux of the premixed@qasces ac — 


Pe ou 
by the following relationship: 
pet = 99's Catan Ee ~ SP e%o\ tacit) | ai eotuamonr (B-1) 


meerecting pressure drop across the plug and noting then that 
u a’, for a stabilized flame front. The positive direction 


Wemin the direction of gas flow. 


G = mass flux, lbm/hr ft? 
u = flame velocity into the unburned mixture, ft/hr 
V, = the velocity of the unburned mixture upstream of 


ehe plug, £t/nr 

Bie — heat flux, Btu/hr fee 

) = density of the unburned gas mixture, including 
Griuent (ut any), lpom/re- 

AN = adiabatic flame temperature of the combustion 


product and diluent mixture, °R 


Te = actual flame temperature, °R 

C,, = mean .specific heat of burned gases and diluent 
between Te and Ee Btu/lbm°R 

i, = enthalpy of the adiabatic combustion product and 


diluent mixture, Btu/lbm 
i, = enthalpy of the combustion product and diluent 
mixture at the actual flame temperature, Btu/lbm 
The actual flame velocity into the unburned mixture is 
related to the adiabatic flame velocity, use by the following 


approximate relationship: 
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Og 7 (B/2RT,-E/2RT_) 


ie (BeZ) 
where 
= velocity of propogation of the flame fronteimieo 
the unburned mixture at the adiabatic flame 
temperature (no heat removal) 


E = an activation energy for the reaction, Btu/°R lb- 
mole or cal/°K g-mole | 

Torts in consistent units 

It is obvious then that, at zero flame velocity, there 
will be a condition of no heat flux. A similar condition will 
exist when Te equals To: At some intermediate temperature, 
Moore (8) shows that heat flux achieves a maximum value. 
Since flame velocity, u, is related to Use adiabatic flame 
velocity, by an exponential function, one finds that an 
increase in Te when Te<<T. results in a very rapid increase 
in u, the effect of which is to increase the heat flux even 
though the enthalpy rejected per unit mass of the burning 
mixture ee It is necessary to keep heat flux to the 
burner within acceptable limits. Reference 3 states that 
gomeal/cm-sec (2.7 x 10° Btu/hr ft*), roughly one third of 
the heat flux to the liner of a typical rocket combustion 
chamber, is a reasonable design figure. Operation of a porous 
plug burner at a flame temperature near the adiabatic flame 
temperature yet w.thin practical heat flux limits may be 


possible, but would occur with a significantly higher pressure 
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drop across the porous plug and with more critical material 
problems downstream of the plug due to the high combustion 
temperature. 

Moore (9) states that at flame temperatures of 900-1000°C 
and flame velocities of 2-3 cm/sec the hydrogen-air flame 
becomes unstable. This imposes a lower limit on flame 
temperature. 

Sregler and Moore (3) found that the plug tendedttouse 
self-regulating insofar as flow area is concerned. At lower 
than design flow rates of the stoichiometric hydrogen-oxygen 
mixtures, condensation formed in low flow regions, thus 
restricting the flow area and producing a nearly constant 
flame velocity. This effect was apparently achieved with a 
constant mass flow rate of coolant, which would produce such 
a phenomenon. Another self-regulating feature was the ability 
of the porous plug to maintain a uniform heat flux and thereby 
a stable flat flame. This can be explained by an increase in 
kinematic viscosity of the premixed gases in areas of high 
velocity (and consequently higher heat flux), which in turn 
reduces the velocity. 

The method of flame stabilization in a porous plug burner 
is briefly described below. The independent variable in the 
combustion/heat transfer problem is the mass flux of the 
unburned gas, ne The self-regulating feature of the flame 
is such as to exactly match the premixed gas velocity, v= 
with a flame velocity, u, into the unburned mixture, thus 


producing a stable flame in space at the exit of the porous 
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plug. This occurs because there is a heat sink available 
(the porous plug). 

In the high velocity regime, at velocities greater than 
the velocity for maximum heat flux, the stabilization process 
can be explained as follows: first assume that the flame is 
Stabilized at a given distance, y, from the exit of the plug. 
Next, consider a small perturbation to the distance, y, moving 
the flame front closer to the plug to some position y'<y. 

This in turn increases heat transfer to the plug, decreasing 
the flame temperature as shown by equation B-l. Since flame 
temperature and flame velocity are related through equation 
B-2, the flame velocity will decrease. This, in conjunction 
with a predetermined velocity of the unburned mixture, V ae 
will blow the flame front downstream to some point y">y, where 
the heat transfer to the plug will be sharply reduced, again 
increasing the flame temperature and flame velocity. 

In the low velocity regime, that of operation of the 
porous plug burner, the stabilizing process is somewhat 
different as explained by Wohl (10). There exists a minimum 
distance from the burner at which the flame can exist. Ata 
lesser distance quenching will occur. As the velocity of the 
unburned mixture, Vor is decreased below its value corresponding 
to this distance, the flame then marches away from the burner 
to a new equilibrium point at a distance greater than the 
minimum value. Thus a flame can actually exist at a given 
distance above the minimum distance from the burner at two 


different flame velocities. The minimum distance actually 
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Corresponds to a velocity slightly greater than that for 
maximum heat flux (7). 

Thus the laminar flame speed, u, is set by the approach 
velocity, Dat Since the flame velocity, u, will equal es as 
long as there is sufficient cooling available to remove the 
tweae from the flame front. This will occur theoretically at 
any coolant exit temperature below the flame temperature 
corresponding to u. Practically speaking, however, the coolant 
temperature, the porous plug geometry, particularly the thick- 
ness and its effective thermal conductivity (based on plug 
Material and porosity), and the heat transfer coefficient to 
the coolant in conjunction with the heat flux absorbed by the 
plug, a" es establish some fixed plug exit temperature. Clearly 
these are design parameters which must be tailored to achieve 


an acceptable plug exit temperature. 


Application to the Semi-Closed Rankine Cycle Propulsion Plane 


Premixing of the reactants in a gaseous state would be 
advantageous for a combustion chamber for the hydrogen-oxygen 
fueled semi~closed Rankine cycle propulsion plant. Of course, 
it would be necessary to demonstrate that such premixing could 
be accomplished without any significant danger of explosion. 
Premixing of the fuel and oxidant in a gaseous state rather 
than injection as liquids will improve the combustion efficiency 
of a combustion device. This will occur because of improved 
homogeneity of the fuel-oxidant mixture in the chamber. An 


additional reason for premixing in the gaseous state is the 
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ease with which such burner can be throttled to lowes 
meactant feed rates. The injection pattern is not so critical 
mowcombustor performance here as it is in a liquid propellant 
rocket type combustor. Therefore, reactant feed rates may be 
throttled with reasonable ease without the danger of producing 
combustor instabilities and associated excessive heat fluxes 
to combustor walls. 

Another advantage of the porous plug burner is its lower 
combustion temperatures, which preclude dissociation and 
require less or perhaps no wall cooling. In general, -it is 
advantageous to lower the combustion temperatures as far as 
possible, since this reduces material problems and improves 
reliability. 

The application of the porous plug burner to stoichio-~ 
metric oxyhydrogen mixtures at high pressures requires 
investigation of several factors. First, very little investi- 
gative work in the field of oxyhydrogen combustion at pressures 
in excess of one atmosphere has been done. However, it can be 
readily seen that such factors as quenching distance, adiabatic 
laminar flame velocity and ignition delay combine to produce 
emmere stringent design than at atmospheric pressure. One 
way of reducing the deleterious effect of pressure on these 
Peameters is to add inerts to the premixed gases. This, 
however, is only possible for the proposed application if 
steam is used as the diluent. Saturated steam at high pressures 
has associated hiyh temperatures (400-500°F), and, if added to 


premixed gaseous oxygen and hydrogen, will preheat the mixture. 
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This effect will negate to some degree the effect of adding an 
inert. 

The effect of high pressures on adiabatic laminar flame 
velocity was investigated by Edse (11), who measured these 
velocities using a bunsen burner-like apparatus at pressures 
up to 14.6 nee oe His is the only work found which 
reported results at high pressures. Edse, however, was forced 
to make estimates of laminar flame speed from turbulent flames, 
Since at higher pressures, the flame speeds were such ae Lo) 
produce turbulent flames. In order to estimate laminar flame 
velocity, Edse measured the angle which the flame front made 
with the rim of the burner tube at the supposed intersection 
of the laminar sublayer with the turbulent core. Such 
estimates cannot be considered to be highly accurate, but 
probably do reveal trends. Edse found that, in an undiluted 
stoichiometric oxyhydrogen mixture, laminar flame velocity 
increased from about 1000 cm/sec to about 3400 cm/sec when 
pressure was increased from 1 to 14.6 atmospheres. With 
undiluted richer or leaner mixtures, the effect of pressure 
became less significant. The effect of pressure on flame 
velocity on very rich and very lean mixtures was not investi- 
gated but appeared to be negligible for rich mixtures (>85% H.) 
and appeared to increase by a factor of 2 from 1 to 14.6 
atmospheres for lean mixtures (<30% Ho). 

More accurate work was performed by Manton and Milliken 
(12), who used the spherical bomb method for determining 


laminar flame velocity. Their work, however, was conducted at 
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pressures of less than 1.5 atmosoheres. It showed, however, 
that adiabatic laminar flame re of all the gases 
tested, which included propane, ethelene, acyetylene, methane 
and hydrogen with various amounts of diluent argon, helium or 
nitrogen could be correlated by use of a pressure exponent, x, 


by the following equation: 


meat (ae) * | : 
) (B-3) 


The unprimed quantities indicate the values at some reference 
pressure, p, and the primed quantities at some other pressure. 

Manton and Milliken's work showed that at adiabatic flame 
velocities between about 50 and 100 cm/sec, there was no 
pressure effect whatsoever. At higher flame velocities, the 
pressure effect became more significant, predicting a laminar 
flame speed of about 2440 cm/sec for undiluted stoichiometric 
hydrogen and oxygen at 14.6 atmospheres, as opposed to the 
3400 cm/sec velocity estimated by Edse at this pressure. 

Flame theory (13) predicts that laminar flame velocity 
is related to pressure and temperature by the following 
expression: 

ee sive oe E/2RT (B-4) 
where n = the order of the reaction (the number of molecules 


which react in a chemical reaction) 
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The hydrogen-oxygen reaction, aithouy.: somewhat complex, 
is probably the best understood of all multistep reactions. 
This sinertle not imply, however, that the reaction is sufficiently 
well understood that one may always make reliable predictions 
as to its behavior. Basically, the hydrogen-oxygen reaction is 
a two body process, which should make it independent of 
pressure. 

This apparently is not the case. Edse explains this by 
the effect of pressure see dissociation phenomenon at high 
temperatures. It is known that dissociation is not significant 
for this reaction at temperatures below 1800°K. When the gas 
temperature is Bee to this value or less by means of 
Gilution, at one atmosphere adiabatic flame speeds on the order 
of 300 cm/sec result, with a less significant pressure effect. 

A further effect of lowering flame speed is produced by 
dilution itself and may be approximated by 

+ 

aa ' Od 

a” tir 16,1 “ 
2 2 

where the brackets, [], represent molal or molecular concentra- 

tions of the reactants atagiven pressure and temperature and 

the primed quantities indicate values with dilution. 

Equation B-5 is a consequence of the dependence of flame 
velocity on the square root of reaction rate. The reaction 
Mee Constant is in turn proportional to the product of the 


concentrations of the reactants (the collision frequency term). 
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Again, since the hydrogen-oxygen reaction is not a simnle 


reaction involving only H and 0., but also many intermediate 


2 
species, this expression is only approximate. 

A comparison of equation B-5 and observed laminar flame 
velocities at standard conditions for undiluted stoichiometric 
Bearogen and oxygen and for stoichiometric hydrogen in air 
under the same conditions would predict a laminar flame speed 
for air of 147 cm/sec, correcting for the change in adiabatic 
flame temperature and dilution. The observed laminar flame 
velocity for stoichiometric hydrogen in air ranges from 193 to 
232 cm/sec with an average of 215 cm/sec (14). Thus equation 
B-5 predicts a laminar flame velocity of about 70% of the 
observed value. 

If equation B-5 is used in conjunction with a constant 
correction factor of 1.46 (based upon the ratio of the observed 
versus the predicted laminar flame speed for air) one can 
probably predict the effect of dilution with steam with some 


reliability. Then 


a {#1 [0,11}4 
— = 1.46 = —“—__+_,— (B-6) 
a Eon i 


is an approximation of the dilution effect by decreased 
collision probability on adiabatic laminar flame velocity. 
The other effect of dilution, of course, occurs through its 


effect on adiabat' c flame temperature, as seen in equation B-4. 
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So there are actually four corrections to make to the 
adiabatic flame velocity of undiluted hydrogen and oxygen at 
standard conditions to arrive at a predition of adiabatic 
flame velocity at elevated pressures with steam dilution. 

a) the dilution effect due to decreased collision probability 

b) the effect of a change in adiabatic flame temperature 

due to dilution 

c) the effect of preheating 

dad) the effect of increased pressure 

The effect of pressure on quenching distance is primarily: 
a function of its effect on the term, Pua: Since according sto 


Williams (13) 


eT 
wi 





g (B=7) 


pPova 


Qs 
ie 
Qt 


Where d = quenching distance, ft (oru) 


a = a Peclet number based upon quenching distance and 
adiabatic flame velocity, about 40 
Ky = average thermal conductivity of the premixed gases 


between the initial temperature and adiabatic flame 
temperature, Btu/hr ft°R 
c., = average specific heat of the premixed gases between 
the initial temperature and the adiabatic flame 
temperature, Btu/lbm°R 
Mm. is found that Ky and C,, are nearly independent of pressure. 


Thus, as a first approximation, one may expect quenching 


distance to vary inversely with pressure. This effect becomes 
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more pronounced if in turn the pressure exponent, x, in equation 
B-3 is positive. Preheating the mixture, of course, may increase 
Kg and cy slightly and will decrease Pee but at the same time, 

us will increase. For reasonably small increases in temperature 
(preheat temperature <500°F) significant changes in quenching 
distance are not expected. 

No data was found on the effect of large amounts of steam 
diluent in oxyhydrogen flames, but it is expected that steam 
will act in a way similar to nitrogen or carbon dioxide. 
Reference 14 states that, insofar as flammability limits are 
concerned, water vapor behaves approximately like carbon 
dioxide. Hence one may use charts of flammability limits of 
References 14 and 15,and, for nearly all mixtures of interest, 
it appears that there is no problem of flammability. Reference 
14 indicates that the pressure effect on flammability limits 
should not be significant. 

Mixing superheated steam at eoipevatiinee in excess of 600°F 
with stoichiometric hydrogen and oxygen at elevated pressures 
might be cause for concern with respect to spontaneous ignition. 
Although the rate equations are somewhat complicated for the 
hydrogen-oxygen reaction, if one ignores reactions at the wall, 


reference 14 predicts the following equation for ignition lag 


l eS 
T= C 25 lige oe (1-— My? (B-8) 
e p 6 
where 
T = ignition lag, milliseconds 
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p = pressure in aporopriate units 


kor Ke = reaction rate constants in the notation of 
Lewis and Von Elbe (16) 

c = constant 

[M] = concentration of all components of the mixture 


(molecules/cc) 


The associated chemical reactions are as follows: ., 


: 
H +0, 2 OH +0 


ne 

fm + 0. Meee HO. + MM 

Ignition lag has been found to depend upon the mole 
fraction of hydrogen in a mixture and not on the relative 
amounts of oxygen and inerts. Charts of ignition lag versus 
mole percent hydrogen at atmospheric pressure and a tempera- 
ture of 649°C (Anagnostou's data shown in Reference 14), were 
used in conjunction with Lewis and Von Elbe's method (16) for 
calculating Ke and ky to determine the constant, c, in equation 
B-8 for an air-hydrogen mixture. Subsequently, equation B-8 
was used to calculate the ignition lag at a lower temperature 
Of interest (440°F) and at 600 lbf/in* total pressure with a 
mixture of the following composition (mole fractions): 
hydrogen - .25, oxygen - .125, steam - .625. An ignition lag 
in excess of 50 hours was predicted. Reference 14 substantiates 


this finding. The data of Zabetakis on H, - water vapor-air 


2 


mixture explosion hazards reports spontaneous ignition temper- 


atures of 580°C for 30% water vapor at a pressure of 7.8 
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atmospheres. Hence it appears that spontaneous ignition 
should not be a problem. 

The method of mixing the steam with the hydrogen and 
oxygen will be important and should be the subject of experi- 
mental work. Because of reported ignition of hydrogen in air 
by static electricity discharges, the steam should first be 
added at low velocities to only one of the two reactants, 
followed by the other reactant after sufficiently uniform 


mixing has occurred. . 


Design of a Selected Combustion Chamber 

It is noteworthy that several basic differences exist in 
the operating conditions of the porous plug burners of Siegler 
and Moore and the high pressure combustion chamber investigated 
in this thesis. First, Siegler and Moore's burners operated 
with effectively unlimited coolant supply. In the proposed 
chamber, the coolant available is limited by the ratio of the 
mass flow rates of diluent (feed) water and the combustion 
products. This ratio is established by the desired turbine 
inlet temperature. While more diluent water at a cooler 
temperature could, of course, be provided, it would represent 
wasted energy with a consequent reduction of overall thermal 
efficiency of the power plant. Furthermore, under reduced 
power conditions, a plant with a regenerator will exhibit in- 
creased coolant inlet temperatures because of decreasing 
turbine efficienc., a condition common to gas turbine power 


plants. With a coolant inlet condition already near that of 
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the saturated liquid state at design conditions, it is probable 
that nucleate boiling will occur very near the entrance of the 
porous plug cooling coil under reduced power conditions. This 
situation virtually eliminates the possibility of any self- 
throttling feature as experienced in the burners of Siegler and 
Moore. Of course, with higher power level plants with high 
expander efficiencies, the coolant inlet temperature will be 
much lower since a regenerator will not be necessary. 

Secondly, if the proposed combustion chamber requires 
dilution with steam, the effluent of the porous plug cooling 
coil must be superheated prior to mixing with the unburned 
hydrogen and oxygen, assumed to be at 75°F. This superheating 
would be necessary in order to produce steam which is at least 
saturated at its partial pressure in the premixed gases. 
Superheating the steam creates material problems with a super- 
heating coil in a high temperature environment downstream of 
the porous plug combustion zone. Even the best high tempera- 
ture nickel-based superalloys are operating at their limits at 
temperatures much in excess of 2000°F. 

Another problem in the proposed combustion chamber is the 
tolerable heat flux established by Siegler and Moore (20 cal/ 
em? sec) and its corresponding flame velocity. For the 
effectively undiluted stoichiometric hydrogen-oxygen burners, 
Operating at near atmospheric pressures, the flame velocity 
corresponding to this heat flux is about 20 cm/sec. As shown 
in subsequent parugraphs, the flame velocity in the high 


pressure burner is about 3-4 cm/sec. Thus the proposed 
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chamber operates very near the unstable flame speed limit 
previously mentioned (2-3 cm/sec). All of these factors, - 
coupled with temperature limits on the porous plug material, 
combine to produce a more stringent design situation than that 
existing at normal pressures. 

To determine the feasibility of a combustion chamber, 
certain assumptions must first be made. The chamber pressure 
selected was 600 lbf/in?, and an exit temperature of 1750°F 
was chosen on the basis of previous studies which had angen 
these values to be suitable for a high overall thermal 
efficiency. 

A 50 kw design shaft power was selected with a capability, 
insofar as the combustion chamber is concerned, to overate 
indefinitely at power levels 25% in excess of design power. 
Table B-1l gives a summary of plant parameters for the combus- 
tion chamber feasibility design. 

Several methods of burning the hydrogen and oxygen are 
available, from burning with no diluent steam to burning with 
all of the diluent water available in the premixed gases. 
Calculations showed that burning with no inert at 600 lbf/in? 
and a flame temperature of 1300°K (the probable lowest flame 
temperature at which the flame will be stable) results in 
excessive heat flux to the porous plug, on the order of 6 x 10° 
Btu/hr ft? (460 cal/cm? sec). 

For this reason, the effect of various amounts of diluent 
steam in the unburned mixture were investigated and for each 


mixture the heat flux to the porous plug was calculated. 
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A RGTSN Dee S Tesi 
Summary of Design Conditions 
for Porous Plug Combustion 


Chamber Feasibility Design Study 


Meember exit pressure, lbf/in* ........ +... «. 600 
Maemoer exit temperature, °F ... .. + «6 « « « » « »L/50 
Mumm rated shaft power, kw. . . «© » « « « «© »« « «© « « 990 
Membustion chamber maximum power rating, kw. ..... 62.5 
Memeeene Cfficiency . 2 « < 2. « « «© 6 s « » « #6 «= « « « 66% 
Mupeerle plant G@LEVCLENCyY 2 si) 666 ee es Vere 315% 
Recirculation ratio, (m diluent/m reactants) ..... 221.8 
Total chamber exit mass flow rate at 125% overload, 

OVA) a ee ee rere eras Soren 22 7) (0) 
Reactant feed rate at 125% overload lb/hr. ..... . 102.0 

Ho: Os eae. a ee See, ie Se See ee ee 

0.- ae ae, ae eee ye ee ee Ae ee mee!) 
Feed water pressure (diluent) Mp flea. ee 2 ere COle 
Peed water inlet temperature, °F . ...... «+... 420.6 
menous plug cooling coil inside diameter, in..... 02236 
Porous plug cooling coil outside acne. I Ieestts ey )e) ot ee 0.250 


Member inside diameter, in. .....«». « « « « «ew. 11.7 
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Because of the desirability of maintaining uniform heat flux 
and consequently uniform temperature at the exit plane of the 
porous plug, a cooling coil geometry in the form of a spiral 
of Archimedes was selected. The diluent water (feed water) 
would enter at the periphery of the plug (near the combustion 
chamber wall) and exit at the center as saturated steam (see 
Figure B-2). The superheating, it was concluded, would probably 
be best accomplished by use of a separate superheating coil 
downstream of the porous plug burner in the path of the hot 
combustion gases or imbedded in the combustion chamber wall. 
Superheat is necessary to produce a saturated steam condition 
in the premixed gases when the steam is added to 75°F hydrogen 
and oxygen. 

Except for a short section of natural convection heating, 
the majority of the porous plug coil would exhibit local and 
bulk boiling. Figure B-3 is a graph of the lengthwise tempera- 
ture distribution of a coil for the design conditions selected. 
Using the Jens and Lottes correlation (17) for nucleate flow 
boiling, it is seen that a uniform wall temperature exists 
Eaeoughout most of the coil. | 

At some flow less than the maximum flow rate, it is 
evident that the instability limit of flame velocity will be 
reached. Depending upon the mixture selected, the temperature 
of this occurrence may vary slightly. Calculations show, 
however, that a reduction in flame velocity (mass flow rate) 
of about 25% is possible before the flame stability limit is 


approached. This information is based upon an initial flame 
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temperature of 2000°F and a corresponding flame velocity of 
about 4.0 cm/sec. 

In the calculation of temperature distributions within 
the plug, it was assumed, for the sake of simplicity, that 
uniform heat flux to the round coolant tube exists, i.e. wall 
thermal conductivity is sufficiently high to permit uniform 
heating of the coolant walls (on the upstream and downstream 
Sides with respect to gas flow). Calculations for the design 
condition selected poveaiea that a uniform round tube heat 
flux is only slightly less than round tube critical heat flux 
predicted by MacBeth's correlation (17). Thus if the design 
is not carefully accomplished, partial film boiling may result 
with a consequent increase in tube wall temperature and a 
disruption of the uniform temperature distribution in the porous 
plug as well as an increase in plug exit temperature. Further 
calculations show that proper sizing of the cooling coil and 
proper spacing (sizing of the gap) of the spiral can adjust the 
Operating conditions sufficiently to maintain an acceptable 
Meegin to critical heat flux. 

Table B-2 shows calculated values of heat flux for several 
mixtures of stoichiometric hydrogen and oxygen with steam as the 
diluent. Laminar flame velocities have been corrected for 
flame temperature (including preheating) dilution and pressure. 
It is seen that the mixture temperature (including preheating) 
dilution and pressure. It is seen that the mixture temperature 
does not vary much from 440°F. The temperatures shown for the 
premixed gases correspond to saturation temperatures for the 
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Mixture 
No. 


1* 


NOTES: 


Mole Fraction 
Steam Premixed 
5 estes 
5075 
765-0 
nOZo 
-600 
275 
355 0 
no2 5 


2200 


fxx 


Bees: 


-876 


Se(iisits 


5 7 le 


~641 


oS 


e525 


~471 


Mixture 
Temperature 


“8 
446.9 
445.9 
442.2 
439.9 
434.5 
430.4 
426.3 
421.9 


417.4 


13 
on/°K 
2480/1830 
2870/1850 
3035/1940 
3210/2040 
3390/2120 
3552/2225 
3740/2334 
3872/2405 


4060/2510 


TABLE 52 


mee 


Of /°K 


1840/1280 
E9027 1310 
20087 LoGe 
2117/1430 
2250/1500 
235977 1560 
2412/1635 
2580/1685 


2722/1768 


u 


a 
cm/sec 


22,0 
25,8 
51.4 
82,3 
102 
172 
239 
312 


439 


*This condition denotes all diluent water added to the unburned mixture 


**Mass fraction of total diluent flow which is used for diluting the usmburned gases 


Assumed E = 48,000 cal/g mole 


Total Pressure ~ 600 lbf/in? 
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u 


cm/sec 


10. 
liye 
25. 
36. 


Sere 


210 


Bi 


> 9, 


MEBs, 


Heat Flux to Burner Versus Unburned Mirture Composition 
G(T Pe) 
Btu/lb 
610. 

a Ailbe 
674. 
726. 
781. 
840. 
901. 
S155 


1020. 


A 


ii 


Tf 


q 2 
Btu/me EG 


cal/cm 


2 


= 


.66x10°/5.76 
44x10°/10. 80 
.41x10°/25.6 
260/472 
.84x10°/74.0 
.81x10°/136.0 
89x 10- 7 2agee 
28x10 732400 


.28x10°/547.0 


/ 


OV 
OO 
lb/hr ft? 


126. 


235. 
507. 
865. 
1DaD- 
2160. 
3210. 
4490. 
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partial pressure of steam in the mixtures. Comparison of the 


laminar flame velocities predicted for these mixtures with 


those observed by Kaskan (1) for air and CO. "air" (nitrogen 


2 
replaced by an equal fraction of co.) are favorable. A some- 
what unexpected effect is the increase in heat rejected from 
the flame (the Ss (To = Te) term in equation B-1) with increasing 
Tee This occurs because adiabatic flame temperature increases 
with decreasing amounts of diluent steam in the unburned mix- 
ture. Assuming that no flame temperature below 1300°K will 
maintain a stable flame and assuming that no heat fluxes in 
excess of 5 x 10° Btu/hr ft? (37.6 cal/cm* sec) can be absorbed, 
mixture 3 was selected for the feasibility design. 

Calculation of the porous plug gas side pressure drop and 
axial temperature distribution requires an accurate assessment 
of quenching distance. Friedman's work (4) gives a method of 
predicting quenching distance but it is based on an ignition 
temperature, itself a somewhat fictitious quantity. His data 
does, however, show the effect of the addition of nitrogen to 
a stoichiometric hydrogen-oxygen mixture at one atmosphere. 
If an equivalent amount of nitrogen is calculated for hydrogen- 
oxygen-steam mixture 3 of Table B-2 (allowing for the preheating), 
one may predict a quenching distance of Oe atmosphere 
Beane mOle fraction of nitrogen of 0.71. 

Ignoring a probable small temperature effect on quenching 
distance and noting that the ratio if Kg/Cp for air and steam 


are not significantly different, one may conclude that the 


pressure effect on quenching distance alone is Significant. 
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Hence a quenching distance of 37u is a reasonable estimate for 
mixture 3. A 1.8 pore size was selected as a conservative 
ieare (possibly overconservative) . 

The temperature distribution within the plug is a function 
of heat flux, mass flux, porosity, plug thickness, material and 
coil placement. If a one-dimensional condition is assumed, the 


heat equation may be written in the axial direction, 


di s= idol. 
So Gis ay Ghee ee a) 


with the boundary conditions at the plane of the cooling coil 
placement. 

If the plug is separated into two regions upstream of the 
plane of the cooling coils (I), and downstream of the cooling 
coils (II), the boundary conditions for region I, (see Figure 
B-4) become at the right hand boundary (x = &5) 


at = T 


11 a 


aa 


mies); = —G C. (T(x ,)-9,) = —ko ae (8) 


ee 


here it is assumed that the temperature of the gas and the 
porous plug temperature at any axial position, x, are equivalent 
and that axial conduction in the gas is negligible. 

For region II, the boundary conditions are also known at 


the plane of the coolant tubes, where 
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pe ited Ve ae 
where 
T, = the outside temperature of the coolant tube, assumed 
to be of zero thickness in the chamber axial direction 
q" = the frontal heat flux absorbed by the porous plug, 


equation B-1l 


c., = average gas mixture specific heat between inlet and 
exit for a given region 
kK. = effective thermal conductivity of the porous plug 


The actual heat flux at the exit plane of the burner 
(x = Xo) is that required to heat the premixed gases to the 


plug exit temperature plus the frontal heat flux, Q" er ie 
eo ry = "2 Conegry “OUD Teg) toe 


recalling that a" ¢ is negative in sign. 

The following conditions were assumed for the porous plug 
Materials (data is that of Laverty (6)): 

Material - carbonyl nickel power - INCO A(Typical analysis 


-1$C, .1% 0 ,00L SS, .0L Fe) 


Oe 
miGtal porosity, %s°- 18.0 
Open porosity, % - 15.9 


Average pore size, - 1.81u 
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Pore size range, -0.5}: 

Sintering temperature, -F —- 19007F (H, atmosphere) 

Assumed thermal conductivity (fully dense), Btu/hr ft°R - 

36.0 
Assumed effective thermal conductivity, Ker Btu htwsine R — 
ee 0 

Using this data, a plug exit temperature of 1290°F was 
calculated for a plug of thickness 3/8" downstream of the plane 
of the coolant tubes with a total thickness of 1/2". 

The marked reduction of thermal conductivity of sintered 
Mtiectials predicted by reference 5 is corroborated by more 
detailed and recent work of Tye (18) and Koh and Fortini (19) 
for sintered stainless steel powders, "rigamesh" porous stain- 
less steel, and sintered copper powders. This is also true for 
Sintered porous nickel (20). For a total porosity of 20% a 
thermal conductivity of about 50% of the solid material is 
expected. Furthermore, the assumed thermal conductivity for 
pure nickel of the fully dense material is perhaps a little 
high for nickel powder. A value of 33.5 Btu/hr ft°R (20) is 
probably more representative of nickel powder in the tempera- 
Gare range of 500-1000°F. 

Several pressure drop calculations were performed on the 
1/2" plug, resulting in pressure drop predictions ranging from 
30-1400 lbf/in*. One pressure drop correlation used was that 
Of Grootenhuis (21), which compared well with a correlation of 
Rose for packed beds of granular materials. Calculations 
using another correlation for packed beds (22), that of Ergun, 
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m2eided a pressure drop about the same as by the method of 
Grootenhuis, about 1400 lbf/in?’. ROCHoEN more recent method 
(23) predicted a pressure drop of about 30 lbf/in?, and this 
figure was verified as being reasonably accurate and represen- 
tative of pressure drop in filtration at high pressures with 
filters of very small pore size for the design flow conditions 
(24). Certainly experimental work will be required to verify 
the pressure drop. 

In the tradeoff of the effective thermal conductivity, 
pressure drop and pore size sufficiently small to prevent 
flashback, it may be possible to increase pore size to reduce 
pressure drop if necessary. Another possibility is to reduce 
the size of the coolant tubes, already desirable from the 
aspect of critical heat flux, and thereby reduce the thickness 
of the porous plug. This would result in a reduced plug exit 
temperature as well as a decreased pressure drop, dependent on 
plug thickness, of course, ym also on exit temperature through 
gas viscosity. This would also benefit in reducing oxidation 
of the porous material. Siegler and Moore (3) used 1/8" OD, 
0.020" wall thickness copper tubing for their application. 
Quenching distance studies must be performed with the design 
mixture of steam, hydrogen and oxygen to determine the maximum 
allowable pore size. Pressure drop studies of the sintered 
nickel and cooling tube compact must be performed to determine 
the actual pressure drop. The strength of the compact must 
also be determined in the proposed application. 

In spite of the questionable character of the pressure 
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Meop crecdictions, a preliminar 
combustion chamber was performed. If the pressure drop can be 
fmept at 100 lbf/in* or less, the burner could probably operate 
successfully. 

The superheating arrangement of this chamber, located 
downstream of the porous plug, can be a very simple device, 
Since the temperature difference between the bulk gas tempera- 
ture and the maximum superheated steam temperature is over 
1000°F. Table B-3 gives a complete description of uae super- 
heating coil arrangement. Several high temperature nickel- 
based alloys are suitable for the superheater tube material. 

An alternate superheating arrangement is placement of the 
superheating coil in the combustion chamber wall. Both 
arrangements are shown in Figure B-5. 

After superheating, a certain fraction of the steam is 
mixed with the unburned hydrogen and oxygen upstream of the 
porous plug. The remainder of the steam is mixed with the 
combustion steam downstream of the superheating coil, reducing 
the combustion gas temperature to a final temperature of 1750°F. 

Figure B-5 is a schematic diagram of the overall combus- 
tion chamber. Detailed arrangements to improve mixing in the 
premixed chamber and downstream burner will require experimenta- 
tion. Throttling must be provided by sectioning the porous 
plug with one pilot section, probably centrally located, which 
will operate continuously. Each section of the porous plug is 
separated from the remaining sections by solid walls, and each 


section would have its own premixing chamber. 
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PABLE B= 
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Superheating Coil 


Peemangement. .« . « « « « « » « « « 


Tube 
Tube 
Tube 


Tube 


Gutside diameter, ins . « . . 
maside diameter, in... . - 
spacing, center to center, in 


agit, ) ciwere tw +» 2 eee es ems 


Gas side heat transfer coefficient 


Gas temperature, °F. .... «es 


Steam inlet temperature, °F... . 


Steam outlet temperature, °F... 


Mem PLreSSure . . « « « © © © 6 
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Startup of the chamber will require some diluent other 
than steam. If excess hydrogen is used, flammability limits 
will restrict flame temperature to no lower than 1940°F (24). 
This may limit turbine speed during the startup period, or 
require diluent water spray and discharge of the mixture to 
the condenser, bypassing the turbine during the startup sequence. 
Fortunately, Edse's work shows virtually no pressure effect on 
adiabatic laminar flame velocity with very rich mixtures of 
undiluted hydrogen and oxygen. The effect of pressure quenching 
distance then should only occur Papo tne Po term in equation 
B-7, and the pore size selected for the porous plug under normal 
operating conditions should be satisfactory in preventing flash- 
back during the startup condition. Ignition may be provided by 
a spark plug or plugs suitably located or by some gas igniter 
such as that described for the liquid propellant rocket-type 
combustion chamber in Appendix A. 

In the throttling process, starting at maximum flow rate, 
the entire burner would be throttled to 75% of its maximum mass 
flow rate, after which a section of the burner would be cut out 
as the flow to the remaining sections was increased to maximum 
in as smooth a transition as possible. Subsequent sections of 
the burner would then be cut out in the same way until the 
design minimum flow rate is achieved in the pilot section. It 
is estimated that about 9 such sections would be required to 
Span the power range from 5 to 62.5 kw, assuming constant 
expander efficiency with mass flow rate. During the stepping 


process some incomplete combustion would probably occur. 
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Cooling steam to the section not in operation might be required 
to regulate the plug exit temperature at an acceptable value. 
Such cooling steam would require more thorough mixing down- 
Stream of the porous plug to produce a uniform chamber exit 
temperature. From a structural standpoint, the sectioning of 
the burner into several separate porous plugs would permit the 
use of thinner porous plugs than for one unsectioned piug for 
the entire burner. 

Control of the combustion chamber exit conditions must be 
accomplished through sensing of Several parameters: 

1) Chamber pressure 

2) Chamber exit temperature 

3) Reactant and cunente (feed) water flow rates 

4) Stoichiometry of any unreacted hydrogen and oxygen, 

probably best accomplished in the condenser 

One may select hydrogen as a governing flow rate, adjusted 
by variation in combustion chamber pressure from some preset 
value, the alaska pressure. The oxygen flow rate would be 
coupled with the hydrogen flow rate, maintaining a stoichio- 
metric ratio and operated with feedback from sensing devices 
in the condenser of the plant. Such sensing devices would 
determine any departure from stoichiometric conditions in the 
condenser, i.e., an excess of hydrogen or oxygen, allowing for 
Other than 100% combustion efficiency. Combustion chamber 
exit temperature would be sensed by a thermocouple and would 
be compared to a preset signal providing feedback to a basic 


diluent water/hydrogen mass flow rate ratio. Stepping of the 
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burner area to produce throttling could be accomplished by 
coupling with the hydrogen flow rate. The maintenance of a 
constant dilution ratio of steam to hydrogen and oxygen in the 
premixed gases can be achieved by coupling premixed steam flow 


with the hydrogen flow rate. 


Assessment 

The porous plug burner with its low operating temperature, 
appears to be a possible alternative to undiluted stoichio- 
metric burning in a rocket-type combustion chamber. Although 
the chamber outlet temperature selected for the feasibility 
design was fairly high, any desired lesser temperature could be 
achieved simply by the addition of increased amounts of diluent 
Steam downstream of the superheating coil. This could be 
accomplished by an evaporating coil placed downstream of the 
porous plug in the stream of the combustion gases or in the 
chamber walls, since addition of the diluent as a water spray 
would unduly increase the size of the combustion chamber. The 
throttling problem is not simple but appears to be solvable. 

The determination of quenching distance for the design 
hydrogen-oxygen-steam mixture is critical for the determination 
Of pore size in the sintered material. Possible alternatives 
to reduce porous plug pressure drop while at the same time 
Maintaining sufficient effective thermal conductivity should 
be explored. The exit temperature of the porous plug should 
be kept as low as possible to reduce possible plugging due to 


Peat zation of the nickel. The use of other materials than 
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nickel does not seem ee since other materials do not 
possess a high resistance to oxidation and at the same time 

a reasonably satisfactory thermal conductivity. Copper would 
not be suitable for the proposed application because of its 
corrosion related temperature limitations. Silver would not 


be suitable for the same reason. 
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ELot OF SYMBOLS 


quenching @iseance, EE or i 

activation energy for a chemical reaction Btu/lb mole 
or Kceal/g mole 

Speciive heat, Seu tom 

mass flux, lbm/hr fFt- 

enthalpy, Btu/lbm 

thermal conductivity of the porous plug material, 
Beu/hnr £e° R 

effective porous plug thermal conductivity, Btu/hr ft°R 
average thermal conductivity of the premixed gases, 
BEu/hr Et°R 

reaction rate constants for the hydrogen-oxygen 
reaction in the notation of Lewis and Von Elbe 

order of a chemical reaction 

Pressure, 1bf/in- | 

heat flux, Btu/hr £t* 

heat flux removed from the flame to lower the flame 


temperature from tT. to T,. (in excess of heat flux 


f 
necessary to heat the incoming premixed gases) 
Universal gas constant, 1.9860 Btu/°R lb mole or 
cal/°K-g mole 

Temperature, °R or °K, as appropriate 

flame velocity of the premixed gases 


pressure exponent for adiabatic laminar flame velocity 


Ot axialy distance through the porous plug, ft 
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distance 
plug; EE 
density, 
ignition 
fraction 


material 


of the flame front downstream of the porous 
Gr cm 

Dom/ ie" 

lag, ms 


of open area in a cross section of the plug 


SUBSCRIPTS 


refers to adiabatic conditions 


refers to actual flame conditions 


refers to the unburned premixed condition upstream 


of the porous plug 
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APPENDIX C 


Re-Entry Turbine Analysis 


General 

The single disk, axial impulse, multistage re-entry tur- 
bine has many of the attributes of the axial impulse partial 
admission turbine, but generally has higher efficiencies at 
lower specific speeds. Figure C-l is a schematic diagram of 
such a turbine. Its inherent simplicity (a single disk) often 
makes it preferable to a series of partial admission stages. 
Its major disadvantage is interstage leakage which can at best 
be minimized. As might be expected, the impact of this leakage 
is worse for turbines with smaller blade heights, simply be- 
cause the tip and axial clearances cannot be set or reliably 
maintained at values much less than .004-.005". The effect of 
leakage on efficiency is to lower the efficiency without leak- 
age, by a multiplicative factor. This factor may be simply 
defined as the mass flow rate which actually does work to the 
total mass flow rate passing through the multistage device. 

The losses, other than normal turbine nozzle and blade 
losses, of re-entry turbines are summarized below. Many of 
these losses are common to partial admission turbines. The 
terminology is that of Baljé, Linhardt and Silvern (1), (2), 


and (3). 
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CUTAWAY SCHEMATIC DIAGRAM OF A THREE 


FL GUIR Es Ge 


STAGE RE-ENTRY TURBINE WITH CROSSOVER BUCS 





a) Scavenging losses (the pumping 
and mixing loss of Stenning (4)) 
b) Disc friction losses 
Common to 
c) Blade pumping losses 
partial admission 
dad) Tip leakage losses 
turbines 
e) Effective reduction of the rotor 
velocity coefficient by Stenning's 
filling and emptying loss factor (4 
£f) A speed dependent, "dynamic", leakage 
Unique to 
g) Interstage leakage composed of radial 
FeEsen ery 
tangential and tip leakage components 
turbines 
h) Return duct losses 
Previous Work 
Baljé and Silvern (5) published a preliminary analysis of 
multistage axial impulse re-entry turbines in 1958. This work 
was summarized in a classical paper on turbine optimization 
published by Baljé (1) in 1962. A two-stage, 300:1 pressure 
ratio, re-entry turbine was built and tested by Sundstrand 
Corporation (6) and the results published in the open literature 
by Linhardt and Silvern (3). An extension of the analysis 
developed in the design of the two-stage turbine was published 
PyeLinhardt (2) in 1962. 
References 2 and 3 provided a good approach to the analysis 
of interstage leakage but neglected interstage leakage through 


gaps on the the downstream side of the rotor. Baljé (7) states 


that the leakage assessments of the earlier papers cited were 
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Bataimictic and that the leakage occurring is in fact quite 
complicated. 

The optimum pressure ratio split for the individual stages 
is not easily analyzed. Linhardt (6) approached this problem 
by a design analysis for two cases, that of equal isentropic 
Work per stage and that of equal pressure ratio. He found that 
the equal pressure ratio case afforded a significant improvement 
in efficiency (~ten efficiency points over a wide range of 
specific speed) for a two-stage re-entry turbine. 

Re-entry duct losses were also analyzed in reference 6 and 
compared with experimental results. For the re-entry duct 
geometry utilized (see Figure C-2), the duct pressure loss was 
found to be less than is commonly assumed (1), i.e. that the 
dynamic head of the previous stage is completely lost in wall 
MePction and turning. 

There does not seem to be much information on re-entry 
turbines of more than two stages which have actually been built. 
Only one four-stage machine was found in the literature (8). 
Baljé (7) is unaware of any re-entry turbines of three or more 
Stages which have been built using the design method developed 
by him and his colleagues. 

Reference 8 reports a four-stage axial impulse re-entry 
turbine of unequal pressure ratios, with an overall pressure 
ratio of 55.7 and a speed of 24,000 rpm. In a probable attempt 
to reduce the interstage leakage from the first stage, the first 
stage pressure ratio was 4.51 compared to pressure ratios of 


fee, 1.63, and 3.14 for the second, third and fourth stages 
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respectively. The high pressure ratio fourth stage may have 
been an attempt to reduce Reynolds number losses in the last 
stage. The design point total to static efficiency achieved 
was 0.432 for the 1/8" blade height machine in comparison to 
the predicted efficiency of .52. The analysis of the design 
showed an appreciation for interstage leakage but avparently 
was done without the knowledge of Baljé's work. As might be 
expected from the high first stage pressure ratio, a super- 
sonic relative inlet velocity to the rotor was ee eee 
there, while that of the second, third, and fourth stages was 
subsonic. This design abc cas apparently was used in an 
attempt to reduce the static pressure in the first stage rotor 
Cavity and hence interstage leakage into subsequent stages. 
However, the predicted interstage leakage as seen by the first 
stage analysis obviously does not take into account the signi- 
ficant tangential leakage on the inlet side of the rotor des- 
cribed first in references 5 and 6. 

It is apparent from Baljé's analysis and the test results 
of reference 8 that there is a trade-off due to reduced effici- 
ency because of partial admission filling and emptying losses 
and reduced stage flow rate because of interstage leakage. It 
is noted as well, that the blade design of this turbine could 
probably be improved, and that elimination of the supersonic 
First stage conditions would probably improve the overall total 
to static efficiency. A significant cause of the low efficiency 
of this machine is its small blade height and relatively large 


Clearances resulting in large values of tip and interstage leakage. 
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Design Analysis 


With the assumption that a turbine with equal pressure 
ratios for each stage represents a good design, although per- 
haps not the optimum pressure split, one may use the basic 
analysis of Linhardt (2) to design a multistage turbine. Cer- 
tain of Linhardt's assumptions must, however, be reconsidered 
and refined, particularly those involving interstage leakage. 
Seeemight simply present in this thesis only the refinements 
made to Linhardt's analysis. But in the interest of clarity 
and continuity, Linhardt's equations are included as well. 

In the terminology. of Baljé and Linhardt specific speed, 
N. is defined as 

NYQ. 
i = 7) (C-1) 


ad 


= 
=! 
¢) 
i 
c¢) 
Z 
tl 


speed in rpm, 
Q,= exit volume flow rate, ft 7sec, 


and He = isentropic total to static head, ft lbf/lbm, 


d 


and specific diameter, Dos another similarity parameter intro- 


duced by Baljé, is defined as 


1 
DH.” 
D = Las (C-2) 
¥Q 
3 
where D = overall diameter of the rotor including blades (tip 
io tip). 


A ratio of blade speed to isentropic "spouting" velocity, 


C,, may be then determined from equations C-1 and C-2 as 
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where = ¥29 Hog . 
The ratio of the arcs of admissions of the various stages 
may be obtained from the assumption of equal pressure ratios 
and an assumed polytropic exponent, A, related to the polytro- 


pic efficiency, a7 by the equation 


=: : 
= T= CRIT a 


The polytropic efficiency is related to the isentropic effici- 
ency by the relation 


ee ical 
Pie. CO Dihats = eran ene; 
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in Pro (Cae) 
where k = ratio of specific heats, 
n= the polytropic or small stage efficiency, assumed 
constant for all stages, 
Np-g = total to static isentropic efficiency for the 
turbine, 
and PY, = the overall total to static pressure ratio for 
the turbine. 
Then the ratio of the density of the fluid entering the 
cml Stage to the density corresponding to the total pressure 


entering the first stage is 


PAIS, 





7 (C-6) 








On 1 Pol 
and the ratio lengths of the arcs of admission (ay, Ass etc.) 
become 
a D _ pal 
=e ca ps ey (C-7) 
m Oya 


neglecting interstage leakage. Another useful relationship is 


the area ratio for a nozzle, A/A*, 


k+1 1 2 ke+1, 31 
ea) ao -—— _ op (=) 42 7 
= | [ (py) k lat 5 1S leer A [Pr. na Pry gues (C-8) 
where A = the exit area of a nozzle, 
A* = the related throat area, 
and n = the number of stages. 


The stage temperature ratio can be conveniently represented 
by the relation, 
aly k-1 


30 = Loe | = 
not Mp Ge RI, (C9) 





where Ti. = the inlet total temperature to stage m 
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-1 the inlet total temperature to stage m-l. 
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Here the exit static temperature of a stage and the exit total 


temperature are assumed equivalent. 
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The ratio of the blade spate to the mean spouting velocieys 


u/co, may be evaluated by considering the relation 





5s. oie 
o) m=] Oo 
where 
on = the overall spouting velocity, corresponding to 
the overall pressure ratio, 
co = the mean spouting velocity for the individual 
| stages, 
and som the spouting velocity for the me stage. 


Two convenient definitions of terms are 


1 k-1 
Yo ee es ) k (C-11) 
O 
and 
1 Ie. 
a ee ag! k (C-12) 


Now, using equations C-9 through C-12 in conjunction with 


the well Known isentropic flow relationship 


ce 2che cet Pre’ * (C-13) 
and 
/ k i k=1 
“o,m a AS 2h oa igo K-11" (px?) See (C-14) 
where Pret and | represent total inlet pressure and specific 
th 


volume to the m stage (assumed equivalent with static exit 


conditions of the previous stage). Noting the following 
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relationship, 





one may then derive the equation below through use of the 


geometric series formula: 








- . oe n/2 
Bel vivdh ® = ee eee ea oe eee | 
Sa (Fe Gy = a(S)? [+12 yy, 
re) Yo m=1 m1 Yo 1-(I-n py 4) * 
(C-16) 
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. Yi y= ae 
y = Ce { ———_—]} ~. (C-17) 
O° 1- (1-npy4)* 


Then the ratio of blade speed to average spouting velocity 
for a re-entry turbine with equal pressure ratio in each stage 


is as follows: 


ay 
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This ratio of blade speed may then be introduced into the 
expression for hydraulic efficiency developed by Baljé (1) for 
partial admission turbines, using the velocity diagram conven- 


Pron in Figure C-3. 
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FIGURE C-3 


Velocity Diagram 


Hydraulic efficiency, is the total to static efficiency 


le 
neglecting leakage, disk friction and partial admission losses 


other than Stenning's filling and emptying loss factor. 
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vay = rotor velocity coefficient W3/W. (Figure C-3) 
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R Z w2 
by = nozzle velocity coefficient, e5/ ce 
M9 = relative Mach number of the jet entering the 


rotor blading 
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The first term in the expression for rotor velocity co- 
efficient is effectively a fit of the blade loss factor of 
Soderberg and others (Figure 3.13 of reference 9). The second 
factor corrects for aspect ratio and the third for Mach number 
losses when the inlet Mach number of the fluid relative to 
rotor blades exceeds Uniey 

t n a 


=—— 2 ahs in equation C-19 is the 
Zia wl a 
1 m=1 m 


The term, (l- 


filling and emptying loss factor of Stenning (4) modified to 
include the ratio of blade pitch, t, to twice the average arc 
@eeeamission for the n stage turbine. It is noteworthy that the 
term, (1-t/2a), is not quite in agreement with Stenning's 
analysis, where he predicted the factor to be (1-t/3a). Baljé 
(5) states that the larger loss is in agreement with test data. 
Equation C-19 must now be modified by the other partial 
admission losses and re-entry losses. These are applied in 
two ways. First those which involve leakage are assumed to be 
mutually independent and are applied as deductions to the power 
developed by use of a mass flow rate deficit. Thus the hydraulic 


efficiency, equation C-19, is multiplied by a factor 


(-or rad ~ bn,tan ~ Su,tip ~ Su,dyn) (C-20) 


284 





and the ratio of actual power developed to ideal power available 
is 


Th (1- Sr rad a “r,,tan - St,,tip 7 Tae (C-21) 


This factor indicates the fraction of total mass flow rate seen 
by the turbine rotor which actually does work. The separate 
leakage factors are explained in detail in subsequent paragraphs. 

mnie second correction factor, common to all partial ad= 
mission turbines, represents the power developed in the turbine 
which is dissipated in doing work other than shaft work. Thus 
these factors represent power losses which are subtracted from 
the power developed to give the resultant shaft power. 

The turbine total to static efficiency then becomes as 


follows: 


=s n, (- SL rad Sr,tan °L,tip~ Sh chen ong 35 ote 


(C=22) 
where Cp = a disk £riction Lactor; the portion of the power 
developed which is dissipated in disk friction 
Sp = a blade pumping loss, first described by Stodola (10) 
Sot Scavenging loss coefficient, due to mixing of the 


high velocity steam with stagnant fluid trapped in 
the blade passages. 
The disk- friction coefficient includes the effect of 
Spacing of the disk and the surrounding walls (diaphragms) 
and the Reynolds number dependence due to disk diameter and 


local fluid conditions (5). The blade pumping coefficient 
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represents power dissipated by drawing of the stagnant fluid in 
the rotor cavity into the blade passages in the sector of the 
blading which is not active (unadmitted arc). The scavenging 
loss, analyzed by Stenning (4), is simply the fraction of work 
done on the rotor which must be expended to pump out stagnant 
fluid in a blade passage as the blade passage enters the active 
region. 

mie derivation of the equations for disk friction blade 
pumping, and scavenging is discussed in references (2), (4), 
(5), and (10). In the terminology of reference (2), they are 


presented below: 
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Xx = blade density factor, usually about 0.6; 
3 2 el 
bp = 2K V2g, (re 8 reqs)? )? -=e yg], (C-24) 
where 
5 = a pumping loss coefficient, 1.4 - 2.9 x oe 
and I, a (oa = 1) ee als 
and NOD. 
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where 
Ky = Gisk £riction cocificient, a function of Reynolds 
moamoer based on disk diameter and fluid conditions in the rotor 


cavity, given in references 1 and 5. 


Leakage Losses 

There are several paths which the working fluid may take 
in an axial impulse re-entry turbine other than through the 
rotor blading. Some of these paths will result in a complete 
loss of the working fluid. In other instances the path may be 
such that only one or more of the several stages will be 
affected by the leakage. 

Tip leakage, present as well in full admission and partial 
admission turbines, is fairly well understood and is predic- 
table. In the re-entry turbine, however, tip leakage may short 
circuit a particular stage. 

A diagram of the possible leakage paths in a re-entry 
turbine using the cross-over duct geometry is shown in Figure 
e-4, 

In addition to the leakage paths shown in the figure is a 
speed dependent, dynamic leakage. Dynamic leakage is caused 
by. that volume of fluid which enters the blade passage yet which 
is prevented from leaving the blade passage in the normal manner 
through the re-entry duct simply because the blade moves into 
the sector of admission for the succeeding stage. Earlier 


analyses (2), (3), (6), neglect any work done by this leakage. 
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Baljé (7) estimates that one nalf of this leakage actually 
eecomplishes work. In the analysis used in this thesis, it 
was assumed that no work was performed by dynamic leakage, which 


is conservative. 
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FIGURE C-4 


Leakage Paths 


Radial leakage may occur through the gaps on either side of 
the rotor. A notable discrepancy in the analyses of references 
2, 3, and 6 is a failure to include leakage on the downstream 
Side. It is evident that no accurate leakage assessment may be 
obtained without inclusion of leakage paths on both the inlet 
and exit sides of the rotor. 

In assessing the radial leakage on the upstream side of 
the rotor, it is easiest to look first at the leakage experi- 
enced by the first stage. Assuming choked flow conditions 
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between stages, the ratio of the leakage flow to the inlet 


Blow to the turbine, m may be expressed through use of the 


a Hes 
choked flow relation 
m RTO 
Pay = constant (C-26) 
Oo 


where Por Ty refer to total conditions and A* is the area of 
the leakage gap or nozzle throat as applicable. 
Using equation C-26, the radial leakage factor for leakage 


upstream of the rotor in the first stage is as follows: 


e -l1/n 

My xad,u,l _ Ce,r *1?st,1 91 _ Ce,r 81° a, 
@ a p A de 
my the OL le" ay Sina, A*/A 

a s,/D Pr 1/n pps 
- = Ky ; (C-27) 

h/D Sina, 
where 


 rad,u,1 = leakage mass flow rate in the radial direction, 


upstream of the rotor, in the lst stage, 


ain = the nozzle throat area, 


a the radial leakage gap upstream of the rotor, 


Po.1= total pressure entering the first stage, 
7 


Per 1 static pressure in the first stage cavity, 


n = number of stages, 
and Cp. =a radiil restriction factor which can be adjusted 
a 


to a value of ~1/3 with proper seal geometry (2). 
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Here it is assumed tnat static pressure in the first stage 
motor cavity is roughly equivalent to the total pressure in the 
leakage gap and that the total temperature existing in the leak- 
age gap is equal to the total temperature at the inlet to the 
Stage. This is contrary to the assumption of references 2, 3, 
and 6 which state that the stage exit static temperature exists 
at this point. This, however, is a small discrepancy since the 
Mieetstage temperature ratio is nearly unity for reasonable 
stage pressure ratios. 

It may be shown as well that the radial leakage on the 
upstream side of the rotor for the second stage is related to 
the mass flow through the nozzle of the second stage by the same 


actor, K Yet when the leakage is compared to the total inlet 


L° 
flow rate, the ratio may be shown to be 


mo ,rad,u,2 


m 


= K, (1-K))- (C=2Z 2) 
a 

The total leakage seen by the second stage is the sum of 
that lost in the first stage and second stage. Thus the second 
stage experiences a total leakage factor of 


Mo rad,u,l ~  e,rad,u,2 


m 


= Ki + K, (1-K,).- (C-29) 


1 


Since there is no leakage from the ‘ae stage itself, the 


_ and uaemy) SS stages see a total leakage factor of 
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Ky + K, (1-K)) + K, (1-ky) Tiassa ats Ky (1-k ) (C=30) 

By use of the geometric series formula, the average radial 
leakage factor for all n stages, for leakage upstream of the 
rotor, may be shown to be 


=e ae a a = 1p ieee ae _ 
Gh radu 7 n (Mtltg ((1-2k)(1-K) ale (C-31) 


On the downstream side of the rotor, the analysis is similar, 
except that a factor for stage temperature ratio must be intro- 
duced into equation C-27. Additionally, it should be noted that 
the leakage here occurs after it has done work in the first 
Stage. Thus the first stage is unaffected by radial leakage at 
the exit side of the rotor. The last stage again experiences no 
radial leakage of its own but sees a mass flow rate deficit due 


to leakage in preceding stages. 


. =1/n 
_ m c S,/DiP rE A/A* 
Sr, rad,d,2 eke sey = fir 2 S z = K, 
™m Jay by Siatietens (lag si 5 )) 
a 2 pe 2: (c-32) 
S5 = axial clearance at the exit of the rotor 


Equation C-32 is the radial leakage factor for first stage 
leakage downstream of the rotor, which affects the second stage. 
The third stage experiences a leakage factor which repre- 


Sents leakage occurring at the exit of the first and second 


Stages, 
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s 
m = P es 
L v IgGl 7 a w 1 l raG ie > 


‘L,rad,a,3 ~ Se a ee 


and the leakage factor for the ae stage is 


ey i-2 
Bmesd,dm =, “gi Ka) 


The average radial leakage factor for all n stages for leak- 
age at the exit side of the rotor may then be shown to be as 
follows: 


ee ee ee ” 
SL, rad,d a Kon [1- (1 Ky) }. (C-34) 


The total radial leakage factor, combining equations C-31l 


ma C—-34, becomes 


= ntl, zoe [(1-2K,) (1-K 


> n-1l if 
SL rad n Kjn 1} +1 


) 
1 Kon 





n 

Pes 

2 K,) ine 
(C=35) 

Tangential leakage may also be conveniently separated 

gnto components at the inlet and exit side of the rotor. Again 

modeling the leakage path as a simple gap between stationary 

walls, one may develop equations similar to those for radial 

leakage. Reference 6 notes, however, that leakage at the inlet 

tan,l ee 

Figure C-4, exceeds that predicted when the assumption is made 


Men oL the rotor in the direction of rotor rotation, L 


that the total pressure at the inlet to the leakage gap is the 


Static pressure in the rotor cavity. Baljé (5) initially 


Zoe 


assumed that this pressure exceeded the static pressure by a 
factor of V2. For higher pressure ratio tUmbines (3), (6);5 
this assumption was found to be inadequate. Reference 6 pro- 
poses that the tangential component of the nozzle exit velocity 
be converted into its equivalent stagnation pressure in order to 
give a more realistic assessment of the total pressure conditions 
ahead of this gap. 

This leads to the following equation for the leakage fac- 


mor through path L for the first stage, assuming the total 


tan,1l 


temperatures ahead of the nozzle and leakage gap to be equivalent: 


Ss. A 


ie capepacnly le Ce ytane, ane 


o1,tan 
’ ,path 1,1 : ; A 
m4 ay Sina, h A*/A Poval 


* 
> tan F s,/D A/A 


a,/D Sino 


= R) (C=26) 


2 


where the last subscript on m designates the stage 


ean 


affected, 
k-1 2 2 2 
= O,l,tan on 2 aN yas eo we = 
m= r.,  ' K-12 ieee 
He eae oe 
M, = the Mach number of the fluid leaving the nozzle 
f 
under isentropic conditions, 
and Ce tan = a tangential restriction factor, assumed to be 
’ 


about 2/73 (3), (6). 
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The tangential leakage through path | 


but augments the mass flow in the 
is reasonably small in comparison 
effect on the second stage nozzle 

The leakage through the same 
be expressed as a fraction of the 


using equation C-7 yielding 


m . 
Pycan,path 1,2 ee 
il! Oo 


ay) 


el 1s not all lost 


second stage. If the leakage 


to the nozzle mass flow, the 
jet may be considered minimal. 


gap in the second stage may 


second stage nozzle flow 


(C-37) 


which, when expressed as a fraction of the total inlet flow to 


the turbine, m becomes 


aie 


or,tan,path 1,2 iy 
i 


- my ,tan,path 1,2 Sea -1/na (1-R.) 
‘ 1 Oo ‘a 


(C=3i8)) 


A similar expression can be shown to exist for succeeding 


tele 


Stages. For the m 


SL.,tan,path 1,m al Pro 


m-1 
ae, 


stage, the leakage factor is 


US oem ole lene 
(C-39) 


In the first stage a tangential leakage exists upstream of 


Pie rotor in the direction opposite to rotor rotation. 


This 


leakage is lost to all stages and may be expressed as follows: 


Ce tan an Se 


SL,tan,path 2,1 ~ 
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a,/D Sina 


~I/n s,/D A/A* 


= R 


3 (C-40) 


2 





The average tangential leakage upstream of the rotor com- 


Persing paths L and L may then be expressed as 


tare tan 
a 1 iam - (824) 
Sr,tanl,tan2 7 a ae a Bhs as PE es, Seen eth emeite! 2s 


(C-41) 
On the downstream side of the rotor tangential leakage 


Occurs through paths L and L 


tan,3 eetin ly The first stage is not 


affected by this leakage, as was the case for radial leakage at 
the exit side of the rotor. The leakage on the downstream side 


of the first stage in path L bypasses the second stage 


tan, 3 


nozzle, and the ratio of this leakage to the mass flow rate 


entering the first stage is as follows: 


: = ln A 
: path gel” Seen eat tc shia’. 
fptan,path 3,2 — . 7 = + 3 
m, a,/D Sina. (1 De y;) 
(C-42) 


Similarly, the leakage seen by the third stage becomes 


m 3-2 
Pe)tan,path 3,3 * Po 
1 


- sl Na = 
= R, PY. Ct R3), (C-43) 


and for the mo stage, 


m m-2 
_ _&,path 3,m-1 a -(——) » 
St,,tan,path S37, it ° 3 oe Bn 
m 
1 
Le (C-44) 


(SG lean eeten 3,m~-2 
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Additionally, there is a tangential leakage at the exit of 
mie rotor, in the first stage only, directed opposite to rotor 


rotation. This leakage, L reduces the mass flow 


tan,path 4’ 
rate in all the stages after the first stage and is equivalent 
to the leakage in the direction of rotor rotation for the first 
stage, equation C-42. 


Thus the average tangential leakage factor for the exit 


Side of the rotor is as follows: 


n 

= IL : 

SL,tan 3, tan 4. ee [(n-1)R, uJ He $t,,tan,path Saal” (C45) 
where 

~ (Mz 2y 

ST,,tan,path seme: Se ae ys ASE hem, pee Semel 
and 

C = 0. 

im,tan,path 3,1 

Tip leakage through path es 1 is the normal tip leakage 


from the "pressure" to the "suction side" of the blade experi- 
enced in axial turbines. Baljé (5) expresses this leakage as 


a fraction of the total mass flow as follows: 


i> 5: ae 
if tip,path 1 Lar h/s 3 ee h/D 
s,/D 
where S3 = tip cle irance 
and C/t = chord/pitch ratio, assumed to equal 2 Sin 28. (a 


particular geometry found by Baljé <= be adequate -see Fig. C-5). 
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PACGURE 1€-5 


Turbine Blade Geometry 


Two other tip leakage paths exist and are unique to re- 


entry turbines. One of these paths,L exists only in the 


tip,;. 
first stage, where leakage occurs over the blade tips into the 


exhaust of the last stage. The other path, L exists in 


tip, 2 
all stages but the last. It, however, is not completely lost, 


in that it re-enters the exhaust duct of the succeeding stage 
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and actually does work in the second stage downstream from 
miter it occurs. 

Using an analysis similar to that for tangential leakage, 
the tip leakage through path 3 for the first stage may be 
expressed as a fraction of the inlet mass flow rate to the 


turbine: 


mo, tip,path oly £, cup = Pst 1 = 
my ha, eee iL 


=l/n a*/D CosB.A/A* 


4c S,/D Pr, 


eee (C-47) 


h/D ae Sino, 


where c =a tip restriction factor, probably about 


12 pe ealge) 
equivalent to ean 
and a* = cutter diameter or rotor blade throat width (Fig. 
C= 5). 
The other leakage path for tip leakage unique to re-entry 


turbines L may be analyzed in the following manner. The 


low 
furst Stage experiences a leakage out of its normal working 
path and into the exhaust duct of the second stage. Thus the 
first and second stages are affected by this leakage. The 
second stage also experiences leakage over the tips of the 
blades into the exhaust duct of the third stage. Thus any 
intermediate stage experiences tip leakage through path 2 in 
two manners, (a) bypassing of the stage nozzle by incoming 


tip leakage and (b) bypassing of the stage rotor by tip 
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Heakage into the exhaust of the succeeding stage. 
The leakage experienced by tip leakage by path 2 in the 


first stage is as follows: 


Bee eae sh Onl 


$1,tip,path Zn : 
m 
1 
4c . $.,4/D F a*/D CosB 
eee Lp. 3 2 E, (C-48) 


h/D a,/D Sina, 
The second stage experiences this leakage plus its own 
leakage into the third’ stage exhaust. 


il 


E, + E PY, nr (1-E,) (C-49) 


o1,tip,path 2,2 “1 1 


The average leakage factor for all n stages for tip leak- 


age paths 2 and 3 becomes 


i 
= ‘ ll 
ot,,tip,path 2-) path 3 5 a {E, a aoe iO, tip,path2 ,m-1 


! (C=50)) 


y OT tip,path 25m i ST fee poke) aie” 


ae Ce cipe earn Zk sy =a 


mel 


and Q = 10 ee nar ) (1 - Q 


b,tcip,;path 2,;m 1 moe 
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Bor the conditions of interest in this thesis, 1.625) low 
pressure ratios, tip leakage through paths 2eand 3 didenoer 
significantly affect the efficiency and was ignored in the 
computer analysis. 

Dynamic leakage may be represented by the following 
expression (2), (6), assuming that the dynamic leakage does 


mo work: 


Neo 


CD x (C/D) Ny De 


SL,dyn ~ Wy (a, 7D) Sino 154 


Z 


II (C-51) 


where Cy = restriction factor, based on exhaust geometry and 


nextstage inlet pressure, probably about 0.35 (6), and 


ne 
1-Pr, An 
I, 2 = 5 ¥y 
ike a 
re) 
Summary 


The model used to analyze leakage in this thesis is that 
used in references 2, 3, 5 and 6. Nowhere in the actual tur- 
bine, however, is there a leakage flow through a gap with 
Stationary walls. The model is probably best then for radial 
leakage, which most closely approximates this condition. A 
possible alternative approach to leakage might be to model the 
tangential and tip leakage unique to re-entry turbines as that 


which occurs through a series of labyrinths comprised of blade 


tip, trailing or leading edge and the stationary wall. Thus 


the leakage would be dependent upon the total number of 
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Mapyrinths in a leakage path and the interstage or totgi 
pressure ratio (depending upon the path) rather than the static 
or total pressure within the rotor cavity of the stage in 
question. Such a model would be speed dependent, since the 
number of "labyrinths" in a leakage path would vary with speed. 
Perhaps this approach might yield a more precise agreement 

with experimental results. It is evident, however, that more 
experimentation is necessary with re-entry turbines of three 

and more stages before any leakage model can be shown to be much 


better than another. 


301 





IFO . 


REFERENCES 


Baljé, O.E., "A Study on Design Criteria and Matching of 
Turbomachines: Part A-Similarity Relations and Design 
Criteria for Turbines", Journal of Engtneerning for Power, 
Paw o4, NOwil, January 19627, 0p S32102. 


Linhardt, H.D., "Re-entry Turbines for Space Power Systems", 
mes JOuunal, Vol. 32, No. 10, October 1962, pp 1152-1560. 


iiennardt, H.D. and D.H. Silvern, "Analysis of Partial 
Admission Axial Impulse Turbines", ARS Journal, Vol. 31, 
ooo, Maren 19617 pp 297-306. 


Stenning, A.H., Design of Turbines for High-Energy-Fuel 
Low-Powen-Output Applications, Dynamic Analysis and 
Bontrol Laboratory Report No. 79, Massachusetts Institute 
of Technology, September 30, 1953. 


ewije, O.E. and D.H. Silvern, A Study of High Energy Levelt 
Low Power Output Turbines, AMF/TD No. 1196, Department of 
the Navy, Office of Naval Research, Contract No. NONR -. 
wee (00). Task No. NR 0944343., 9 April 1958 (AD 161323). 


Linhardt, H.D., Study of Turbine and Turxrbopump Desagn 
Parameters, Finak Report-Volume 1, A Study of High Pressure 
Ratio Re-Entry Turbines, Department of the Navy, Office of 
Naval Research, Contract No. NONR-2292 (00), Task No. 

NR 094-343, 30 January 1960 (AD 232635). 


Baljé, O.E., personal communication, March 8, 1972. 

Mong, Robert Y., D.L. Darmstadt and D.E. Monroe, Invedatiga- 
tion of a 4.0-Inch-Mean-Diameter Four Stage Re-entry 
Turbine for Auxikiany Power Drives, NASA TM X-152, 1960. 


Merlock, J.H., Axial Flow Turbines, London, Butterworth 
and Company, 1966. 


Stodola, A., Steam and Gas Turbines, New York, Peter Smith, 
1945. 


S02 





LEST OF SYMBOLS 


are of admission, £t. 
mozzile thteais area, ft: 
cutter diameter (throat width of rotor blade passage), ft 
nozzle area ratio 
motor blade chord length, £t. 
absolute velocity ft/sec 
Bestriction factor 
tip to tip rotor diameter 
specific diameter, D Ba" /¥Oq: for the entire turbine 
leakage factor fae re-entry tip leakage for the first 
Seagc 
ratio of total pressure in the tangential direction up- 
stream of the rotor to total pressure ahead of the nozzle 
conversion factor, 32.2 ft lbm/lbf sec? 
blade height, ft. 
isentropic head, ft lbf/lbm 
ratio of specific heats 
radial leakage factor for first or second stage, or loss 
coefficient for power dissipating losses 
mass flow rate lbm/sec 
Mach number 
number of stages 
specific speed 25 for the entire turbine 
, a 7 


stage or total pressure ratio 


pressure, lbf/in? or lbf/ft? 
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tangential leakage factor for first to second stage 
clearance | 

axial clearance npseeean Of the rotor 

axial clearance downstream of the rotor 

tip clearance 

PplLade pitch, £t. 

blade speed, ft/sec 

welocity “relative Ee jehe blade, ft/sec 

pressure ratio function 

angle of absolute entering or leaving velocity, measured 
from the tangential direction 

relative blade angle, measured from the tangential 
direction 

loss factor 

efficiency 

polytropic exponent 

density of working fluid lbm/ft? 

blade density, ratio of open volume to blade volume in 
the rotor blading 

velocity coefficient 

ratio of tip leakage in the direction of rotor rotation 


from a particular stage to total turbine mass flow rate 


304 





SUBSCRIPTS 


downstream of the rotor 

disk friction 

refers to dynamic leakage 
hydraulic 

refers to stage conditions 
refers to leakage 

refers to general stage 
refers to total conditions or isentropic conditions 
refers to nozzle 

polytropic 

blade pumping 

~ refers to radial leakage 
refers to rotor 

refers to static conditions 
scavenging 

refers to total conditions 
total to static 

refers to tangential leakage 
refers to tip leakage 


upstream of the rotor 


305 


APPENDLX .o 


Condenser Analysis 


General 

Since the condensing system of the proposed propulsion 
plant must allow for the presence of non-condensables in sign- 
ificant concentrations, the condenser was analyzed in detail. 
The product of this investigation is a computer program capable 
of predicting surface condenser area, weight, volume and fluid 
exit conditions for specified condenser inlet conditions and 
geometry. The program has the capability of analyzing cases of 
up to three non-condensable gases, hydrogen, oxygen and nitro- 


gen, with superheated, saturated, or wet vapor steam. 


Theory 


The classical approach to condensation with the presence 
of significant quantities of non-condensables present is that 
of Colburn and Hougen (1). Colburn's approach uses the so- 
called Chilton-Colburn (2) analogy between heat and mass trans- 
fer. Votta (3) presented a more direct method of obtaining the 
condensation rate using the principles developed by Colburn and 
Hougen. Neither one of these approaches, however, is directly 
capable of analyzing condensation of superheated steam. In 
both methods, the solution is obtained by assuming a set of 
temperatures less than the temperature of the entering steam 
and greater than the bulk coolant inlet temperature, each temp- 
Beature corresponding to a specific location in the condenser. 
For saturated steam conditions, this fixes the bulk (free stream) 
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steam partial pressure at the point in the condenser where the 
assumed temperature exists. At this point one may determine 
the condensation rate by an iterative procedure. One 

assumes different vapor-ligquid (condensate) interface tempera- 
tures and corresponding vapor partial pressures until the 
energy transfer due to heat and mass transfer through the gas 
film matches the heat transfer to the bulk coolant. By Votta's 
method, the interface temperature is calculated directly. By 
the classical method, or by Votta's method, one risvaettee deter- 
Mine the condenser area required (based upon a previously 
selected condenser geometry) to condense enough steam to lower 
the steam partial pressure to that corresponding to the assumed 
bulk gas-vapor temperature. Using this method with appropriately 
selected temperatures, one may then estimate the required con- 
denser surface area. 

Neither the previously described classical iterative 
method nor the method of Votta is directly applicable to mix- 
tures of superheated steam and non-condensables. In the proposed 
application of this thesis it was anticipated that the state of 
the steam entering the condenser would be superheated in nearly 
every conceivable case, with between 200-300 degrees of super- 
heat likely. For this reason a method of analyzing the super- 
heated steam case was devised. 

It was also anticipated that various mixtures of non- 
condensable gases could exist in the condenser at any one time. 
Therefore, it was determined that the method devised should 


permit the analysis of condensation rates for the simultaneous 
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presence of both hydrogen and oxygen in the condenser (the case 
of other than 100% combustion Saeucioneny Since the condenser 
would be placed in some sort of pressure vessel, itself having 
an "atmosphere", provisions were made for the additional 
simultaneous presence of a third gas, nitrogen. Inclusion of 
nitrogen in the analysis would also permit testing of the 
method against existing data where air is the non-condensable. 

It was also determined that the method of analysis should 
include the most significant factors which affect the heat/mass 
transfer process. It is common, for instance, in the problem 
of condensation in the presence of a non-condensable gas, to 
neglect the condensate film heat transfer coefficient (4) as 
well as sensible heat transfer. The condensate film heat trans- 
fer coefficient becomes increasingly more important as the 
condensate layer is augmented by dripping from tubes above a 
particular row in question. Coolant film resistance to heat 
transfer as well as scale and tube wall resistance are also of 
importance and cannot be neglected. 

Certain other aspects of condenser design are not easily 
adapted to analytical treatment and are more the subject of 
detailed condenser layout and design. Among these are baffle 
Placement, steam lane location and the location of the suction 
for the removal of non-condensables. Each of these items has 
an effect on pressure drop. Condenser pressure drop is complex 
and is not well understood at present. Silver (5) and Griffith 
(6) have described the difference between the frictional pressure 


drop in condensers and that in normal heat transfer tube bundles. 
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This is essentially an effect similar to that of a diffuser 
since condensation delays boundary layer separation and thereby 
reduces form drag. Silver estimates that under normal conden- 
sing conditions, the frictional pressure drop in a condenser 
may be as low as 68% of the pressure drop of a non-condensing 
gas of the same properties. The small pressure differences 
fetoting in a condenser make experimental determination of 
Pressure variations difficult. 

Sebald (7) gives a method of determining a mean condenser 
pressure which uses friction data derived apparently from non- 
condensing situations. This method requires iteration of steam 
lane and bundle arrangements for minimum pressure drop between 
the periphery of the condenser tube bundle and the air ejector 
suction. 

While pressure drop is of significance in detailed conden- 
ser design, it is apparent that it is by no means as important 
as the effect of the presence of a significant amount of non- 
condensables. Reference 8 indicates that as little as 5% non- 
condensables by weight can result in an increase in condenser 
area of 300% over that of the pure vapor case. For this 
reason, pressure drop has been neglected in this analysis, but 
May be included, if desired, after the condenser surface area 


and basic geometry have been determined. 


Method of Analysis 


The best approach to the superheated steam case was 
determined to be a row-by-row analysis of the steam conditions 
in the condenser, assuming a uniform downward flow of steam 
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across a staagered tube bank, rectangular in lavout, with a 
Beare pitch rotated. Such a bank is shown in Figure D-l. In 
such an analysis it would be necessary to write conservation 
equations for a given tube row, assuming that, except for the 
first and second rows, a drip of condensate would occur onto 


the row in question from the second row above it. 


FIGURE D-1l 


Staggered Tube Bank with Square Pitch Rotated 


Pcontrol volume for a typical tube in the first or second 
rows is shown below. The dotted lines indicate the boundaries 


of the control volume. 
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Aak 
Ame tvib 
AL 
CONDENSATE 
INDICATES S--7-- FILM 
ee CONTROL. | AUite 
BOUNDARY 
FIGURE D-2 
Control Volume I 
Am. = condensation rate for a particular tube, lbm/hr 
A, = outside area of the tube per unit length, fies 7 ft 
A, = inside area of the tube per unit length, EEC7EE 
a = enthalpy, Btu/lbm 
i, ot enthalpy of the bulk vapor in the free stream Btu/lbm 
iy 
us = average enthalpy of the condensate film, Btu/lbm 
Iget = sensible heat per tube transferred to the condensate 
from the bulk vapor, Btu/hr 
. = sensible heat per tube transferred to the bulk 


coolant, Btu/hr 
From the basic Nusselt analysis (8) the following energy 


equation may be written for control volume I: 


e a Ss = eS = % G = _ 
Am {i- cee 4 Tol 2. q., oe Am. lob 0 (D-1) 


Side 





a = saturated liquid enthalpy at the gas vapor interface 
partial pressure, assuming no subcooling of the 
condensate at the interface 


= average condensate film specific heat, Btu/lbm°R 


pre 
T. = vapor-condensate interface temperature, °R 
TO an outside wall temperature, °R 
a 


From equation D-1l one may determine the condensation rate 
for a single tube, 


F ome | 
Am = ine ce a ree eo (D-2) 


Cc : ae a 3 
17,b egg costly ToL 
The heat transferred to the bulk coolant, qj, May be deter- 
Mined from the overall heat transfer coefficient, in this case 
selected on the basis of outside tube area, Uo and the log 
mean temperature difference of the condensate-vapor interface 


and the bulk coolant. 
f= U_ AT (D323) 


U_ = overall heat transfer coefficient based upon outside 


tube area, Btu/hr fE- °R 


a oe use e 

ATM = ar =, , the log mean temperature difference, °R 
i L,e 
Tey 
2 eye 


| 
ll 


bulk coolant temperature entering a tube row, °R 


wal? 





ae) 
I 


bulk coolant temperature leaving a tube row, °R 
ns = condensate-vapor interface temperature, °R 
The overall heat transfer coefficient, Use is obtained 
from the various resistances to heat transfer which comprise 


the path from interface to bulk coolant, 








a It & 
Bo Ca d din d_/d, Oe, 
Oo Oo fe) On iL 
d 4 . pee Ae 
a 1 Se W C 
d, = inside tube diameter, ft 
d, = outside tube diameter, ft 
hy = condensate heat transfer coefficient, Btu/hr ft? °R 
hy, = coolant film heat transfer coefficient, Btu/hr £e- eo 
Dee = coolant side scale heat transfer coefficient, 
Wey /oke sme> —1R 
k = thermal conductivity of the tube wall, Btu/hr it “kK 


The coolant film heat transfer is obtained from McAdams 


correlation: 


hed. God: li ¢ 
2 L Lyin EO | - Bel) 0.4 eas 
i D ULb if 





rate 
lI 


thermal conductivity ef the coolant, Stuynr £t F 


GCeoolane viscosity, Von/ire ce 
c = coolart specific heat, Btu/lbm°R 


and the subscript, b, denotes bulk properties. 
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The condensate film heat transfer coefficient, Aas 1s 
Obtained from Chen's correlation (8), modified to allow for 
Superheated vapor. For Control Volume I, Chen's correlation 


may be approximated by 


ee) 
9 02 ae 
oe 0.720 (= ee: (D-6) 
5 Eg te o) 
where 
G@e= gravitational constant, 4.17 x 10° £t/hr? 
e. = average density of the condensate film, lbm/ft? 
Re = average thermal conductivity of the condensate film, 
Bicwd/ rece ok 
OT) = See Le = * 
ae = ea + BCG tt Vo) oe T,) with Ty ,b the 
bulk vapor temperature, and e - the specific heat 
a 
of the condensing vapor based upon the arithmetic 
average of bulk row inlet temperature and condensate- 
vapor interface temperature 
Uo = the average condensate film viscosity, lbm/hr ft 


For subsequent tube rows (3, 4, 5, etc.) the control 
volume must be modified to include dripping of the condensate 
from previous rows. Figure D-3, Control Volume II, gives the 
conservation equation and the resulting condensation rate for 


the third and succeeding rows: 
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FLEURE D=3 


Control Volume II 


e P ia 2 - ° 2 _ 2 e - e 
dwn tim fig Bowe od eee Lele dof n (Ami. bln 0 
(Dam) 
Am = E 
c,n =. 5 SSS. ee 
fay pn leg Bey (TAT In) 
fa, nM nna te78S5 6 Ti wo! In 7 (m3) n-2 Agen! 
(D-8) 
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where 


Am. aoe condensation rate for an aoe row tube, lbm/hr 


f 


The condensate heat transfer coefficient for a particular 
row below the second row may then be approximated from Chen's 
correlation, which assumes, for the row in question (row n), 
that all previous rows exhibited the same interface temperature 
as row n. This assumption, of course, is not strictly correct 
Since the average condensate film temperature decreases with 
succeeding rows. Thus,this method predicts a somewhat higher 
condensation rate than would actually occur. However, since a 
real condenser would have baffles to remove condensate drip, 
this assumption is probably not overly optimistic. With this 
assumption, equation D-8 reduces to equation D-2. 

Another aspect of Chen's correlation is that in its form 
presented in reference 8, it estimates the average heat trans- 
fer coefficient for an n-row, in-line bank of tubes. This may 


be applied to the ace row of a staggered bank by the approxima- 


Lon: 
en r a a row bank ™ 7 WS ee row baa? (D-9) 
where 
n-l . 
m = a it pn odd 
n= > if n even 


Then, substitating Chen's correlation, one may write 


a6 





140" 2(r aa (p aa APES! au 


) aie 
ho, = mio. 728 [05 Ce ee 
“fg . oh Hee! 
Oe 2c ela eee) Gi Goals a1 1 
B(m-1) {0.728| Pe 4 Wio j;__¢ ¢ fg ____}*} 
EG (m-1)d. Oe ee 
(D-10) 


Computation of the mass transfer coefficient requires the 
use of the Chilton-Colburn analogy and certain tube bank heat 
transfer data. Tube bank data was obtained from Colburn's 1933 
data (9), based upon staggered tube banks data from as early as 
1911. There was some doubt as to the number of rows in the 
staggered tube bank data observed by Colburn. It was evident 
from investigation of the American data referred to by Colburn 
that each of these banks were four rows deep. This is of con- 
cern since data of Kays and Lo presented in McAdams (10) shows 
that for succeeding rows in a staggered bank up to the 6th row, 
the heat transfer coefficient improves. After the soe row 
there is no further improvement. Applying this information to 
data of Colburn, assumed to apply to banks of 4 staggered rows, 
One may then determine the heat transfer coefficient for a 
given row in an n-row staggered bank. Colburn's method of 


data presentation is given below in Figure D-4. 
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FIGURE D-4 


Colburn j Factor Vs. Reynolds Number 


ime term, j, 1s Colburn's j factor 





j = 9£ (Spam bom) # cca Bes (wn yt (D-11) 
p,vm vm max’ mean vm? 
where 
hoe megdas film heat transfer coefficient Btu/hr £t*°R 
Co ,vm = film average specific heat of the gas vapor 
Maxture, Btu7lbm°R 
 s = mass flux of the gas vapor mixture based on minimum 
€ress sectional area, Ibm/hr £t- 
om = film average thermal conductivity of the gas vapor 
mazture, Btu/hr ELoR 
Pot = eee , the log mean partial pressure 
Nee difference of the non-condensing 
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gases across the gas film: p-total pressure, 


IbEJ in oe pee ae partial pressure of the condensing 
Vu Vip) 


vapor at the condensate-vapor interface and bulk 


conditions respectively 


K, = mass transfer coefficient lb-mole/hr-ft?-lbf/in? 
D =NdiPritcion Cocenlotent mimasscmaLEEUSTY1ty tt. im 
oe, — gas mixture density, lbm/ft? 


= mean molecular weight of the gas vapor mixture 
in the free stream 

For heat and mass transfer, Colburn defines film average 
temperature as the arithmetic mean of the bulk and wall 
temperature. 

Some confusion does exist in the definition of minimum 
cross sectional area, from which Goa is determined. Through 
examination of Colburn's references it was determined that the 
Minimum area for a staggered bank is that shown in Figure D-5 


for the assumed square pitch rotated configuration. 









MINIMUM 


CROSS SECTIONAL 
AREA 
2xL2te 


S 
FIGURE D-5 
Tube Layout 
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he 1s evident from Figure D=5 that the ratzo of mMininum 
area to frontal area is Y2 - 2 d/s- 
Using equation D-11l, one may write an expression for the 


condensation rate Am, for a given tube, 


Am, = ALk Ka a (DF oe (D-12) 
where M. is the molecular weight of the condensing vapor 
(lbm/lb-mole), and the sensible heat transfer rate, 

ee = At hoe (Ty py Ty): (D-13) 


The determination of the properties of the gas vapor mixtures is 
discussed in a subsequent section. 

Using equations D-3, D-12, and D-13 one may then write 
for a single tube the condition which must be met for a 


solution of the heat/mass transfer problem. 


Ee ) 


Bea LAT, th AQ’ (TL Be CSS a aa 


LM ger oO me a 


ts G M. A,% ub 


(D-14) 
One then may determine the enthalpy of the condensate 
leaving a particular row. 


Sec) 1 oe ) 


ea 2 PAU ESS ooh My el a) 


,O 
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Mae determination of bulk vapor exit conditions fromed 
tube row requires the solution of the energy equation across 
still another control volume. 


Mve Lv.e y Mai Cert Lge, j 


cg . 
| 
| 
| 
| 


i 





ee Ytieg,j Ug, 2, i 
tA Ld FIGURE D-6 tA, £ 


Control Volume III 


r 
» oe 0 : e : © » oe 
Mg a ee eae te fae + ee eee! 
r se e 
ry mo . iL ot a = 0 (D-16) 
sy 9rd Gres) v,e v,e 
J 
where 
ma, @ = Vapor mass flow rate entering the tube row, lb/hr 
a 
mi g = vapor mass flow rate leaving the tube row, lbm/hr 
a 
m, 5 = mass flow rate of non-condensable gas j, lbm/hr 
a 
iL, Ba. We enthalpy of the bulk vapor entering/leaving 
a a 


the tube row, Btu/lbm 
1 a es ee = enthalpy of non-condensable gas j entering/ 
leaving the tube row, Btu/lbm 


Se. 


154 number of tubes in a tube row 


i number of distinct non-condensable gases present 
From equation D-13, one may derive an expression for the 
temperature of the bulk gas vapor mixture leaving the tube 


row, uSing the familiar relationship, 


ove 60fV,R eee y bv en be 
where 
T 5 the bulk gas-vapor temperature entering the tube 
v Ul ; 
row, assumed identical with Th ee of .equaetioneb—1 3. 
a 
es e me rt e 
- i “ - + n. owes 
Bye SOLE ere! oe 50 ee. ai Msg Sane aS 
ey, 220! + Am, ee ae ome = me Loe 
r e 
+ £ : CE =-7 = Dain 
son Me “5.g1 5 b,v,e o : ! : 


= condensing vapor specific heat based upon the 
arithmetic average of row bulk inlet temperature 
and bulk exit saturation temperature, Btu/lbm°R 
i = enthalpy of the saturated condensing vapor at 
—,sat,% 
the row exit free stream partial pressure, 
Btu/1lbm 
oO . = specific heat of the non-condensable gas jj, 


PrGrJ 
based on row bulk inlet temperature, Btu/lbm°R 


a2 





i = saturated vapor enthalpy of the bulk vapor at row 


.@ | 
exit conditions, Btu/lbm 
To = the reference temperature for condensing vapor 
enthalpy, °R 
eat = saturation temperature of the condensing vapor 


at the row exit partial pressure, °R 
The partial pressure of the condensing vapor at the exit 
Seeene row 1s obtained from the molal flow rate of the conden- 
Sing vapor and the total molal flow rate, knowing the inlet 
molal flow rates and the molal condensation rate. Thus from 
the exit partial pressure of the condensing vapor and equation 
D-17, the row exit temperature for superheated conditions may 


be obtained. 


a AE a cle ya a 
is,V,2 . oe - v,e v,e 
m =—EAM oe 2 Me! eee : 
[ ( Vv, ec Se p,v j=1 QrJ bg,3! 


¥ (mi -tAm_) (i ) 


m ae ek = 
ma «694d Pr9r) b,v,e Vi; 


UMAR 


g,sat,2 “p,visat 


-Am tli,te ( } (D-18) 


ee een eo 


GE 
In both equations D-17 and D-18 it is noted that the 
average non-condensable specific heat ug based upon the temper- 
ature of the bulk gas-vapor mixture entering the tube row and 
that the sensible heat transferred to the condensate film is 

based upon this same temperature. Actually, the tube row 
experiences some average bulk temperature between inlet and 
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outlet conditions for the transfer of sensible heat. However, 
Since the difference in bulk inlet and outlet temperatures is 
seldom much more than 10 degrees R,the effect of assuming the 
bulk inlet temperature for sensible heat transfer and bulk gas- 
vapor specific heat is negligible. 

For the saturated and wet vapor conditions, of course, 
bulk gas-vapor exit temperature is assumed at all times to be 
Boaquely determined by the partial pressure of the vapor 
leaving the tube row. A further assumption is necessary in 
the wet vapor case, and this involves the deposition of wet 
vapor droplets on condensing surfaces. If the droplets are 
sufficiently large, they will be deposited high in the conden- 
ser, i.e. in the first few rows. However, if they are 
sufficiently small, they may not be deposited until the latter 
few rows of the condenser are reached, and they may even pass 
completely through. Dalton's law of partial pressures, of 
course, can only be applied to the vapor portion of the two- 
phase mixture and even then it only approximates the actual 
Situation in a two-phase flow. In the development of the 
computer program for this thesis, it was assumed that no 
droplets are deposited on condenser tubes, which is optimistic 
With respect to condensation rates yet simpler to solve. The 
assumption was also made in the energy equation (for the wet 
vapor case) that equilibrium is maintained between droplet 
temperature and bulk vapor temperature. Equation D-16 may 


now be modified to include the vapor droplet enthalpy flux, 
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e 


m 1 a 
a Vk 


a a a ++ a ae as F one ene, 


we a 


e e 


ieee! eo ee 
1 Sled MSIge cs Vre@ vre d,2 d,k 


= 0 (Dai) 


I 
Le OI 


where 

m /m = the mass flow rate of the wet vapor droplets 

d,%° d,e 
leaving/fentering the tube row assuming none 
are deposited on the tubes of the tube row in 
question 

1 aah = the enthalpy of the wet vapor droplets 

a2 d,e 
leaving/fentering the tube row, assumed to 
exist in equilibrium with the bulk condensing 


vapor 


Knowing the condensation rate, one may then solve equation 
D-19 iteratively for the bulk outlet temperature and the 
corresponding wet vapor droplet flow rate leaving the tube row. 
This, of course, assumes that the possibility exists for some 
evaporation of the wet vapor droplets to maintain the equilibrium 
condition. This assumption admittedly is somewhat unrealistic 
Since equilibrium between droplets and the bulk gas vapor 


mixture will probably not exist. 
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Vavor-Non Condensable Gas Mixture Prope ries 


Tt was necessary to determine properties of mixtures of 


steam, 


a polar gas, and the three possible non-condensable 


gases, hydrogen, oxygen, and nitrogen. The mixture properties 


were determined as follows, using methods set forth by Reid 


and Sherwood (11), as modified or amplified below: 


1) 


Zz) 


3) 


4) 


5) 


Viscosity - by the method of Buddenburg and Wilke (12) 
Thermal conductivity - by the method of Lindsay and 
Bromley (13) as modified by Bennett and Vines (14) for 
polar-non-polar gas mixtures 

Density - by Amagat's Law (reference 11, pp. 321-322) 
Since no reliable method exists for determining 


densities of mixtures containing polar components: 


= le oy: ae where y is the molal density and y,. 
Y deny a 
m a i 
the mole fraction of component i 
Specific heat from the molal average value: Co aS 

id 
bey. CC. ., where © is the mollal specific heat, 
according to reference ll, pp. 323-325. Again no 
reliable method exists for determining the specific 
heat of polar mixtures. 
Diffusion coefficient, reference 11 (pp. 523-528) and 


Wilke (15) 


Computation of a mean diffusion coefficient for a mixture 


of steam and several non-condensable gases in a diffusion/ 


momentum/temperature boundary layer is an extremely complicated 


Miteration. lf it is assumed that the Lewis number and Prandtl 
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Mumber are both unity, the problem is tractable. This, of 
course, assumes that the rates of momentum, heat and mass trans- 
fer in the boundary layer are all equal. 

Wilke (15) has developed a method for determining the 
average diffusion coefficient for the diffusion of a single 
gas through a mixture of stagnant gases: Applying Wilke's 


method to a situation where there exists a temperature gradient, 











dy NeOpmeRn(c) Cy y y 

- Woe ee ov (5° ~ > + = j 
G P v-a v-b vc 

d nn 7 Ne Sn Rene) es 

e = 


d a N, OF: R P(e) Yy 
Do ao) 
v-b 


dy NM Oar Roe) we 


cee Vv OT 
dt ee (D=20) 
where 
etc) = the mole fraction of the diffusing (condensing) 
vapor 
wets). Yp,'o)s yee) = mole fractions of non-condensables 
ay Dye 
¢ = non-dimensional boundary layer thickness, zero at the 


free stream extremity and 1 at the condensate-vapor 


interface 
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Noo = molal flux of the ditiusing. (condensing) evanceu, 
Ib moles /hr £t- 


R = Universal gas constant 


E@talweressure, bert. 


1 Dopp: Plug = binary diffusion coefficients for the 


vr-a Vv 


condensing vapor and each of the non- 
condensables a, b, c, each functions 


of temperature, T(Z) 


T = temperature in the thermal boundary layer as a 
PInCELen Of Gp. 
6,, = thickness of-the temperature boundary layer, ft 


The multicomponent diffusion coefficient Oe is definedwas 


follows: 
a Ny Sep i Tmean 
‘ae oe a 
mp) on (= itm 
y,, 9) 
where Bean is the mean temperature in the boundary layer, °R. 


The solution of equation D-20 is based upon known boundary 
conditions at zt = 0 (the bulk conditions) and the assumed 
condensate vapor interface temperature, which fixes ae at C = as 
y,(1), Yp)/ (1), y,@) are all unknown as is the combination of 
terms ae Ome 

If the temperature distribution in the boundary layer is 
known, equation svstem D-20 may be solved numerically for the 


eigenvalue Ne K/p, with the condition at the condensate 


Sp 
mMapor interface, t = l: 
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y,, (1) 1 y, (1) + Yp (1) a YA) = 1 (= 22) 


For purposes of this thesis, it was assumed that Le=Pr=1. 
Hence a 1/7 power turbulent momentum boundary layer profile was 
assumed for the temperature profile in the turbulent region, 
and a parabolic profile in the laminar region. In the region 
of boundary layer detachment, of course, neither of these pro- 
files represents the existing situation. The laminar parabolic 
profile, strictly applicable only in plane Poisuille flow, has 
been assumed because of its simplicity. Admittedly this is a 
somewhat crude ae ee ination to the existing Situation, which 
probably most closely resembles the developing boundary layer 
maea £lat plate, at least prior to the point of separation. To 
use the Blasius profile, dependent upon the distance downstream 
of the stagnation point, would excessively complicate the 
solution to equation system D-20. To be sure, any assumption 
of a general boundary layer velocity and temperature profile 
for flow normal to a round tube is a crude approximation at 
best, and any reasonable one is almost as good as another. Any 
Significant refinements to a general temperature/velocity pro- 
file will require detailed understanding of the angular 
dependence of mass/momentum/heat transfer for condensing flow 


normal to a round tube. 


Fluid EEOpereles 


Table D-l gives a summary of references for the properties 


of fluids utilized in the condenser analysis. In each case, 
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Property 
Density 


Viscosity 


Specific Heat 


Thermal 
Bonductivity 


TABLE Da 


References for Fluid Properercs 


ioe. 


Liquid Water 
Salt Water 
steam 
Hydrogen 
Oxygen 
Nitrogen 
Liquid Water 
Salt Water 


Steam 


Hydrogen 
Oxygen 
Nitrogen 
Liguid Water 
Salt Water 
Steam 
Hydrogen 
Oxygen 
Nitrogen 


Liquid Water 
Salt Water 


steam 
Hydrogen 
Oxygen 
Nitrogen 


Reference 


iG; 2g 22 
24, 26, 27 
BG, 17 

ES 

18 

18 

22 jelous 
23526 


20 


20 
20 
20 
10s) 
26 
15S) 
Ue) 
i) 
Lys) 
20 


26 


AES) 
19 
19 
19 
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Comments 


pressure dependence 
from fresh water 


data, 


reference 25 


extrapolated below 


Zoe 


pressure dependence 
from reference 21 


except for steam, a least squares curve fit of the best 


available data was performed. A computer curve fit for steam 


table data, based on Keenan and Keyes data (16) was obtained 


from the Naval Ship Engineering Center (17). 


Description of the Computer Program 


A computer program in FORTRAN IV was developed for the IBM 


System 370 to estimate condenser surface area, weight and 


volume for an assumed condenser geometry, inlet steam and non- 


condensable conditions and sea water inlet temperature and 


eepth. 


a) 


b) 


e) 


dq) 


e) 


£) 


g) 


The detailed input data are listed below: 

mass flow rates of steam and non-condensable gases 
(hydrogen, oxygen, nitrogen) 

bulk temperature and total pressure of the entering 
gas-vapor mixture 

enthalvy of the vapor portion and quality of the 
entering steam 

sea water depth, inlet temperature, and velocity 
through the condenser tubes 

tube thermal conductivity, active length, inside and 
outside diameters, outside area per unit length, 
internal flow cross sectional area, and number of tubes 
per row 

condenser shell (steel) thickness 

header (waterbox) density (assumed identical wane tube 
density), yield stress and safety factor (cylindrical 


tube headers) 
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h) 


Results 


a) 


b) 


e) 


dq) 


e) 


an initial estimate of the overall coefficient of 
heat transfer, Uo: 

of the program are listed below: 

row by row condensation rate, fluid inlet and exit 
temperature and enthalpy data, if desired 

condenser exit fluid conditions, including mass flow 
rate and temperature of uncondensed steam, and average 
condensate temperature and enthalpy 

total number of tubes 

surface area 

condenser dimensions, volume and weight for deep 


submergence applications, if desired 


The condenser exit conditions were arbitrarily set at that 


point where the condensation rate for a particular row became 


less than one thousandth of the entering steam mass flow rate. 


The condenser program was tested under a variety of 


conditions as well as against examples in references 4 and 28. 


some of the test conditions are given in Table D-2. This table 


Shows the increase in surface area with increasing amounts of 


non-condensables under saturated and superheated steam condi- 


ELOnsS. 
rch an 


fom /hr , 


While the surface area does not increase drastically 


increase of oxygen flow rate from 0.605 lbm/hr to 8.0 


it should be noted that there is still a significant 


amount of uncondensed steam in the 8.0 lbm/hr case. This is 


indicative of the cutoff method used in the computer program, 


i.e. computation ceases when the condensation rate for a row 


in question is less than some fraction of the inlet steam mass 


SZ 


Entering) Temperature, °F Petce gomelsee?s.) 9138.75 $138.75 i38.75mmeIS.75 $138.75 138.75 38.75 138-75 138./5 138-975 138.75 leu reece 
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Demet ebtuy/ tomer ii2l.4 1Ti21.4 1121.4 1121.4 1121.4 9121.4 1121.4 Hi21.4 It2i.4 UN2ZT.a Vi2t.4 ll2l.4 1121.4 Vio ee care eer 
My | oualit. : , q 5 i - ‘ : 5 s 5 ; 21 4 ce ia 1211.4 i 
Non-Condensable | H O- 0 0 0 0 0 0 0.0754 0.125 0.1875 0.250 0.315 9.375 1.00 0 ; er me ey 
. O, +605 io nyo 2=0 oe 3.0 8.0 0 0 0 0 0 0 0 0.605 an oo 
Mass Flow Rate, lbm/hr) N 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
Row Number Condensation Rate -~ lbm/hr 
1 41.1 4 ieee 40.8 40.9 40.5 56.5 39.4 41.0 Ad: 40.8 40.6 49.4 40.3 39 34 SiGiecl: 34.9 SS 3 4s 
2 40.9 40.9 40.8 40.8 40.7 40.6 Bio 6.5 41.0 41.0 40.9 40.6 a0<5 NO) 22 39.49 36.4 es pel 3660 34.2 
5 Sy re ce ee BO ed Bs 0 36.9 56.9 Says 2) Be ioe So Jk 2 6) pious BiG ec 36.6 3 5 ae 33 4A 29 30 ee 2. Se 7 
4 Bee Balak 518) Aas, 56)./6 3652 36.4 Boies a a | SB ee S119) 5 2, 36.8 SOS 36.4 34.9 aa 5 323g) 34.0 J 2a 
5 346) 405 Chere. Sree 34.8 34.6 34.4 B26 Boe 3526 3573 35.0 34.8 34.6 32.8 55 ...8 Sl A 33.4 0, 
6 553.10 34.6 34.0 33.5 33/0 B2.7 30.3 35.2 34.8 Hal, © 33.9 33.4 Bo. 0 30.4 le: 30.4 34.0 2002 
7 30.6 29.0 28.4 272.7 1 27 26.8 24.5 320 30.9 BOel 29.4 26.7 28.2 25.6 B29 yg es 330 28.0 
g 12.9 14.0 14.5 MS i2 150 15.9 7 ao TV 7 12.6 13.6 Lab 15.1 15.7 417.9 24.2 21.7 27.2 23.0 
9 Loe 22 S604 4.24 5.08 5.59 9.26 0398 Lo Zi Se) 4.13 Ao 9,9 5.33 138 Be 14.3 
10 Oa 0c Om Za. 0.51 Oa3 5 1. 2g #.59 4.20 O92 OA ate C257 Orsi 1 0G 3 .9C Onoda Gavel j-28 76s 
iid: 0.08 0.14 0. 2a 2.235 e728 0.08 0) IES) 0.22 1.56 Nee 22 01 DIRS) 
12 ».10 Oe ial 0:.62 1.26 Le 
13 0.32 0.24 0.50 0.45 
14 Ons Oreck 
a. Eas, lpm/he 997 Opies 0.16 0.24 0.34 ).34 0.91 Oe A 0.20 0.34 0.42 M55 W387, 1.74 0.06 0.98 6.13 ths BE 
Surface Area, ft* ile ee eet 14.4 14.4 14.4 15.7 18.3 11.8 13.1 13.1 14.4 re oe ile) OL 14.4 Lees Loe Eye 
Cooling water temp/velocity, °F/ft/sec 75/7 
Total pressure - 2.8 Teste ey 
Buloeetlength/OD/iD. in 12.0/0.25/0.19 
ioe On) tLUbes per row - 20 
Tube spacing, in - 0.50 
Tube thermal Sonaucecivity, BSeuyhr tt °F — 16.5 
Steam rate, lb/hr - 272.0 TABLE D~-2 


Benlemeoetrelcient, Btu/hr ft* °R - 2000 Condenser Data 
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Memeerate (0.001). An interesting observation in the super- 
heated cases is that the second, Pee and sometimes sixth 
rows exhibit slightly higher condensation rates than the first, 
third, and fifth rows, respectively. An explanation for this 
phenonmenon is the combination of reduction of sensible heat 
transferred in succeeding rows, tending to increase condensa- 
tion rates, and an increase of condensate film resistance, 
tending to reduce condensation rates. Of course, when the 
thickening of the condensate film becomes significant (after 
the first few rows) and as superheat is decreased, this effect 
diminishes. 

A test of the program against an example given in Kern's 
Meecess Heat Transfer (4) resulted in excellent correlation in 
the first few rows, but correlation was not possible in latter 
rows since Kern neglected condensate film resistance. This 
comparison tested the program for cases where large amounts of 
non-condensables are present. A test was also made for very 
Small amounts of non-condensables (.05% by weight) by comparison 
with an example given in the U. S. Navy's condenser design data 
Sheet (28). The result was very conservative, predicting a 
condenser area approximately 150% of that in the Poe data 
Sheet example. This significant discrepancy can be explained 
by the fact that condensate drip, if allowed to remain on the 
tubes in a very large condenser (steam rate 247,000 lbm/hr), 
results in a Serious degradation of heat transfer in the lower 


parts of the condenser. For the relatively small condensers 
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Bensidered in this thesis, baffles to remove condensate drip 
will not be so important. 

A considerable amount of difficulty was experienced in 
consistently achieving convergence in the iterative solution 
scheme utilized. The scheme used throughout the program is 
first a bracketing method which then "homes in" on the solution 
by repetitively bisecting the interval in which the solution is 
known to lie. Iteration of the following variables was required 
mer a SOlution: 

1) Outside tube wall temperature 

2) Vapor-condensate interface temperature 

3) Bulk coolant temperature 

4) Condensation rate 

The number of required interface temperature iterations 
was decreased by the assumption of an initial overall coefficient 
of heat transfer, Uo and this was updated for succeeding rows 
after the first row, using that of the previous row. Condensa- 
tion rate proved to be the most troublesome, particularly in 
the latter few rows. It was found here that computation of 
condensation rate based upon row inlet conditions invariably 
resulted in predicting condensation rates in excess of the 
entering steam mass flow rate, which is obviously impossible. 

It was concluded that the Chilton-Colburn analogy, if based upon 
tube bank data, is not ees eanie to a Single row in the bank 
unless some average value of the inlet and exit steam flow rates 
from the tube row is used for computing the condensation rate. 


Iteration using the arithmetic mean of inlet and exit steam 
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mass tlow rates was also unsuccessful. Iteration was successful, 
however, when condensation rate was computed using a steam mass 
flow rate nearly that of the exit of the tube row. While this 
may, at first glance, seem no more than a trial and error scheme 
of getting just any solution, the dependence of condensation 
rates on row exit conditions is defensible. One must remember 
that the Chilton-Colburn analogy is only accurate when no 
Significant reduction in mass flow rate occurs across the mass 
transfer element in question. The mass transfer analogy is 
based on heat transfer data, in which there is no reduction in 
Mass flow rate. Hence, one may not expect the analogy to give 
accurate results when reductions in mass flow rate as great as 
50% or more occur across a given row. The Chilton-Colburn 
analogy may be used, however, if some average flow rate of the 
condensing vapor, probably one near the exit conditions, is 
assumed for computing mass transfer. In any event, computation 
of the condensation rate for a given row based upon a steam 

Mass flow rate nearly that of the row exit conditions is con- 


servative. 


Condenser Program Input Form 


The following data must be supplied to the condenser 
program in the FORTRAN format given below: 
Sara A - FORMAT 315 
KDATA - Number of condensers to be analyzed 


NOPT - Option for deep submergence weights and volumes: 


NOPT 
NOPT 


0, weights and volumes not desired 


1, weights and volumes desired 
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KOPT - Ovtion for row~by~row thermodvnamic and 


condensation rate data: KOPT = 0, do not desire 


row-by-row information; KOPT = 1, row-by-row 


information desired 


Card A is followed by 5 data cards for each condenser to 


be analyzed. 


Condenser data cards 


Caca 1 - 
MDOT ( 
MDOT ( 
MDOT ( 
MDOT ( 


Card 2 - 


FORMAT 4F10.4 


1) - mass flow rate 
2) - mass flow rate 
3) - mass flow rate 


4) - mass flow rate 


FORMAT 215, 6F 10.4 


of 
OF 
oe 


ona 


entering steam, lbm/hr 
hydrogen, lbm/hr 
oxygen, lbm/hr 


nitrogen, lbm/hr 


NROW - the initial tube row encountered, 1 


NTUBES - the number of tubes in a row 


oe rN — 


ee N = 


bulk temperature of the entering gas vapor 


mrexcure, °R 


bulk temperature of the entering sea water, °R 


HVENT - enthalpy of the entering steam vapor (use 


saturated vapor enthalpy for wet-vapor), Btu/lbm 


SPACE - the space between tube centers in a tube row 


(see figure D-5), 


ft 


LENGTH - tube active length, 


BOUT — 


Gata 3 - 


BIN = 


Eube>o.D., fe 
FORMAT 7F 10.4 


GUDe 2.0.7, she 


£t 


HSC - inside tube scale heat transfer coefficient, 


Beu/nr ft- °R 
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KTHWAL - tube: wall therraMrconductivitv, .5tU/nn fee f 
DEPTH - seawater depth, ft (use design depth if 
NOPT = 1 is desired) 
AXN - tube internal flow cross sectional area, ft? 
AOUT - tube external surface area per foot of length, 
sn (me 
acd 4 - FORMAT - 3F 10.4 
Bet = total condenser anlecl pressure (may be the 
average condenser pressure if inlet pressure 
is known and pressure drop can be estimated), 
1e7in- 
UOVRL - initial estimate of overall heat transfer 
coefficient from vapor-condensate interface 
to bulk coolant, based upon outside tube area 
(500-900), Btu/hr ft” °R 
QUAL - entering steam quality, 0 to 1.0 
Card 5 - FORMAT - 4F 10.6 
SHLTHK - thickness of the condenser steel shell, in 
DENS - density of the tube/header material, lbm/ft? 
SF - safety factor for the tube headers (1.5, 2.0 etc.) 
Wo —- yield stress of the tube header material (lbf/in“) 
The program requires approximately 80,000 bytes of storage. 
128,000 bytes of storage are sufficient when the FORTRAN G 
compiler is included. The IBM Scientific Subroutine Package 
Subprogram HPGC (Hamming's modified predictor corrector method 
for the solution of a simultaneous set of linear first order 


differential equations) must be available for operation of the 
condenser analysis program. 
Soa 





fg 


nan OW 


LIST OF SYMBOLS 


outside area of a condenser tube per foot of tube 
lengthyeeeea/ ft 

inside area of a condenser tube per foot of tube 
Henge, £t-/ft 

specific heat, Btu/lbm °R 

average condensing vapor specific heat, Btu/lbm °R 
molal specific heat, Btu/lb mole °R 

average condensate film specific heat based upon 
Nusselt's average film temperature, Btu/lbm °R 
outside condenser tube diameter, ft 

inside condenser tube diameter, ft 

G@eefusion coefficient (mass diffusivity), £t*/hr 
mass flux, lbm/hr ft? 

@eavitational constant, 4.17 x 10° f£t/hr- 

heat transfer re ee Beu/hre £t° OR 

enthalpy, Btu/lbm 

heat of vaporization of the condensing vapor at the 
partial pressure at the vapor-condensate interface, 
Btu/lbm 

heat rejected to the bulk coolant by the vapor 
condensed, Btu/lbm 

Colburn jJ-facter 

gas film mass transfer coefficient lb-mole/hr-ft?-lbf/in? 
ehermal conductivaty, Btu/hr ft OR 


effective tube length, ft 


Boe 


AT EM 


molecular weight, lbm/lb-mole 

mass flow rate lbm/hr 

condensation rate per tube in a tube row, lbm/hr 

molal velocity flux, lb moles/hr ft? 

pressure, Ibf/in- or libre" 

gas film log mean partial pressure of the non-condensing 
gases, Ibf/in- 

heat transfer rate, Btu/hr 

sensible heat per condenser tube transferred through 
the gas film, Btu/hr 

Universal gas constant, 1545.43 £t Ibf/°R Ib mole 
total number of distinct non-condensable gases present 
temperature, °R 

log mean temperature difference between the condensate- 
vapor interface and the bulk coolant, °R 

number of tubes in a row 

overall heat transfer coefficient based upon outside 
tube area and condensate-vapor interface and bulk 
vapor temperatures, Btu/hr ft* °R 

mole fraction 

molal density, lb-moles/ft? 


Wascosity, lbm7hr ft 
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sc 


refers 
refers 
refers 
refers 
refers 
refers 
refers 
refers 
refers 
refers 
refers 
refers 


rorers 


Le) 


16) 


to 


to 


to 


to 


ia @ 


to 


to 


Lo 


EO 


to 


EO 


SUBSCRIPTS 


bulk fluid corditions 


the 


condensate film or condensation rate 


droplets in a wet vapor condition 


conditions entering a tube row 


the 
the 
gas 
the 
the 


log 


saturated liquid condition 
saturated vapor condition 
film 

vavor condensate interface 
coolant 


mean conditions 


mean conditions 


a particular row 


outside conditions (tube wall) or reference 


conditions for enthalpy 


refers 
refers 
tube 

refers 
refers 


refers 


to 


tO 


EO 


to 


to 


constant pressure conditions 


scale on the coolant side of the condenser 


the condensing vapor 
the gas vapor mixture 
tube wall 
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